Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,557 G→C 7.6% intergenic (+23/‑173) hrpB → / → mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5570GC 7.6% 478.9 / 29.0 145intergenic (+23/‑173)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base G (44/88);  new base C (2/9);  total (46/99)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 152 total aligned reads displayed.

GAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTT  >  NC_000913/164408‑164702
                                                                                                                                                     |                                                                                                                                                 
 aGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCtgt                                                                                                                                                  >  1:1089149/1‑150 (MQ=255)
 aGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCtgt                                                                                                                                                  >  1:1089084/1‑150 (MQ=255)
  gCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACAGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCtgtg                                                                                                                                                 >  1:3760663/1‑150 (MQ=255)
  gCGACTTCTGGAAAGGAGAGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCtgtg                                                                                                                                                 <  1:345554/150‑1 (MQ=255)
       ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                            <  1:1948390/150‑1 (MQ=255)
       ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                            <  1:942712/150‑1 (MQ=255)
       ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                            >  2:1151366/1‑150 (MQ=255)
       ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                            <  2:2290270/150‑1 (MQ=255)
       ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                            <  2:500611/150‑1 (MQ=255)
       ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                            <  1:2747324/150‑1 (MQ=255)
       ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGGTATCCCAAACATGTCTGGCCGGACGACCCGGCACATACTGCACCGACGCGACGGCCGAACCAGGCTTCGTCCGGTTGGTTTTCTCCCGCTCCCTGTGGgagc                                                                                                                                            >  1:4337836/1‑149 (MQ=255)
       gtCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                            <  1:3949442/149‑1 (MQ=255)
        tCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGagag                                                                                                                                           <  2:2976259/150‑1 (MQ=255)
         cTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAACCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGAGAgg                                                                                                                                          <  1:161102/150‑1 (MQ=255)
         cTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAgg                                                                                                                                          <  1:4095175/150‑1 (MQ=255)
          tGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAggg                                                                                                                                         >  2:3935136/1‑150 (MQ=255)
          tGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAggg                                                                                                                                         >  1:528834/1‑150 (MQ=255)
           ggAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGAGAggga                                                                                                                                        <  1:4628956/150‑2 (MQ=255)
            gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGAGAgggag                                                                                                                                       <  2:2675858/150‑3 (MQ=255)
            gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGAGAgggag                                                                                                                                       <  2:2675715/150‑3 (MQ=255)
            gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCCCCCTCTCCCACAGGGAGAgggag                                                                                                                                       <  2:2805720/150‑3 (MQ=255)
            gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCCCCCTCTCCCACAGGGAGAgggag                                                                                                                                       <  2:4953961/150‑3 (MQ=255)
              aaGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCGCCCTCTCCCACAGGGAGAgggagaa                                                                                                                                     <  1:625621/150‑5 (MQ=255)
               aGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                    <  2:264500/150‑1 (MQ=255)
               aGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                    <  2:264501/150‑1 (MQ=255)
                 gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                                                                                                                                  >  2:3255889/1‑150 (MQ=255)
                    cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                               >  2:2668451/1‑150 (MQ=255)
                    cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                               >  2:2668468/1‑150 (MQ=255)
                     gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                              <  2:1954640/150‑1 (MQ=255)
                     gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                              <  2:1954597/150‑1 (MQ=255)
                        ccGTGAGGGGCAACCAGAGATGAAAGGGCGTTATCCCACACATGTCGGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                           <  2:121098/150‑1 (MQ=39)
                         cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                          >  2:1200040/1‑150 (MQ=255)
                         cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                          >  2:3088161/1‑150 (MQ=255)
                           tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGc                                                                                                                                        >  2:4081029/1‑134 (MQ=255)
                                  cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGTCGGGGTGAGGGCATCAGCGCGCa                                                                                                                 >  2:4836622/1‑150 (MQ=255)
                                   aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                >  2:3332490/1‑150 (MQ=255)
                                     aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                              >  2:4038042/1‑150 (MQ=255)
                                         cgaTGAACGGGCGTTATCCCAACCATGTCTGGCCGGACGCCCCGGCAAATACGGCACCGACGCGACGGACGAAAAAGTTTGCTTAAGGTTGTTTTTCTCCCTCTCCCTGTGGTTGAGGGCCGGGTTGAGGGCATCAGCGCGCACGTtcac                                                                                                          <  2:4649209/149‑1 (MQ=255)
                                           aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccc                                                                                                        >  1:4637082/1‑150 (MQ=255)
                                             gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTtcaccctc                                                                                                      >  1:3843570/1‑150 (MQ=255)
                                             gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctc                                                                                                      >  1:4512819/1‑150 (MQ=255)
                                             gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacccac                                                                                                      >  1:4283554/1‑148 (MQ=255)
                                               aaGGGCGTTACCCCAACCAGGCCGGGCCGGACGACCCGGCAAATCCTGCCCCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                    <  2:1174453/150‑1 (MQ=37)
                                                 gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcaccc                                                                                                  <  2:1187937/150‑1 (MQ=255)
                                                  ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccacacccc                                                                                                 >  2:3030252/1‑149 (MQ=255)
                                                    cGTTATCCCAACCATGTCTGGCCGGACGACCCGGCAAACACTGCACCGACGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                               <  2:3288066/150‑1 (MQ=255)
                                                     gTTAGCCCAACCAGGTCGGGCCGGACGACCCGGCAAATCCTGCACCGACGCGCCGGACGAAAAAGTATTCGTAAGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                              <  2:711507/150‑1 (MQ=37)
                                                      ttttCCCAAACTTTTGTGGCCGGACGACCCGGCAAAGACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                             <  2:834657/147‑1 (MQ=38)
                                                      ttttCCCAAACTTTTCTGGCCGGACGACCCGGAAAATACGGCACCGCCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCGTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                             <  2:834817/147‑1 (MQ=38)
                                                      ttATCCCACAGCTGTGGGGCCGGACGCCCCGGCAAATACTGCACCGACGCGCCGGAGGAAAATGTTTGGGTATGGTTGGTTTTCGCCCTCGCCGGGGGGGTGAGGGCGGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                             <  2:4534264/150‑1 (MQ=255)
                                                      ttATCCCAAACATGTCGGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                             <  2:1237082/150‑1 (MQ=255)
                                                          cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATGCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCGCCCGCACCCAAGcccaca                                                                                         >  1:3248782/1‑147 (MQ=37)
                                                           ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGGCCCGCACCCCCAcccacac                                                                                        >  1:351196/1‑146 (MQ=37)
                                                           ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCACCCGCCCCCTAAaccaaaa                                                                                        >  1:2030030/1‑146 (MQ=38)
                                                           ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCAACCAAAcccacac                                                                                        >  1:3695009/1‑146 (MQ=38)
                                                             aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCGCACCCCAACCCCCGccc                                                                                      >  2:4890369/1‑150 (MQ=38)
                                                              aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGGGCCCCCCCCCCCCAACCCCCAcccc                                                                                     >  1:4156677/1‑149 (MQ=37)
                                                               aCATGTCTGGCCGGACGCCCCGGCAACTACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                                                    <  2:830440/150‑1 (MQ=255)
                                                                 atttgtGGCCGGACGACCCGGAAAATACTGCACCTACGCGACGTACTAAAATTTATTCTTATGGTGGGTTTTCTCCCTCCCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                  <  2:1103105/145‑1 (MQ=255)
                                                                 aggTCGGGCCGGACGACCCGGCAAATCCTCCCCCGACGCGCCGGACTAAACCGTATTCGTATGGTGGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                  <  2:3300421/148‑1 (MQ=25)
                                                                 aTGTCTGGCGGTACGCCCCGGCAAATACTGCACCGAGGCGACGGACTAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                  <  2:1699713/150‑1 (MQ=38)
                                                                 aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCCCCCCAACCCACGcccgccc                                                                                  >  2:4446160/1‑146 (MQ=37)
                                                                 aTGTCGGGCCGGACGACCCGGCAAAAACTCCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                  <  2:2709110/150‑1 (MQ=255)
                                                                  tgacTGGCCGCACGACCCGGAAAATACTCCACCGACGCGACGCACCAAAAACTATACGTAAGGTTGTTTTTCTCCCTCACCCCGTGGGAGAGGCCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                 <  2:2545958/147‑1 (MQ=37)
                                                                  tGTCTGGCCGGACGACCCGGCAAATCCGGCACCCCCGCGCCGGACCAAAATTTTTTCGAACGGTTGGTTTTCCCCCTCCCCCGTTGGGAGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                 <  2:3506668/150‑1 (MQ=255)
                                                                         ccGGCCGACCCGGCAAATACTGCACCGACGCGTCGGACTAAAAAGTATTGTTATGGTGGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                          <  2:3625495/150‑1 (MQ=37)
                                                                         ccGCACGACCCGCCAAATCCGGCACCGACGCGACGGACGAAAATTTATTCTTACGGTTGGTTTGCCCCCCCCCCCGGTGGGAGTGGGCCGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                          <  2:2399012/150‑1 (MQ=255)
                                                                           ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGTGGCCCCACACCCGAACCCACCCCCCCACgggggggggg                                                                        >  2:4931599/1‑150 (MQ=25)
                                                                            gacgccCCGGCAAATACGGCACCGACGCGACGGACGACAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCCACAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                       <  2:3351156/150‑1 (MQ=39)
                                                                            gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                       <  2:5081377/150‑1 (MQ=255)
                                                                            gacgacCCGGCAAACACGCCCCCCCCGCGACGTACGAAAATGTATTCTTATGGTGGGTTTTCACCCCCTCCCTGTGGTTGCGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                       <  1:4932994/150‑1 (MQ=255)
                                                                               gacCCGGCAAATACTGCACCGACGCGACGTACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                    <  2:2906863/150‑1 (MQ=255)
                                                                               gacCCGCAAAACACTGCACCGCCGCGCCGGAAGAAAATGTTTTCGTATGGTTGTTTTCCTCCCTCTCCGTGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                    <  1:3214391/150‑1 (MQ=255)
                                                                                  ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCCAACCCACTCCCGCAAGGGCGGGGGGGcccttc                                                                 >  2:1733204/1‑146 (MQ=37)
                                                                                   cGGCAAATACTGCACCGACGCGACGGACTAAAATTTTTTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                <  2:4501213/150‑1 (MQ=38)
                                                                                     gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCGCACCCGAACCCACACCCGCAAGGGCGGGGGGGcccagagag                                                              >  2:2938907/1‑143 (MQ=25)
                                                                                       aaaaACTCCACCCACCCCACGCACTAAAAACAATACTTACGTTTGGTTTACACCCACCCCCTGGGGGTGAGGGCCGGGGTTAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                            <  2:908887/150‑1 (MQ=255)
                                                                                                 ccgccgCGCCGGCCGAAACTGTATGCGTATGGTGGGATTCCGCCCACTCCCTGTGGGCGAGGGCCGGGGTGTGCGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGTGGACCGATCGAGCACAAATTTTga                                                  <  2:3091221/150‑1 (MQ=255)
                                                                                                 ccgccgCGACGGACTAAAAAGTATTGGTACGGTTGTTTTTCTCCCTCTCCCGGGGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                  <  2:3685889/150‑1 (MQ=38)
                                                                                                 ccGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                  <  2:2734210/150‑1 (MQ=255)
                                                                                                   gacgcgacgGACTAAAATTTTTTCTTATGTTTGGTTTTCCCCCTCACCGTGGGGTTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                <  1:776862/150‑1 (MQ=255)
                                                                                                   gacgcgacgGACGACCACGTATTCGTATGGTTGGTTTTCCCCCACTCCCTGTGGGTGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                <  2:2712852/150‑1 (MQ=255)
                                                                                                   gacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGTGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                <  2:4250193/150‑1 (MQ=255)
                                                                                                         acgtacgAAAAATTATTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                          <  2:528834/150‑1 (MQ=38)
                                                                                                         acggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCAAACCCCAACCCACACCCACCAGGGCGACGGGGCCCATAGAGCACCAAATTTgtgatgtata                                          >  1:1647881/1‑144 (MQ=16)
                                                                                                            gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                       <  2:4700384/150‑1 (MQ=255)
                                                                                                              cgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCCAACCAACTCCCACAAGGGTGAGGGGAGCGATAGAGCCCAAATATTTGGAGATAtattaat                                     >  2:3232253/1‑147 (MQ=25)
                                                                                                                aaaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                   <  2:5104183/150‑1 (MQ=255)
                                                                                                                     gTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                              <  1:3824204/150‑1 (MQ=255)
                                                                                                                       tttCGTACGGTTGTTTTTCTCCCTCTCCGTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                            <  1:5139795/149‑1 (MQ=39)
                                                                                                                           gtttggttgtttTGCCCCCCCGCCGTGTGGGAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                        <  2:3897113/146‑1 (MQ=255)
                                                                                                                           gTATGGTTGGTTGCCTCCCTCGCCGTGTGGGTGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGACTTGTAACa                        <  1:2538153/150‑1 (MQ=255)
                                                                                                                              aGGTGGGTTTCCCCCCTCTCCCTGTGGGTGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                     <  2:4970512/150‑1 (MQ=255)
                                                                                                                               gtttggttTTCTCCCTCTCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                    <  1:2938907/148‑1 (MQ=255)
                                                                                                                                gttggttTTCTCCCTCTCCCTGTGGGAGAGGGTGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                   <  1:4836622/150‑1 (MQ=255)
                                                                                                                                   gtttttCCCCCTCCCCGGGTGGGTGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                <  1:3935136/148‑1 (MQ=255)
                                                                                                                                       ccggccctctcccaCAGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCCAACCCGCTCCCGCAAGGGAGTGGGGGCCGATCGAGCCCAAAATTTGGGAGATATATTATTATGTCTTGTAACAGGCGAACAGaaa            >  2:823657/5‑150 (MQ=255)
                                                                                                                                          tccctctccctGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc         <  1:929950/150‑1 (MQ=255)
                                                                                                                                           cccacgccctGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg        <  2:1242200/144‑1 (MQ=255)
                                                                                                                                                 ccctGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt  <  1:1854698/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTT  >  NC_000913/164408‑164702

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: