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

TGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  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)
                                                             aaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCACTCCCCCAAGGGAGAGGGGACCGATCGAGca‑c                                                                                    >  1:1799595/1‑150 (MQ=255)
                                                                 cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCAAACCCACACCCCCAACGGGGGGGGGGCCGATAGAGCA‑Caaaa                                                                                >  1:3560967/1‑149 (MQ=25)
                                                                    cACCGACGCGACGTACTAAAATGTATTCTTACGGTTGTTTTCCCCCCTCCCCCTGTGGTTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAAtttt                                                                             <  1:2538491/150‑1 (MQ=255)
                                                                    cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCAACCACACCCAAACCCACACCCACAAGGGAGAGGGGACCCAACGAGCA‑CAAAattt                                                                             >  2:3878760/1‑150 (MQ=37)
                                                                     cccGACGCGACGGACTAAAATGTTTTCTTATGGTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTg                                                                            <  2:411475/149‑1 (MQ=38)
                                                                      cggacgggacgTACGAAAATGAATTCGTATGTTTGGTTTTCTCCATCGCCGTGTGGGTGTGGGCGGGGGTTTGGGAATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTga                                                                           <  2:3008704/148‑1 (MQ=255)
                                                                      ccgccgCGAGGGACGAAAAAGTTTTGGTAAGGTTGGTTTTCGCCCTCTCCGTGTGGGCGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTga                                                                           <  2:989074/150‑1 (MQ=255)
                                                                      ccGACGCGCCGTACCAAATATTATTCTTAGGGTTTGTTTTCACCCTGCCCGTGTGGTTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTga                                                                           <  1:3398326/150‑1 (MQ=255)
                                                                      ccGACGCGCCGGACTAAAATGTATTGTTTATGTTGTTTTTCCCCCTCTCCGTGGGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTga                                                                           <  2:3251727/150‑1 (MQ=255)
                                                                      ccGACGCGACGGACTAAAAATTATTCGTATGGTTGGTTTTCTCCCTCTCCGTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCCCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTga                                                                           <  1:3157682/150‑1 (MQ=38)
                                                                       cgacgcgacgGACTAAATTTTTTTCTTATGTTTTTTTTTCTCCCTCTCCCTGTGTTTGTGGGCGGGGGTGGGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTgag                                                                          <  1:656023/150‑1 (MQ=255)
                                                                       cgacgcgacgGACGAAATAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGGACCCGCCCCCAAACCCACAACCACAACCGCGAGGGGGACGCACGGGCA‑Caaaaatttag                                                                          >  1:604481/1‑143 (MQ=16)
                                                                       cgacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGTCCCACCCCCCAACCCACACCCCCAACGGCGTGGGGGCCGAGCGAGCC‑CAAAAATTgag                                                                          >  2:1151080/1‑150 (MQ=21)
                                                                       cgacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACACCCACAAGGGCGAGGGGGCCGCACGAGCACCAAAttttag                                                                           >  1:4133127/1‑148 (MQ=25)
                                                                          cgcgacgGACTAAAAAGTATTCTTATGGTTGGTTTTATCCCTCTCCCTTTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTgagaga                                                                       <  2:3073845/150‑1 (MQ=38)
                                                                          cgcgacgGACGAAAAAGTATTCTTAAGGTTGGTTTACTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTgagaga                                                                       <  2:262084/150‑1 (MQ=255)
                                                                           gcgacgGCGGAAAATTTTTTCGTTTGGTTGTTTTTCTCCTTCTCCCTGGGGGAGTGGGCGGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGat                                                                      <  2:4044595/150‑1 (MQ=255)
                                                                             gacggactAAAATGTATTGGTAAGGTTGGTTTTCTCCCTCTCCGTGTGGGATTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGatat                                                                    <  2:3421885/150‑1 (MQ=255)
                                                                              acggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGCCCCACCCCCCAACCCACACCCCCCACGGCGCCGGGCCCCCTCGCGCC‑CCAAATTTGAGAGAGAtc                                                                   >  1:2057285/1‑150 (MQ=12)
                                                                                 gacgAAAAATTATTCGTATGGTTGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATAtcttc                                                                <  1:2011807/150‑1 (MQ=255)
                                                                                   cgAAAACGTATACGTCAGGTTGCTTTACTCCCGCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATAtcttctt                                                              <  1:3438748/150‑1 (MQ=255)
                                                                                    gAAAATGTTTTCGTACGGTGGGTTTTCGCCCTCTCCGTGGGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATAtcttcttc                                                             <  1:2687260/150‑1 (MQ=255)
                                                                                        aaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTc                                                         <  2:2484298/150‑1 (MQ=255)
                                                                                              tCGTATGGTTGTTTTTCACCCTCTCCCTGTGGTTGTGGGCGGGGTTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                   <  1:4023742/150‑1 (MQ=255)
                                                                                                  acggttggttTTCGCCCTCTCCCTGTGGTAGTGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                               <  2:4088023/148‑1 (MQ=255)
                                                                                                   tgtttgtttTCCTCCCCCTCCGGGGGGTTGCGGGCGGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                              <  1:1285360/147‑1 (MQ=255)
                                                                                                   tggttgtttTTCTCCCACCCCCGGTGGTAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                              <  1:27338/149‑1 (MQ=255)
                                                                                                   cggttggttTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                              <  1:2033135/149‑1 (MQ=255)
                                                                                                     gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCGAACCCACTCCCGCAAGGGGGAGGGGCCCGCTCGAGCA‑CAAAGTTTGGGAGATATCTTCTTATGTCGTGGTACagacga                                            >  2:4207704/1‑150 (MQ=25)
                                                                                                     ggtggttTTCCCCCTCCCCCTGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                            <  1:1212575/148‑1 (MQ=255)
                                                                                                      tgggttTCCCCCCCCGCCGTGGGGTAGTGGGCCGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAAACCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                           <  1:3216906/148‑1 (MQ=255)
                                                                                                       tggttTTCTCCCTCTCGCTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                          <  1:2055821/150‑1 (MQ=255)
                                                                                                       tggttTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                          <  2:2065972/150‑1 (MQ=255)
                                                                                                       tggttTTCCCCCCCGCCCTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                          <  1:2372591/150‑1 (MQ=255)
                                                                                                          ttttCGCCCGCGCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                       <  2:3652795/150‑1 (MQ=255)
                                                                                                                 cctctcccaCAGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCTCCCCCACCCGAAACCACTCCCGCCATGGCGCGGGGCCCGTACGATCC‑CCAAGTTTGGGGGATTACGTCTTAGTGCTTGTATCCGCAGCACAGGAAATCgg                                >  2:1374367/1‑150 (MQ=255)
                                                                                                                 cctccccctGGGGGTGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                <  1:2929938/150‑1 (MQ=255)
                                                                                                                  ctctccctGTGGGAGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               <  2:604481/150‑1 (MQ=255)
                                                                                                                  ctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               <  2:3507321/150‑1 (MQ=255)
                                                                                                                  ctcgccctGTGGGAGCGGGCCGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCAACCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               <  2:406781/150‑1 (MQ=255)
                                                                                                                           tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                      <  1:3206629/150‑1 (MQ=255)
                                                                                                                           ggggAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                      <  2:258570/149‑1 (MQ=255)
                                                                                                                             ggAGAGGGCCGGGTTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                    <  1:3878259/150‑1 (MQ=255)
                                                                                                                                   ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCACACCCCCAAGGGAGAGGGGGCCCAACCACCA‑CAAAAATTGAGAGATATCATATTAAGTCTTGTAAAACAAGAACAAAAAAAAGGGCTTTTTCGCCGGAatat              >  2:978566/1‑150 (MQ=255)
                                                                                                                                    gCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATTTCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt             <  1:4133466/150‑1 (MQ=255)
                                                                                                                                     ggggggTTCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg            <  2:3492/148‑1 (MQ=255)
                                                                                                                                     cgggggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg            <  1:161722/148‑1 (MQ=255)
                                                                                                                                       ggggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:10832/150‑1 (MQ=255)
                                                                                                                                       ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:2129769/150‑1 (MQ=255)
                                                                                                                                       ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg          <  1:1241899/150‑1 (MQ=255)
                                                                                                                                         ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga        <  2:1880648/150‑1 (MQ=255)
                                                                                                                                              gggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa   <  1:3757864/150‑1 (MQ=255)
                                                                                                                                              gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa   <  1:2225198/150‑1 (MQ=255)
                                                                                                                                              gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa   <  1:2225224/150‑1 (MQ=255)
                                                                                                                                               ggCCTCAGCGCGCAAGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATCTCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa  <  1:1666054/150‑1 (MQ=255)
                                                                                                                                               ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑AAAATTTTGAGAGATATCTTCTTCTGTCTTGAAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa  <  2:1755540/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  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: ‑