Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 22.6% intergenic (+47/‑149) 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,5810AC22.6% 132.8 / 15.3 99intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (24/51);  new base C (4/20);  total (28/71)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  NC_000913/164435‑164727
                                                                                                                                                  |                                                                                                                                                   
tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCTAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                  <  2:1436991/150‑1 (MQ=255)
 gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAg                                                                                                                                                 >  1:809722/1‑150 (MQ=255)
             tagaTTAATGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAAAATTATTCGTATGTTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTtca                                                                                                                                     <  2:2044239/149‑1 (MQ=37)
             gagaTGAACGGGCTTTAGCCCACAAATGTCGGGCGGGACGACCCGGAAAATACTGCACCGCCGCGACGGACGAAAAATTTTTCGTACGGTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGCCATCAGCGCGCACGTtca                                                                                                                                     <  2:233831/150‑1 (MQ=25)
             gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtca                                                                                                                                     <  2:2776201/150‑1 (MQ=255)
             gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtca                                                                                                                                     <  2:2776278/150‑1 (MQ=255)
              agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggc                                                                                                                                    >  1:411475/1‑146 (MQ=255)
                aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAATAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAACGCGCAGGtggaggc                                                                                                                                  >  2:3280071/1‑144 (MQ=38)
                aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccc                                                                                                                                  >  1:1413019/1‑150 (MQ=255)
                  gAATGGGGGTTTTCCAAAATTTGTCTGGCCGGAGGCCCGGGAAACTACTGCACCGAGGCGCGGGACTAAAATTTTTTGTTATTGTTGTTTTTCTCCCTCTCCGTGTGGGATTGGGCGGGGGTTTGGGAATCAGCGCGCACGTtcaccctc                                                                                                                                <  2:1127772/150‑1 (MQ=255)
                   aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacccaca                                                                                                                               >  1:3876526/1‑150 (MQ=255)
                      gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggacccgcaccc                                                                                                                            >  2:89373/1‑150 (MQ=255)
                          gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCCaaa                                                                                                                        >  2:1474040/1‑149 (MQ=38)
                           ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGACCCACACCCCAAcc                                                                                                                       >  1:4211319/1‑150 (MQ=38)
                           ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcgcccgcacccccaac                                                                                                                       >  1:1624298/1‑145 (MQ=255)
                           ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACACCCaaaac                                                                                                                       >  1:1624297/1‑148 (MQ=38)
                           ttATCCCAAACAGGGCTGGCCGGACGACCCGGCAAAAACTGCCCCGACGCGCCGTACTAAAAAGTATTCTTACGGTGGTTTTTCTCCCTCCCCCTGTGGGCGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                       <  2:291053/150‑1 (MQ=255)
                                ccAAACTTGTGGGGCCGCACCCCCCGCAAAATACGCCCCCGCCGCGCCGGACCAAATAGTTTTCGTATGGTTGGTTTTCCCCCCCCCCCTGGGGGTGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                  <  1:3066291/150‑1 (MQ=255)
                                  aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGCGCCCGCCCCCGCACCCCCGccc                                                                                                                >  2:2982695/1‑150 (MQ=37)
                                      aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCACCCGCCCCCAAACCCGCAcccaccc                                                                                                            >  2:2510068/1‑146 (MQ=37)
                                       tgtggggCCGGACCCCCCGCCAAATACTGCACCGACGCCCCGCCCTAAACAGTTTTGGTATGGTGGGTTTTCGCCCTCCCCGGGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                           <  2:3943832/145‑1 (MQ=255)
                                        gtgtGGCCGGACGACCCGGAAAATACTGCACCTCCGCGACGGACTAAAAAGTTTTCTTACGGTTGGTTTTCTCCCACTCCCTGTGGTAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                          <  2:4068293/147‑1 (MQ=37)
                                         tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCAAACCCAAAACCACACCCAAAAggg                                                                                                         >  2:244875/1‑150 (MQ=37)
                                          cgggCCGGACGCCCCGGCAAATACGGCACCGACGCGCCGGACGAAAAAGTTTTCGTAAGTTTGGTTTTCGCCCTCTCCGTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                        <  2:547689/148‑1 (MQ=255)
                                          cgggCCGGACGACCCGGCAAATCCTGCACCGACGCGAGGGACGAAAAAGTTGTCGTATGGTTGGTTTTCTCCCTCTCCCGGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTTACCCTCTCCCTCAATGga                                                                                                        <  2:3482144/148‑1 (MQ=37)
                                             gcgggACGACCGGGCAAATACTGCACGGCCGCGCGGTACGAAAAAGTTTTGTTATGGGTGTTTTTCTCCCTCTCCCTGTGGGCGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                     <  2:1307330/147‑1 (MQ=25)
                                             gCCGGACGACCCGGCAAATACTGCCCCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCGCCTTCGCCCTGTGGGTGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                     <  2:2574445/150‑1 (MQ=255)
                                             cccGGCCGACCCGTCAAAAACTGCACCGCCGCGACGGACCAAAATGTATTCGTACGGTTGTTTTTCACCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                     <  1:3767623/149‑1 (MQ=255)
                                              ccGGACGACCCGTAAAATCCTGCACCGAGGCGACGGACGAAAAATTTTACTTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                    <  1:3545120/150‑1 (MQ=37)
                                              ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAAGTTTCCGTACGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                    <  2:3696314/150‑1 (MQ=38)
                                              ccGGACCACCGGCCAACTCCTGCACGGAGGCGAGGTACGAAATAGTTCTCTATTGTTTGGTTTTCGCCCTCCCCCTGGGGGTGTGGGCGGGGTTTTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                    <  2:1627914/150‑1 (MQ=255)
                                                 gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAAACCAAAACCACAAGGGAggggggg                                                                                                 >  2:3426600/1‑149 (MQ=37)
                                                 gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCGCACCCAAACCCACTCCCGCACGGGAGAggggc                                                                                                 >  2:1029593/1‑149 (MQ=37)
                                                 cacgacCCGCAAAATCCTCCACCTACGCCCCGGACTAAAAATTTTTCTTATGTTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGGCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                 <  1:2708890/149‑1 (MQ=21)
                                                  ccgacCCGGCAAATTCTGCCCCGCCGCGACGCACGAAAAAGTATTCTTACGGTTGGTTTTCTCCCACTCCGTTTGTGAGCGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                                <  2:2009846/149‑1 (MQ=255)
                                                      cccGGAAATTAGTGCACCGCCGCGACGGACGAAAAAGTTTTCGTATGGTTGTTTTTCCCCCCCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                                                            <  2:441466/150‑1 (MQ=38)
                                                       ccGGCAAAGACTGCCCCGCCGCGACGGACTAACAAGTATTCGAACGGGGGGTTTTCTCCCGCACCCTGGGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                           <  2:645723/150‑1 (MQ=255)
                                                       ccGGAAAACACGGCACCGCCGCGCCGGACGAAAAAGTTTGCGTACGTTTGGTTTACGCCCTCTCCCTGGGGGAGCGGGCCGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                           <  2:1634120/150‑1 (MQ=255)
                                                         ggCAAATACTGCCCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                         <  2:2411474/150‑1 (MQ=255)
                                                          gCAAATACTGCACCGACGCGCCGGACGAACAAGTATTGGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                        <  2:3876526/150‑1 (MQ=255)
                                                            aaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCAACCGAACCCACTCCCACCCGGGCGAGGGGCCCGATCgcgcc                                                                                      >  1:3125608/1‑149 (MQ=25)
                                                            aaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCCAAACCACAACCACAACGGGGTGGGGGCCGATAGAGcc                                                                                      >  2:291223/1‑149 (MQ=25)
                                                            aaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCCCCCCAACCCACACCCCCAACGGGGGGGGGGCCCAACGaccc                                                                                      >  1:3697680/1‑147 (MQ=21)
                                                             aaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCACTCCCCCAAGGGAGAGGGGACCGATCGAGcac                                                                                    >  1:1799595/1‑150 (MQ=255)
                                                                 cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCAAACCCACACCCCCAACGGGGGGGGGGCCGATAGAGCACaaaa                                                                                >  1:3560967/1‑149 (MQ=25)
                                                                    cACCGACGCGACGTACTAAAATGTATTCTTACGGTTGTTTTCCCCCCTCCCCCTGTGGTTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                             <  1:2538491/150‑1 (MQ=255)
                                                                    cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCAACCACACCCAAACCCACACCCACAAGGGAGAGGGGACCCAACGAGCACAAAattt                                                                             >  2:3878760/1‑150 (MQ=37)
                                                                     cccGACGCGACGGACTAAAATGTTTTCTTATGGTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                            <  2:411475/149‑1 (MQ=38)
                                                                      cggacgggacgTACGAAAATGAATTCGTATGTTTGGTTTTCTCCATCGCCGTGTGGGTGTGGGCGGGGGTTTGGGAATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:3008704/148‑1 (MQ=255)
                                                                      ccgccgCGAGGGACGAAAAAGTTTTGGTAAGGTTGGTTTTCGCCCTCTCCGTGTGGGCGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:989074/150‑1 (MQ=255)
                                                                      ccGACGCGCCGTACCAAATATTATTCTTAGGGTTTGTTTTCACCCTGCCCGTGTGGTTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:3398326/150‑1 (MQ=255)
                                                                      ccGACGCGCCGGACTAAAATGTATTGTTTATGTTGTTTTTCCCCCTCTCCGTGGGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  2:3251727/150‑1 (MQ=255)
                                                                      ccGACGCGACGGACTAAAAATTATTCGTATGGTTGGTTTTCTCCCTCTCCGTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCCCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                           <  1:3157682/150‑1 (MQ=38)
                                                                       cgacgcgacgGACTAAATTTTTTTCTTATGTTTTTTTTTCTCCCTCTCCCTGTGTTTGTGGGCGGGGGTGGGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                          <  1:656023/150‑1 (MQ=255)
                                                                       cgacgcgacgGACGAAATAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGGACCCGCCCCCAAACCCACAACCACAACCGCGAGGGGGACGCACGGGCACaaaaatttag                                                                          >  1:604481/1‑143 (MQ=16)
                                                                       cgacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGTCCCACCCCCCAACCCACACCCCCAACGGCGTGGGGGCCGAGCGAGCCCAAAAATTgag                                                                          >  2:1151080/1‑150 (MQ=21)
                                                                       cgacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACACCCACAAGGGCGAGGGGGCCGCACGAGCACCAAAttttag                                                                           >  1:4133127/1‑148 (MQ=25)
                                                                          cgcgacgGACTAAAAAGTATTCTTATGGTTGGTTTTATCCCTCTCCCTTTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga                                                                       <  2:3073845/150‑1 (MQ=38)
                                                                          cgcgacgGACGAAAAAGTATTCTTAAGGTTGGTTTACTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga                                                                       <  2:262084/150‑1 (MQ=255)
                                                                           gcgacgGCGGAAAATTTTTTCGTTTGGTTGTTTTTCTCCTTCTCCCTGGGGGAGTGGGCGGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                      <  2:4044595/150‑1 (MQ=255)
                                                                             gacggactAAAATGTATTGGTAAGGTTGGTTTTCTCCCTCTCCGTGTGGGATTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                                    <  2:3421885/150‑1 (MQ=255)
                                                                              acggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGCCCCACCCCCCAACCCACACCCCCCACGGCGCCGGGCCCCCTCGCGCCCCAAATTTGAGAGAGAtc                                                                   >  1:2057285/1‑150 (MQ=12)
                                                                                 gacgAAAAATTATTCGTATGGTTGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                                <  1:2011807/150‑1 (MQ=255)
                                                                                   cgAAAACGTATACGTCAGGTTGCTTTACTCCCGCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                              <  1:3438748/150‑1 (MQ=255)
                                                                                    gAAAATGTTTTCGTACGGTGGGTTTTCGCCCTCTCCGTGGGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                             <  1:2687260/150‑1 (MQ=255)
                                                                                        aaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                         <  2:2484298/150‑1 (MQ=255)
                                                                                              tCGTATGGTTGTTTTTCACCCTCTCCCTGTGGTTGTGGGCGGGGTTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                   <  1:4023742/150‑1 (MQ=255)
                                                                                                  acggttggttTTCGCCCTCTCCCTGTGGTAGTGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                               <  2:4088023/148‑1 (MQ=255)
                                                                                                   tgtttgtttTCCTCCCCCTCCGGGGGGTTGCGGGCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                              <  1:1285360/147‑1 (MQ=255)
                                                                                                   tggttgtttTTCTCCCACCCCCGGTGGTAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                              <  1:27338/149‑1 (MQ=255)
                                                                                                   cggttggttTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                              <  1:2033135/149‑1 (MQ=255)
                                                                                                     gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCGAACCCACTCCCGCAAGGGGGAGGGGCCCGCTCGAGCACAAAGTTTGGGAGATATCTTCTTATGTCGTGGTACagacga                                            >  2:4207704/1‑150 (MQ=25)
                                                                                                     ggtggttTTCCCCCTCCCCCTGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                            <  1:1212575/148‑1 (MQ=255)
                                                                                                      tgggttTCCCCCCCCGCCGTGGGGTAGTGGGCCGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAAACCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                           <  1:3216906/148‑1 (MQ=255)
                                                                                                       tggttTTCTCCCTCTCGCTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                          <  1:2055821/150‑1 (MQ=255)
                                                                                                       tggttTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                          <  2:2065972/150‑1 (MQ=255)
                                                                                                       tggttTTCCCCCCCGCCCTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                          <  1:2372591/150‑1 (MQ=255)
                                                                                                          ttttCGCCCGCGCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                       <  2:3652795/150‑1 (MQ=255)
                                                                                                                 cctctcccaCAGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCTCCCCCACCCGAAACCACTCCCGCCATGGCGCGGGGCCCGTACGATCCCCAAGTTTGGGGGATTACGTCTTAGTGCTTGTATCCGCAGCACAGGAAATCgg                                >  2:1374367/1‑150 (MQ=255)
                                                                                                                 cctccccctGGGGGTGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                <  1:2929938/150‑1 (MQ=255)
                                                                                                                  ctctccctGTGGGAGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               <  2:604481/150‑1 (MQ=255)
                                                                                                                  ctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               <  2:3507321/150‑1 (MQ=255)
                                                                                                                  ctcgccctGTGGGAGCGGGCCGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCAACCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               <  2:406781/150‑1 (MQ=255)
                                                                                                                           tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                      <  1:3206629/150‑1 (MQ=255)
                                                                                                                           ggggAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                      <  2:258570/149‑1 (MQ=255)
                                                                                                                             ggAGAGGGCCGGGTTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                    <  1:3878259/150‑1 (MQ=255)
                                                                                                                                   ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCACACCCCCAAGGGAGAGGGGGCCCAACCACCACAAAAATTGAGAGATATCATATTAAGTCTTGTAAAACAAGAACAAAAAAAAGGGCTTTTTCGCCGGAatat              >  2:978566/1‑150 (MQ=255)
                                                                                                                                    gCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATTTCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt             <  1:4133466/150‑1 (MQ=255)
                                                                                                                                     ggggggTTCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg            <  2:3492/148‑1 (MQ=255)
                                                                                                                                     cgggggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg            <  1:161722/148‑1 (MQ=255)
                                                                                                                                       ggggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:10832/150‑1 (MQ=255)
                                                                                                                                       ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:2129769/150‑1 (MQ=255)
                                                                                                                                       ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:1241899/150‑1 (MQ=255)
                                                                                                                                         ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga        <  2:1880648/150‑1 (MQ=255)
                                                                                                                                              gggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa   <  1:3757864/150‑1 (MQ=255)
                                                                                                                                              gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa   <  1:2225198/150‑1 (MQ=255)
                                                                                                                                              gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa   <  1:2225224/150‑1 (MQ=255)
                                                                                                                                               ggCCTCAGCGCGCAAGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATCTCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa  <  1:1666054/150‑1 (MQ=255)
                                                                                                                                               ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCAAAAATTTTGAGAGATATCTTCTTCTGTCTTGAAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa  <  2:1755540/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  NC_000913/164435‑164727

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: