Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 26.0% 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,5810AC26.0% 100.5 / 32.1 89intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (26/39);  new base C (2/22);  total (28/61)
Fisher's exact test for biased strand distribution p-value = 4.29e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164433‑164729
                                                                                                                                                    |                                                                                                                                                    
cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                     >  1:726916/1‑150 (MQ=255)
        gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcgc                                                                                                                                             >  2:2806940/1‑150 (MQ=255)
        gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcgc                                                                                                                                             >  2:2703663/1‑150 (MQ=255)
        gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcgc                                                                                                                                             >  2:2543747/1‑150 (MQ=255)
        gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAgcgcgc                                                                                                                                             >  2:903570/1‑150 (MQ=255)
             aatatttGAAAGGGCGTTATCCCAAACTTGTCTGGCCGGACGACCCGGCAATTAGTGAACCGACGCGACGGACGAAAATTTATTCTTATGTTTGTTTTTCTCCCTATCCGTGTGGGAGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGtt                                                                                                                                        <  2:2161623/144‑1 (MQ=255)
             aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtg                                                                                                                                        >  1:2739404/1‑149 (MQ=255)
                agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcac                                                                                                                                     <  2:3105792/150‑1 (MQ=255)
                           cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACACCCCaa                                                                                                                          >  2:2370942/1‑150 (MQ=38)
                            gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGACCCCCCCCCCAAc                                                                                                                         >  2:1715536/1‑150 (MQ=38)
                            gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCaaaa                                                                                                                         >  2:1171247/1‑149 (MQ=255)
                            gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACCCCCCAAc                                                                                                                         >  2:1171102/1‑150 (MQ=38)
                            gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCACACCCCAAc                                                                                                                         >  2:1715533/1‑150 (MQ=255)
                                 cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCCCCCCAAcccgca                                                                                                                    >  1:1435639/1‑147 (MQ=37)
                                 cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCCCCCCAAcccaca                                                                                                                    >  1:2976447/1‑147 (MQ=37)
                                 cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCCAAaccaca                                                                                                                    >  2:2892620/1‑147 (MQ=38)
                                   cAAACATGTGTGGCGGGACCCCCCGGAAAATACTGCCCCGACGCGAGGTACTAAAATTTATTGTTATGGTTGGTTTTCTCCCTCTCCCTGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                  <  2:1433083/150‑1 (MQ=25)
                                    aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGACCCGCCCCCACACCCACTccc                                                                                                                 >  1:102937/1‑150 (MQ=37)
                                        aggccgggCCGCACGCCCCGCAAAATACTCCACCCACGCCACGGACTAAAAATTATTCGTAAGTTTGTTTTACCCCCTCACCCTGTGGGAGCGGGCCGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:842551/144‑1 (MQ=255)
                                        aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                             <  2:777758/150‑1 (MQ=255)
                                         tttgtGGCCGGACGACCCGCAAAATACGCCCCCGACGCTACGCACTAAAATTTTTTCTTATGGTTGTTTTTCTCCCTCTCCGGGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                            <  1:1202608/146‑1 (MQ=255)
                                         tGTCTGGCCGGCCGCCCCGGCAAATCCTGCCCCGACGCGCCGGACGAAAAAGTATTCGTCCGGTGGGTTTTCTCCCTCTGCCTGGGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                            <  1:1383597/150‑1 (MQ=37)
                                         tGTCTGGCCGGACGCCCCGGAAAAAACTGCCCCGCCGCGACGGACTAAAAATTTTTCGTTGGGTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                            <  1:1202612/150‑1 (MQ=25)
                                          gTCGGGCCGGACGACCCGGCAAATCCTGCACCGCCGGGACGGACGAAAATGTTTTCGTATGGTTGGTTTTCCCCCTCTCCCTGTGGTAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                           <  2:524068/150‑1 (MQ=255)
                                           tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACACCCACAAggg                                                                                                          >  2:748453/1‑150 (MQ=38)
                                              ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                       <  2:488247/150‑1 (MQ=255)
                                               gCCGGACTACCCGGAAAATCCTCCCCCGACGCCCCGGACTAAAAATTATCCGTATGGTGGTTTTTCCCCATCTCCCTGTGGGAGTGGGCCGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                      <  1:2777577/150‑1 (MQ=255)
                                               gCCGGACGACCCGCAAAATACTGCACCGACGGGCCGGCCTAAAATTTTTTGTTATGTTTGGTTTTCCCCCTCGCCGGGGGGGTGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                      <  1:1018106/150‑1 (MQ=255)
                                                ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                     <  2:473237/150‑1 (MQ=255)
                                                ccGGACGACCCGGAAAATACTGCACCGCCGCGCGGGACGAAAAATTTTCCGTAAGGTTGTTTTGCCCCCCCTCCCTGTGGGCGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                     <  2:2110434/150‑1 (MQ=255)
                                                ccGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAATTATTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                     <  2:326941/150‑1 (MQ=38)
                                                ccGGACCACCCGGAAAATACTGCCCCGACGCGCCGGACTAAAAAGTATCCGAAGGGTTGGTTTCCCCCCTCCCCCTGGGGGTGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                     <  2:2382440/150‑1 (MQ=255)
                                                 cGGACGACCCGGCACATACGGCCCCGACGCGACGGCCGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCACTCCCTGTGGGAGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTATCCCTCTCCCTCAAGGGAGAggg                                                                                                    <  2:3076425/150‑1 (MQ=255)
                                                      gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCTCTCCCTCAAGGGCGTGGGGACCg                                                                                               >  2:2042547/1‑150 (MQ=255)
                                                      gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCTCACCCGAACCCACTCCCTCACGGGAGAGGGGACCg                                                                                               >  2:2042457/1‑150 (MQ=255)
                                                         ccgcaaaaGACTCCACCGCCGCGACGGACTAAAATTTTTTCGTATGTTTGTTTTACTCCATCACCCTGTGGGAGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                            <  2:906870/145‑1 (MQ=255)
                                                          cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCACCCAAAACCACAACCACAAGGGAGAGGGGGcccaaag                                                                                           >  2:2248587/1‑145 (MQ=25)
                                                           ggaaaaTCCTGCACCGACGCGACGGACGAAAAATTATGCTTAAGTTTGGTTTTCTCCCTCTCCGTTTGGTTTTTGGCCGTGGTTTGGGAATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                          <  2:966958/147‑1 (MQ=255)
                                                           ggCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                          <  2:1939994/150‑1 (MQ=255)
                                                               aaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACCCCCGATCCCGCACCCGCAAGGGGGAGGGGGCCGGTAGAGcac                                                                                      >  1:1029459/1‑150 (MQ=25)
                                                                     gCACCGCGGCGCCGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCCCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAttt                                                                                <  2:2172718/150‑1 (MQ=38)
                                                                      ccccgccgCGCCGGACTAAAATGTATTCGTATGGTTGGTTTACCCCCGCTCCGGGTGGGAGCGGGCGGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                               <  1:489050/148‑1 (MQ=255)
                                                                      cACCGCCGCGCCGGCCAAAAAAGTATTCGAAAGGTTGGTTTTCCCCCCCTCCCGGGGGTAGCGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGGGGGGACCGATCGAGCACAAAtttt                                                                               <  1:1098747/150‑1 (MQ=255)
                                                                        ccgccgCGACGGACTAAAATGTATTCGTATGGTTGTTTTTCCCCCTCTCCCGGGGGTTGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  1:1731093/150‑1 (MQ=255)
                                                                        ccGACGCGACGGACTAAAATTTATTCGTATGTTTGTTTTTCTCCCTCCCCCGGGGGTTTTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  1:2914994/150‑1 (MQ=255)
                                                                         cgccgcgacgGACGAAAAAGTCTGCGTCAGGTTGGTGTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                            <  2:2020666/147‑1 (MQ=38)
                                                                         cgacgcgacgGACGAAAAATTATTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                            <  2:2020773/150‑1 (MQ=255)
                                                                          gaggcgccgGACGAACATGTTTTGGTATGGTTGGTTTTCGCCCTCGCCGTTGGGGTGCGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                           <  2:1981507/147‑1 (MQ=255)
                                                                          gacgcgacgGACGAAAATTTTTGCGTACGTTTGGTTTTCTCCCTCCCCCTGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                           <  1:2002200/150‑1 (MQ=255)
                                                                          gacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                           <  2:3092554/150‑1 (MQ=255)
                                                                            cgcgccgGCCGCAAACGTAGTGGTACGGTTGGTTTGCGCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga                                                                         <  2:2466561/150‑1 (MQ=37)
                                                                               gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGAAACCAAACCCACTCCCCCAAGGGAGAGGGGACCCAACGCGCACAAAATTTTAGAGataa                                                                      >  2:2289299/1‑149 (MQ=25)
                                                                                   gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCACCGCAACCAAAAACGCTCCCACACCGGAGATGGGCCCCAACGCCCACAAAATTTGAGAAATAactag                                                                  >  2:1936951/1‑148 (MQ=16)
                                                                                    tcgAAAATGTATTCGTTCGGTGGGTTTTCTCCCTCTCCGTGGGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                                                                 <  2:1829483/149‑1 (MQ=255)
                                                                                    acgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCACCCCAACCCACACCCCCCCAGGGGGTGGGACCGCACCGGCAAAAATTTTTAGAGataaattat                                                                 >  2:930466/1‑144 (MQ=12)
                                                                                          atttATTCGTCCGGTGGGTTTCCCCCCTCTCCCGGTGGGTGCGGGCGGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                           <  1:456982/147‑1 (MQ=255)
                                                                                                 cGTATGGTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                    <  1:3020901/150‑1 (MQ=255)
                                                                                                 cGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                    <  2:2565012/150‑1 (MQ=255)
                                                                                                 cGTACGGTTGGTTTTCCCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                    <  1:3020888/150‑1 (MQ=255)
                                                                                                    aaGGTGGGTTTCCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                 <  1:1314987/150‑1 (MQ=255)
                                                                                                     cgggtggttTCCCCCCACTCCCTGGGGGCGTGGGCCGGGTGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCACGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTGCTTCGGTCTTGTAACagaa                                                <  2:2107629/146‑1 (MQ=255)
                                                                                                     aGGTTGGTTTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                <  2:772287/150‑1 (MQ=255)
                                                                                                       gttggttTTCTCCCTCTCCGTGTGGTAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                              <  1:458252/150‑1 (MQ=255)
                                                                                                       gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGGGCACCCGCACCCGACCCCGCTCCCGCCAGGGCCGGGGGGCCGAGCGAGCCCCAAGTTTGTGAGCTATCTTATTCAGTCGTGGTACagacga                                              >  2:1311431/1‑150 (MQ=12)
                                                                                                       gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCTCACAATACCACTCCCGCAAGGGCGTGGGGGACGATCGAGCACAAATGTTGGGTGATTTCTTCTTCGGTCGTGGTACagacgg                                              >  2:2854924/1‑147 (MQ=12)
                                                                                                       gttggtgTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                              <  1:347243/150‑1 (MQ=255)
                                                                                                        ttggttTTCTCCCTCGCCCTGTGGGATAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                             <  1:2247169/150‑1 (MQ=255)
                                                                                                            ttttCTCCCTCGCCGTGGGGGTTTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                         <  1:1220756/150‑1 (MQ=255)
                                                                                                                   cctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCTCCCCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGCGCGAGCCCAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGGAAAGCgg                                  >  2:2318858/1‑150 (MQ=38)
                                                                                                                   cctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCGCACCCAAAACCTCTCCCGCAAGGGGGAGGGGACCGATCGAGCACAAAATTTGTGAGATATATTCTTCGGTCTTGTAAAAGAAAAACAGaaaaaagg                                  >  2:2318863/1‑146 (MQ=25)
                                                                                                                         ccggtgGTAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                            <  1:1754687/147‑1 (MQ=255)
                                                                                                                           tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACAACCAAAACCACAACCCCAAGGGAGGGGGGGCCCCTCCAGCCCAAAATTTTTGAGAAAAATTATTAAGTGTTTTTACAGAAAAAAACAAAAAAGGGCgttttg                          >  2:2541207/1‑148 (MQ=255)
                                                                                                                              gggAGCGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                       <  1:2661276/150‑1 (MQ=255)
                                                                                                                              gggAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                       <  1:2661307/150‑1 (MQ=255)
                                                                                                                                      gCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt               <  2:703508/150‑1 (MQ=255)
                                                                                                                                      gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTGGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt               <  2:1967481/150‑1 (MQ=255)
                                                                                                                                       cgggggggCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg              <  1:2995738/148‑1 (MQ=255)
                                                                                                                                       ccGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg              <  1:591509/150‑1 (MQ=255)
                                                                                                                                        gggggTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc             <  2:2976393/149‑1 (MQ=255)
                                                                                                                                        gggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc             <  2:1435639/149‑1 (MQ=255)
                                                                                                                                        gggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc             <  2:2976447/149‑1 (MQ=255)
                                                                                                                                        cGGGGTGCGCCCCTCAGCGCTCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc             <  1:2892620/150‑1 (MQ=255)
                                                                                                                                        cGGGGGGCGGGCCCCCGCGCGCCCGTGCACCCTCACCCGACCCCACTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGCCGGGTAACAGCAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc             <  2:388949/150‑1 (MQ=255)
                                                                                                                                         ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGATAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  1:2032075/150‑1 (MQ=255)
                                                                                                                                          gggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg           <  1:436427/150‑1 (MQ=255)
                                                                                                                                                gggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa     <  1:635072/150‑1 (MQ=255)
                                                                                                                                                 ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa    <  2:1918131/150‑1 (MQ=255)
                                                                                                                                                  gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg   <  1:2370942/150‑1 (MQ=255)
                                                                                                                                                   cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc  <  1:2731443/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164433‑164729

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: