Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,595 T→G 11.6% intergenic (+61/‑135) 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,5950TG11.6% 606.7 / 23.6 217intergenic (+61/‑135)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (21/166);  new base G (27/2);  total (48/169)
Fisher's exact test for biased strand distribution p-value = 4.00e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 240 total aligned reads displayed.

GAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGAC  >  NC_000913/164448‑164744
                                                                                                                                                    |                                                                                                                                                     
  gaTGAATGGGCGTTATCCCAAACATGTCTGGCCGGCCGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTTTTGGTATGGTTGGTTTTCTCCCTCGCCCTGGGGGTGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcacc                                                                                                                                                   <  2:9591041/150‑1 (MQ=255)
   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccc                                                                                                                                                  >  2:8380001/1‑150 (MQ=255)
      aacgggCTTTATCCCAAACTTGTCGGGCCGCACGACCGGTAAATTACTGCACCGCGGCTAGTTATTAAAATTTTTTGGTATGGTTGTTTTTCTCCCGCTCCGGGGGGGAGTGGGCGGGGGTTTGGGCCTCAGCGCGCACGTtcaccctca                                                                                                                                               <  2:7901307/147‑1 (MQ=255)
       aaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGAGGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                                              <  2:7373564/150‑1 (MQ=255)
        aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTgcccccccccc                                                                                                                                             >  2:5954314/1‑145 (MQ=255)
        aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGGCCGGtcaacctcacc                                                                                                                                             >  1:200063/1‑150 (MQ=255)
         gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAAGTATTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcaccc                                                                                                                                            <  2:6434925/150‑1 (MQ=255)
          ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGgcacccacaccca                                                                                                                                           >  1:8204305/1‑149 (MQ=255)
            ggTTATCCAAAACATGTGTGGCCGGACGCCCCGGCAAATACTGCCCCTACGCGACGGACTAAAATGTTTTCTTATGGTTGTTTTTCTCCCTCTCCCTGTGGTTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                                         <  2:6929517/149‑1 (MQ=25)
                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTGGCCCacacccacacccaca                                                                                                                                   >  1:9688402/1‑147 (MQ=38)
                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAAcccaca                                                                                                                                   >  1:7561783/1‑147 (MQ=255)
                   ccAAACATGTCTGGCCGTACGACCCGGCAAATACTGCACCTAAGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTATCCCTCTCCCTGTGGGAGAGGGCCGCGGTGAGGGCATCAGCGCGCACGTTCACCAACACCCCAAcacacta                                                                                                                                  >  1:7883847/1‑144 (MQ=37)
                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCACAcccacac                                                                                                                                  >  2:9507566/1‑146 (MQ=38)
                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCCCCCCCGAACCCCCTcc                                                                                                                                 >  2:9015910/1‑150 (MQ=38)
                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCACACCCAAAACCACAccc                                                                                                                                >  2:5443308/1‑150 (MQ=37)
                         agggcTGGCCGGACGCCCCGGAAAATCCGGCACCGACGCGAGGGCCTAAAAATTTTTCGTAAGGTTGGTTTTCCCCCACTCCGGGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                            <  1:9029265/146‑1 (MQ=255)
                         aTGTCTGGCCGTACCCCCCGTAAAATACTCCACCGACGCGACGTACTAAAATTTTTTCGTATGGTTGTTTTCCTCCATCTCCCTGGGGGTGTGGGCGGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                            <  1:6548703/150‑1 (MQ=255)
                         aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCACTCCCACaa                                                                                                                            >  1:9155738/1‑150 (MQ=255)
                          tgtgtgGCCGGACGACCCGCAAATTCCTCCACCGACGCCAGTCACTAAAAATTATACGTAAGTTTGTTTTTCACCCTCACCCTTGGGTTGTGGGCCGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                           <  2:346953/150‑1 (MQ=255)
                          tGTCGGGCGGGACGACCCGGCAAATACGGCCCCGACGCGCCGGACTAAAAATTTTTCTTATGGTTGTTTTTCCCCCACACCCTGTGGTAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                           <  2:2073364/150‑1 (MQ=255)
                             cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGAGCCCGCACCCCAACCCACTCCCACAAcggg                                                                                                                        >  1:1206537/1‑149 (MQ=37)
                               gcccGGACGACCCGCCAAACACGGCCCCGCCGCGCCGGACGAAAAAGTATGCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGgaga                                                                                                                      <  2:5261615/148‑1 (MQ=37)
                               gcccGCACCACCCCGCAAATACAGCACCGACCCGACGCACCAAAAACTATTCGTAAGGTTGTTTTACACCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                                      <  2:9878060/148‑1 (MQ=37)
                                gCCGGACGCCCCGGAAAAACCGGCCCCGCCGCGACGGACGAAAATTTTTTCGTACGGTTGTTTTACCCCCTCCCCCGGGGGGAGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                                     <  2:3803005/150‑1 (MQ=255)
                                gCCGGACGACCCGGCAAATACCGCCCCGACGCGAGGGCCGAAAAAGTATTCGGATGGTTGGTTTTCGCCCTCCCCGGGGGGTAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                     <  2:8745086/150‑1 (MQ=255)
                                 ccgcccgacCCGCAAAATACTGCACCGACGCGACGGACCAAAAATTTTTCGTATGGTTGTTTTTCCCCCGCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                                    <  1:4182843/145‑1 (MQ=37)
                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTAACCCCCACCCAAACCCACACACGCAAGGGAGAggg                                                                                                                   >  2:3243707/1‑150 (MQ=37)
                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAAACCGCACCCACAAGGGGGAggg                                                                                                                   >  2:2452351/1‑150 (MQ=37)
                                  cGGACCACCCGCAAAAAACTCCACCGACGCCACGCAATAAAAATTATTCTAATGTTTGTTTTACACCCTCTCCCTGTGGTTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                                                   <  2:7938728/150‑1 (MQ=255)
                                      cgccccGGCAAATACGGCACCGCCGCGCCGGACGAAAAAGTATTCGTACGGTGGGTTTTCTCCCTCGCCCTGGGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                               <  2:1620714/147‑1 (MQ=255)
                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCCCCCAAACCCCCACCCACCAAGGGGAGGGGGCCgcgc                                                                                                           >  1:4725496/1‑147 (MQ=25)
                                          ccGGCAAACACTCCACCGACGCGACGCACGAAAAAGTATTCTTACGTTGGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                                           <  2:2893290/150‑1 (MQ=37)
                                             gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                                        <  2:6547352/150‑1 (MQ=255)
                                                 aTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCACCCAAACCCACACCCTCAAGGGAGAGGGGGCCCAGCGAGcacc                                                                                                    >  1:8375064/1‑149 (MQ=37)
                                                      gccccGACGCGCGTGACGAAAAATTATGCGTAAGGTTGTTTTTCTCCCTCTCCGTGGGGGAGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAttt                                                                                               <  1:4038195/147‑1 (MQ=255)
                                                       cACCGACGCGCCGGACTAAAAATTTTTCTTATGTTTGTTTTTCTCCCTCTCCGTGTGGTAGTGGGGCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                              <  1:7215339/150‑1 (MQ=255)
                                                        cccgccgGGACGGACTAAAATGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGGGGGTGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                                             <  1:7339851/149‑1 (MQ=255)
                                                        aCCGACGCGAGGGACTAAAATGTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGTAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                                             <  2:202093/150‑1 (MQ=38)
                                                         ccgcggcgccgGACGAAAAATTATTCGTATGGTGGGTTTTCCCCCGCCCCGTGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                            <  2:8896777/145‑1 (MQ=255)
                                                         ccgccgCGTCGGACGAAAAATTATTCGTATGGTTGTTTTCCTCCCACTCCGTGTGGGAGAGGGCCGGGGTTCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                            <  1:665824/150‑1 (MQ=255)
                                                         ccgccgCGACGTATTAAAATGTATTCTTCTGGTTGGTTTTCCCCCTCGCCGGGTGGGTTTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                            <  1:1087648/150‑1 (MQ=255)
                                                         ccGACGCGCCGGACTAAAATTTTTTCTTAAGGTTGGTTTTCTCCCACTCCCTGTGGTAGCGGGCCGGGTTGAGGGAATCAGCGCGCACGTTCACCCTCACCATAACCCTCTCCCACAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                            <  1:6988797/150‑1 (MQ=255)
                                                         ccGACGCGACGGACTAAAAAGTATTCGTATGGTTTTTTTTCCCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                            <  1:6184702/150‑1 (MQ=38)
                                                         ccGACGCCACGTACCAAAACGTTTTCTTACGGTTGTCTTTCTCCCTCCCCCTGTGGGAGTGGGCGGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                            <  2:4295399/150‑1 (MQ=255)
                                                           gacgcgccgGCCTAAAAATTATTCGTATGGTTGGTTTGCGCCCTCTCCCTGTGGGCGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                                          <  1:8473619/150‑1 (MQ=255)
                                                                gagggacgAAAAAGTTTTCGTATGTTTGGTTTACTCCCTCTCCCTGTGGTAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGGCCGATCGAGCACAAATTTTGAGAGatat                                                                                     <  1:4880652/147‑1 (MQ=255)
                                                                          caattATCCGTACGGTTGGTTTCCTCCCTCGCCCTGTGGGAGCGGGCCGGGGTGCGGTCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                                           <  2:1293290/146‑1 (MQ=255)
                                                                          aaaGTAGTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                                           <  1:8969281/150‑1 (MQ=255)
                                                                           atgttttCGTATGGTTGGTTTTCTCCCTCTCCCTGGGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                                          <  1:1568133/145‑1 (MQ=255)
                                                                             ttATTGGTATGGTTGTTTTTCTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                        <  2:5544567/149‑1 (MQ=39)
                                                                             gTATTCTCACGGTGGGTGTCCTCCCGCGCCCGGGGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                        <  1:5219806/150‑1 (MQ=255)
                                                                                 tCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                                    <  2:7407727/150‑1 (MQ=255)
                                                                                   gTACGGTTGGTTTTCTCCCTCGCCGTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                                  <  1:5020531/150‑1 (MQ=255)
                                                                                    tAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCCCCCCCACCCAAACCCACGCCCACCAGGGCGAGGGGGCCGATCGCGCCCAAAAGTTGGGTGGTATCTTCTTCTGTCTTgggaccg                                                                 >  2:2477422/1‑144 (MQ=16)
                                                                                     aaGGTTGGTTTCCGCCCTCTCCCTGTGGGAGAGGGCCGGGGTTCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                                <  2:731860/150‑1 (MQ=255)
                                                                                      tggttggttTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                               <  1:9934610/149‑1 (MQ=255)
                                                                                      tggttggttTTCCCCCTCCCCGGGTGGTAGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                               <  2:7610678/149‑1 (MQ=255)
                                                                                          tgtttttCGCCCACCCCGTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                                           <  1:504938/147‑1 (MQ=255)
                                                                                           ggttTTCGCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                                          <  1:5699035/150‑1 (MQ=255)
                                                                                           ggttTCCTCCCCCGCCGGGGGGGTGTGGGCGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                                          <  1:532533/150‑1 (MQ=255)
                                                                                            gttTGCTCCCTCGGCGGGGGGGAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAg                                                         <  1:8687361/150‑1 (MQ=255)
                                                                                              tttCTCCCGCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaa                                                       <  1:8870519/150‑1 (MQ=255)
                                                                                                   tccctcgccctGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                                   <  2:3328099/150‑1 (MQ=255)
                                                                                                     catctccggGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                                 <  1:3996844/148‑1 (MQ=255)
                                                                                                      ctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCTAACCCAAACCCTCACCCGCCAGGGGGAGGGGACCGAACGTGCCCAAACTTTGAGTGATATATTCTTATTTAGTGTTAAACAAGAAACGCAAATCggg                                                >  2:6042933/1‑150 (MQ=16)
                                                                                                             tgtgGGTGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                                         <  1:4051376/150‑1 (MQ=255)
                                                                                                             ggtgGGTGAGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                                         <  1:9742327/149‑1 (MQ=255)
                                                                                                             ggtgGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                                         <  2:7424935/149‑1 (MQ=255)
                                                                                                             ggtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                                         <  2:723858/149‑1 (MQ=255)
                                                                                                                gggcgcgggCCGGGGGGCGGGCCTCAGCGCGCCCGGTCACCCTCACCCGACCCCTCTCCCTCACGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATCTCTTCGTCTGTCTTGTAACAGAAGAACAGAACATCGGGCTTTTGCGcc                                      <  2:5676684/144‑1 (MQ=255)
                                                                                                                gggAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGCACCCCCACCCCAACCCACACCCTCAAGGGAGAGGGGCCCCAACCAGCCCAAAATTTGAGCGGTATATTCTTAGGTGTTGTGACCGAAGAACAGCAAAAAGGGCGGTTGGGcc                                      >  2:4110753/1‑150 (MQ=255)
                                                                                                                  gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTg                                    <  2:8309527/150‑1 (MQ=255)
                                                                                                                   agagGGCCGGGGGGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTGCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                                   <  1:9488717/150‑1 (MQ=255)
                                                                                                                    gtgggCGGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                                  <  1:3604377/148‑1 (MQ=255)
                                                                                                                      gggCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata                                <  2:6185699/150‑1 (MQ=255)
                                                                                                                       ggCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                               <  1:4400475/150‑1 (MQ=255)
                                                                                                                        gCCGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt                              <  2:3412935/150‑1 (MQ=255)
                                                                                                                          cGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                            <  2:5044886/150‑1 (MQ=255)
                                                                                                                           ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                           <  1:2452351/150‑1 (MQ=255)
                                                                                                                             ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga                         <  2:1184919/150‑1 (MQ=255)
                                                                                                                             ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga                         <  1:6302630/150‑1 (MQ=255)
                                                                                                                             ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga                         <  2:9155588/150‑1 (MQ=255)
                                                                                                                              gTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgag                        <  2:1844673/150‑1 (MQ=255)
                                                                                                                                   ggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                   <  1:3239988/150‑1 (MQ=255)
                                                                                                                                   ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                   <  2:4238710/150‑1 (MQ=255)
                                                                                                                                   gccATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCACGTGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                   <  1:2936518/148‑1 (MQ=255)
                                                                                                                                     cctCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGACTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc                 <  1:4047671/148‑1 (MQ=255)
                                                                                                                                     cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc                 <  2:4684237/150‑1 (MQ=255)
                                                                                                                                      ctCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa                <  2:804815/149‑1 (MQ=255)
                                                                                                                                       tCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCAt               <  1:6606622/150‑1 (MQ=255)
                                                                                                                                        cAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATg              <  2:7128917/150‑1 (MQ=255)
                                                                                                                                             cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg         <  1:4932462/150‑1 (MQ=255)
                                                                                                                                             cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg         <  1:4252161/150‑1 (MQ=255)
                                                                                                                                              gcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGCTCGAGCACAAATTTTGAGAGATATCTTCGTCTGTCTTGTAACCGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCggg        <  1:7681242/150‑1 (MQ=255)
                                                                                                                                              gcACGGGCACCCTCCCCCAAACCCAATCCCTCAAGGGAGGGGGGACCGCTGGAGCACAAAATTTGAGCGATATATTCATATGTCTTGTAACAGAAGAACAGAAAAACGGGCTTTTGAGCCAGAAGATTTCGGAGAAAAAGGAGGGCCggg        >  2:2351667/1‑150 (MQ=255)
                                                                                                                                                aCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGGGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGaa      >  2:4026196/1‑150 (MQ=255)
                                                                                                                                                aCGTACACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGaa      <  1:6939409/150‑1 (MQ=255)
                                                                                                                                                aCGATCACCCACACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCATTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGaa      <  2:190292/150‑1 (MQ=255)
                                                                                                                                                 cGTTCACCCTCACCCGAACCCTCTCCCTCACGGGAGAGGGGACCGCTCGAGCACCAATTTTGAGAGCTATCTTCTTCTGTCTCGTAACAGAAGTACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAt     <  2:5513558/150‑1 (MQ=255)
                                                                                                                                                    tcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGAc  <  1:726669/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGAC  >  NC_000913/164448‑164744

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: