Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,558 T→A 7.9% intergenic (+24/‑172) 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,5580TA 7.9% 366.8 / 21.7 145intergenic (+24/‑172)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (44/83);  new base A (2/9);  total (46/99)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 151 total aligned reads displayed.

AGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTT  >  NC_000913/164409‑164702
                                                                                                                                                     |                                                                                                                                                
aGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCtgt                                                                                                                                                  >  1:1089149/1‑150 (MQ=255)
aGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCtgt                                                                                                                                                  <  2:282388/150‑1 (MQ=255)
aGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCtgt                                                                                                                                                  >  1:1089084/1‑150 (MQ=255)
 gCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACAGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCtgtg                                                                                                                                                 >  1:3760663/1‑150 (MQ=255)
 gCGACTTCTGGAAAGGAGAGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCtgtg                                                                                                                                                 <  1:345554/150‑1 (MQ=255)
  cGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTgg                                                                                                                                                >  2:2538153/1‑150 (MQ=255)
      ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAACGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                            <  2:3682190/150‑1 (MQ=255)
      ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                            <  1:2747324/150‑1 (MQ=255)
      ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                            <  1:598184/150‑1 (MQ=255)
      ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                            <  2:500611/150‑1 (MQ=255)
      ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                            <  1:1948390/150‑1 (MQ=255)
      ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                            >  2:1151366/1‑150 (MQ=255)
      ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                            <  1:942712/150‑1 (MQ=255)
      ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGGTATCCCAAACATGTCTGGCCGGACGACCCGGCACATACTGCACCGACGCGACGGCCGAACCAGGCTTCGTCCGGTTGGTTTTCTCCCGCTCCCTGTGGgagc                                                                                                                                            >  1:4337836/1‑149 (MQ=255)
      gtCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                            <  1:3949442/149‑1 (MQ=255)
       tCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGagag                                                                                                                                           <  2:2976247/150‑1 (MQ=255)
        cTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAACCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGAGAgg                                                                                                                                          <  1:161102/150‑1 (MQ=255)
        cTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAgg                                                                                                                                          <  1:4095175/150‑1 (MQ=255)
         tGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAggg                                                                                                                                         >  1:528834/1‑150 (MQ=255)
           gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGAGAgggag                                                                                                                                       <  2:2675858/150‑3 (MQ=255)
           gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGAGAgggag                                                                                                                                       <  2:2675950/150‑3 (MQ=255)
           gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCCCCCTCTCCCACAGGGAGAgggag                                                                                                                                       <  2:2805720/150‑3 (MQ=255)
           gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCCCCCTCTCCCACAGGGAGAgggag                                                                                                                                       <  2:4953961/150‑3 (MQ=255)
              aGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                    <  2:264550/150‑1 (MQ=255)
              aGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                    <  2:264501/150‑1 (MQ=255)
                gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                                                                                                                                  >  2:3255889/1‑150 (MQ=255)
                gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                                                                                                                                  >  1:4970512/1‑150 (MQ=255)
                   cGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                               >  2:2668468/1‑150 (MQ=255)
                    gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                              <  2:1954640/150‑1 (MQ=255)
                     tACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAggg                                                                                                                             <  2:1348900/150‑1 (MQ=255)
                     tACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAggg                                                                                                                             >  2:4499318/1‑150 (MQ=255)
                       ccGTGAGGGGCAACCAGAGATGAAAGGGCGTTATCCCACACATGTCGGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                           <  2:121098/150‑1 (MQ=39)
                        cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                          >  1:451544/1‑150 (MQ=255)
                         gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                         >  1:480864/1‑150 (MQ=255)
                         gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTACTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                         >  2:3139414/1‑150 (MQ=255)
                                 cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGTCGGGGTGAGGGCATCAGCGCGCa                                                                                                                 >  2:4836622/1‑150 (MQ=255)
                                  aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCAc                                                                                                                >  2:3332189/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)
                                              aaGGGCGTTATCCCAAACATGTCGGGCGGGACGACCCGGAAATTACTGCACCGACGCGAGGGACTAAAAAGTTTTCGTTTGTTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                    <  2:1174350/150‑1 (MQ=25)
                                                gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctcaccc                                                                                                  <  2:1187937/150‑1 (MQ=255)
                                                gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGAACCGACGCGACGGAGGAAAATGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcaccc                                                                                                  <  2:1187996/150‑1 (MQ=38)
                                                 ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccacacccc                                                                                                 >  2:3030252/1‑149 (MQ=255)
                                                    gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                              <  2:1975511/150‑1 (MQ=255)
                                                    gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGCCCCACACCCAAAc                                                                                              >  2:3993118/1‑150 (MQ=38)
                                                    gTTAGCCCAACCAGGTCGGGCCGGACGACCCGGCAAATCCTGCACCGACGCGCCGGACGAAAAAGTATTCGTAAGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                              <  2:711507/150‑1 (MQ=37)
                                                     ttttCCCAAACTTTTGTGGCCGGACGACCCGGCAAAGACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                             <  2:834657/147‑1 (MQ=38)
                                                     ttttCCCAAACTTTTCTGGCCGGACGACCCGGAAAATACGGCACCGCCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCGTGTGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                             <  2:834817/147‑1 (MQ=38)
                                                     ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCAGGTTCACCCGCACCCGAAcc                                                                                             >  2:3269308/1‑150 (MQ=39)
                                                         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)
                                                                aggTCGGGCCGGACGACCCGGCAAATCCTCCCCCGACGCGCCGGACTAAACCGTATTCGTATGGTGGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                  <  2:3300421/148‑1 (MQ=25)
                                                                  gTCGGGCCGGACCACCCGCCAAACACGCCCCCGCCGCGCCGGCCGAAAACGTTTTCTTATGGTTTGTTTTCCCCCTCGCCGTGTGGGAGAGGGCCGGGGTTTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                <  2:3034367/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)
                                                                            acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCCGACCCGCGCCCCCACGGGGGAgggggg                                                                      >  1:989918/1‑148 (MQ=25)
                                                                              gacCCGCAAAACACTGCACCGCCGCGCCGGAAGAAAATGTTTTCGTATGGTTGTTTTCCTCCCTCTCCGTGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                    <  1:3214391/150‑1 (MQ=255)
                                                                                 ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCCAACCCACTCCCGCAAGGGCGGGGGGGcccttc                                                                 >  2:1733204/1‑146 (MQ=37)
                                                                                 ccGCCAAATACGGCACCGACGCGCCGGACTAAAAAGTATTCGTAAGGTTGGTTTCCTCCCACTCCCTGTGGGAGAGGGCGGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                 <  2:122378/150‑1 (MQ=255)
                                                                                  cGGCAAATACTGCACCGACGCGACGGACTAAAATTTTTTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                <  2:4501213/150‑1 (MQ=38)
                                                                                    gtaaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCACACCCAAACCCACACCCACCACGGGGAGGGGGCCCATCGac                                                              >  2:1907412/3‑149 (MQ=25)
                                                                                    gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCGCACCCGAACCCACACCCGCAAGGGCGGGGGGGcccagagag                                                              >  2:2938907/1‑143 (MQ=25)
                                                                                      aaaaACTCCACCCACCCCACGCACTAAAAACAATACTTACGTTTGGTTTACACCCACCCCCTGGGGGTGAGGGCCGGGGTTAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                            <  2:908887/150‑1 (MQ=255)
                                                                                          aCTGCACCTACGCGACGGACTAAAAAGTATTCGTAAGTTTGGTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACaaa                                                        <  2:3760663/150‑1 (MQ=255)
                                                                                                cgtacgcgactGACTAAAATTTATTCTTATGGTTGTTTTTCTCCCTCTCCCTTTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                  <  1:796438/147‑1 (MQ=255)
                                                                                                ccgccgCGCCGGCCGAAACTGTATGCGTATGGTGGGATTCCGCCCACTCCCTGTGGGCGAGGGCCGGGGTGTGCGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGTGGACCGATCGAGCACAAATTTTga                                                  <  2:3091221/150‑1 (MQ=255)
                                                                                                ccGACGCGACGGACGAAAATGTATTGGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                  <  2:2591858/150‑1 (MQ=255)
                                                                                                  gacgcgacgGACTAAAATTTTTTCTTATGTTTGGTTTTCCCCCTCACCGTGGGGTTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                <  1:776862/150‑1 (MQ=255)
                                                                                                  gacgcgacgGACGACCACGTATTCGTATGGTTGGTTTTCCCCCACTCCCTGTGGGTGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                <  2:2712852/150‑1 (MQ=255)
                                                                                                  gacgcgacgGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                <  2:2713088/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)
                                                                                                             cgAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                     <  2:2669266/150‑1 (MQ=255)
                                                                                                              gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCACACCCACAAAGGAGAGGGGACCGAACGAGCACAAAATTTGATAGAAAtattaatc                                    >  2:3410738/1‑146 (MQ=25)
                                                                                                               aaaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                   <  2:5104183/150‑1 (MQ=255)
                                                                                                                    gTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                              <  1:3824204/150‑1 (MQ=255)
                                                                                                                    gTATTCGTAAGGTGGGTTTTCCCCCTCGCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                              <  2:4199474/150‑1 (MQ=255)
                                                                                                                          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)
                                                                                                                                  gtttttCCCCCTCCCCGGGTGGGTGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                <  1:3935136/148‑1 (MQ=255)
                                                                                                                                      ccggccctctcccaCAGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCCAACCCGCTCCCGCAAGGGAGTGGGGGCCGATCGAGCCCAAAATTTGGGAGATATATTATTATGTCTTGTAACAGGCGAACAGaaa            >  2:823657/5‑150 (MQ=255)
                                                                                                                                          ccctctccctGGGGGCGTGGGCCGGGGTGAGGGCATCCGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAGCg        <  2:4555377/150‑1 (MQ=255)
                                                                                                                                          cccacgccctGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg        <  2:1242200/144‑1 (MQ=255)
                                                                                                                                                ccctGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt  <  1:1854698/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
AGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTT  >  NC_000913/164409‑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: