Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 17.1% 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,5810AC17.1% 150.4 / 14.7 96intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (27/51);  new base C (4/14);  total (31/65)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164433‑164730
                                                                                                                                                    |                                                                                                                                                     
cGTGAGGTGCAAAAAGAGATGAACGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      <  2:912589/150‑1 (MQ=255)
cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      >  2:1283904/1‑150 (MQ=255)
cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      <  2:3723365/150‑1 (MQ=255)
cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      >  2:3412140/1‑150 (MQ=255)
      gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                                >  1:252046/1‑150 (MQ=255)
      gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                                >  1:2439683/1‑150 (MQ=255)
            aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGt                                                                                                                                          >  1:2071665/1‑150 (MQ=255)
              ctagaTGAAGGGGGGTTAGCCCAAACCTGTCTGGCCGGACGCCCCGGCAAATACTGCCCCGAGGGGAGGGAGTAAAATTTTTTGGTATGTTTGGTTTTCGCCCACGCGCTGTGGGTTTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtc                                                                                                                                        <  2:2097147/148‑1 (MQ=255)
              agatATGACAGGGCTTTATCCCAAACAGGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGATTGGGCCGGGGTGCGGGCATCAGCGCGCACGTtc                                                                                                                                        <  2:2250362/150‑1 (MQ=38)
                  aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTtcaccc                                                                                                                                    >  2:3097106/1‑150 (MQ=255)
                  aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGtcaccc                                                                                                                                    >  2:2997873/1‑150 (MQ=255)
                  aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTggaccc                                                                                                                                    >  1:3667663/1‑150 (MQ=255)
                    gatcgggCGTTATCCCAAACTTGTGTGGCGGGACGACCCGGCAAATACTGCACCTACGCGAGTGACTAAAATTTATTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTTAGGGCCGGGGTGTGGGCATCAGCGCGCACGTtcaccctc                                                                                                                                  <  2:915362/146‑1 (MQ=25)
                     aatgggCTTTATCCCAAACATGGCTGGCGGGACGACCCGGCAAATACTGCACCGACGCGAGGGACGAAAAATTTTTCGTATGGTTGTTTTTCTCCCTCCCCCTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcaccctca                                                                                                                                 <  2:461854/147‑1 (MQ=255)
                      aaGGGCGTTATCCCAAACTTGTCTGGCCGGACGACCCGGAAATTACTGCACCGACGCGACGGACGAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                                <  2:2085075/150‑1 (MQ=38)
                            gTTATCCCAAACATGTCTGGCCGGACGACCCGGCACATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGGCACCCCCCCCCCAAc                                                                                                                          >  2:2172603/1‑150 (MQ=38)
                             ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCACCCaaaac                                                                                                                         >  1:1628778/1‑148 (MQ=255)
                             ttACCCCAACCTGGTCGGGCGGGACGCCCCGGAAAATACTGCCCGGACGCGACGGACTAAAATGTTTTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                         <  2:2293466/150‑1 (MQ=25)
                                 cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGACCCACCCCCAAAcccaca                                                                                                                     >  1:1387512/1‑147 (MQ=38)
                                 cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCACCCACACCCAAAaccaca                                                                                                                     >  1:2362583/1‑147 (MQ=38)
                                  ccAAACTTGTGGGGCGGGCCGACCCGCAAAATACTGCACGGACGGTACGCACTAAAAATTATTCGTAGGGTTGTTTTTCTCCCTCCCCGGGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                    <  2:377983/150‑1 (MQ=255)
                                  ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGACCCCCACCCAAAaccccac                                                                                                                    >  1:1772929/1‑146 (MQ=38)
                                  ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTGCCCCCCCCCCCACCCCtctc                                                                                                                    >  2:3283247/1‑150 (MQ=37)
                                  ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGACCCCCACCCCAAcccacaa                                                                                                                    >  2:1470220/1‑146 (MQ=38)
                                  ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCTCCCGCACCCAAAaccacaa                                                                                                                    >  1:2566152/1‑146 (MQ=38)
                                  ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCCCCCACAcccccgc                                                                                                                    >  1:2127143/1‑146 (MQ=37)
                                   cAACCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCCCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                   <  2:979399/150‑1 (MQ=255)
                                   cAAACTGGTCTGGCCGGACGACCCGGAAAATACTGCCCCGACGCGACGGACTAAAATTTATTCGTACGGTTGTTTTTCTCCCTCCCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                   <  2:400521/150‑1 (MQ=37)
                                    caccATGCCGGGCCGGCCGCCCGGCCAACTACGGCACCTACGCGAGGGCCGAAAATTTATTCTTAGGGTGGGTTTTCGCCCTCTCCGTGGGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                                                                                  <  2:387806/147‑1 (MQ=255)
                                    aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                                                                                  <  2:3069829/150‑1 (MQ=255)
                                    aaaCAGGTGGGGCCGGCCGACGCGGAAATTCCGGCACCTAGGCTTGGCACTAAAATGTATTCGTTTGTTTGTTTTTCGCCCTCCCCTTGTGGGTGCGGGCGGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                                                                                  <  2:387813/150‑1 (MQ=255)
                                        aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCAACCTAAACCACACCCACaa                                                                                                              >  2:2401631/1‑150 (MQ=37)
                                          gtgtGGCGGGACGCCCGGGCAAATCCGGCCCCGACGGGAGGGACGAAAATGTTTTCTTAGGGTTGGTTTCCTCCCCCGCCCGGTGGGCGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                            <  2:2854694/147‑1 (MQ=255)
                                          gTCTGGCCGGACGACCCGGCAAATAGTGCACCGACGCGACGGACGAAAAAGTTTTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                            <  2:1916043/150‑1 (MQ=255)
                                          gTCGGGCCGGCCGACCCGGAAAATCCTGCCCCGACGCGCCGGACTAAAAATTATGCGTTTGGTTGGTTTTCTCCCTCGCCCTGGGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                            <  1:3420204/150‑1 (MQ=255)
                                             gggCCGGACGACCCGCAAAATCCTGCCCCGACGCGACGGACCAAAAATTATTCTTATGGTTGTTTTTCCCCCTCGCCGTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                                         <  1:1972141/149‑1 (MQ=255)
                                              ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAACCCACACCCACCAGGggga                                                                                                        >  1:3088393/1‑150 (MQ=37)
                                              ggCCGGACGACCCGGAAAATACGCCCCCGACGCGACGGACGAAAAAGTTTTGGTTTGGTTGGTTTTCACCCTCCCCCTGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                        <  2:3206587/150‑1 (MQ=255)
                                               gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGCGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGTCCCGCACCCGAAACCACTCCCGCCAggggggg                                                                                                       >  1:1913013/1‑146 (MQ=37)
                                               gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCTAACCCTCACCCACCAGGGagag                                                                                                       >  2:3081870/1‑150 (MQ=255)
                                               gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCCCCCGAAACCTCTCCCACACgggggcg                                                                                                       >  2:3081373/1‑146 (MQ=37)
                                                cggtacgccCCGCAAAATACTGCACCGACGCGAGGTACTAAAATTTTTTCGTATGGTTGGTTTTCTCCCTCTCCTTGTGGGTGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                      <  1:2755702/146‑1 (MQ=255)
                                                 cGGCCGCCCGGGCAACTCCGGCCCCGCCGCGCCGGACGACACCGTTGGCGTAGGGTGGGTTTTCGCCCGCCCCGGGGGGGGGCGGGCGGGGGGGCGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                     <  1:2957533/150‑1 (MQ=255)
                                                   gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGCACCCACACCCAAACCCTCTCCCGCCAggggggggggg                                                                                                   >  1:1606071/1‑149 (MQ=37)
                                                   gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACCCCCCAACCCACAACCACAAGGGAggggggg                                                                                                   >  1:2015819/1‑149 (MQ=25)
                                                     cgccccGGCATATACGGCACCGACGCGCCGGACGAAAACGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                 <  2:3746142/147‑1 (MQ=38)
                                                         ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCACCCCCCAACCCCCGGCCGAAAGGGCGCGGGGACCgcgc                                                                                             >  2:1587480/1‑147 (MQ=25)
                                                            gCAAATCCGGCCCCGCCGCGCCGGACGAAACGGTATGGGTCTGGTGGGTTTTCGCCCCCGCCGTGGGGGCGCGGGCGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                          <  2:772522/150‑1 (MQ=255)
                                                            gCAAATACTGCCCCGCCGCGCCGGACGAAAACGTTTTCGTCCGGTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                          <  2:3544166/150‑1 (MQ=37)
                                                            gCAAAGACGGCCCCGCCGCGACGGACGAAAAAGTATTCGTAGGGTTGGTTTCCCCCCTCGCCCTGTGGGAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                          <  2:2945254/150‑1 (MQ=255)
                                                              aaaTCCGGCACCTCCGCGACGGACTAAAATTTTTGCTTATGGTTGTTTTTCCCCCCCGCCGTGTGGGTGCGGGCGGGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                        <  1:576923/150‑1 (MQ=255)
                                                              aaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCACACCCCCAACCCCACCCCCCAGGGGGGGGGGGCCGCTCgcgcc                                                                                        >  1:2162665/1‑149 (MQ=21)
                                                                      ccccgccgCGCCGTACGAAAAAGTTTTCTTACGGTTGTTTTACCCCCTCTCCCGGGGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                <  2:1926968/148‑1 (MQ=255)
                                                                      cACCGACGCGCCGGACTAAAATTTATTGTTAAGGTTGGTTTTCTCCCTCTCCGTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                <  2:3570652/150‑1 (MQ=255)
                                                                      cACCGACGCGAGGTACTAAAATTTTGTCGTATGGTTGGTTTTCTCCCTCCCCGTGGGGGTTTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                <  1:2580415/150‑1 (MQ=255)
                                                                      cACCGACGCGACGTACTAAAAATTTTGCTTATGTTGGGTTTTCTCCCTCTCCCGGTGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                <  1:2466674/150‑1 (MQ=37)
                                                                       aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCGAACCCACTACCACCAGGGAGAGGGGCCCCAGCGCGCACaaaaattt                                                                               >  1:926401/1‑145 (MQ=25)
                                                                        ccGACGCGCCGGACTAAAAAGTATTCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  2:252046/150‑1 (MQ=38)
                                                                        ccGACGCGCCGGACGAAAAATTATTCTTAAGTTTTGTTTTCTCCCTCTCCCTTTGGGCTTGGGCCGGTGTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACACTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  2:182005/150‑1 (MQ=255)
                                                                        ccGACGCGACGTCCGAAAAATTATTCTTACGGTTGTTTTTCTCCCTCTCCGTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  2:3858420/150‑1 (MQ=38)
                                                                        ccGACGCGACGGACTAAAAATTATTCGTATGGTTGGTTTTCTCCCCCTCCGTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  1:2966773/150‑1 (MQ=38)
                                                                        ccGACGCGACGGACGAAAAAGTTTTCGTATGGTTGGTTTTCGCCCTCGCCGTGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  2:2439683/150‑1 (MQ=255)
                                                                        ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                              <  2:2071665/150‑1 (MQ=255)
                                                                           acgcgacgGACGAAAATGTATTCGTACGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                           <  2:3726567/150‑1 (MQ=255)
                                                                               gacggactAAAATTTTTTCGTTTGGTTGGTTTTCCCCCTCCCCCTGTGGTAGAGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                                       <  1:267055/150‑1 (MQ=255)
                                                                                tcggacgAAATAGTTTTCGTTTGGTGGGTTTTCTCCCTCTCGGGGGGGGTGTGGGGGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                      <  1:2622196/149‑1 (MQ=255)
                                                                                     cgAAAATGTATTCGTAAGGTTGTTTTACTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                                 <  1:2342949/150‑1 (MQ=255)
                                                                                      gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCACACCCAAAACAAAAACCACAAGGGAGATGGGGCCCCTCGAGGCCAAAAATTGAGAGatataataaaa                                                                >  1:1723538/1‑143 (MQ=16)
                                                                                           aGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACACAAACCCACACCCGAAAGGGAGGGGGGACCCATAGAGGACAACTTTTGTGAGATATCTTATTATGTCt                                                           >  2:2655425/1‑150 (MQ=25)
                                                                                                tCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                      <  1:3412140/150‑1 (MQ=255)
                                                                                                  gtcgggttgggtTTCGCCCTCTCCGTGGGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                    <  1:3593298/146‑1 (MQ=255)
                                                                                                     aGGTTGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                 <  1:2757479/150‑1 (MQ=255)
                                                                                                       tttgtttTTCTCCCTCCCCCTGTGGGTGCGGGCCGGGTTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGTGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                               <  1:2112707/149‑1 (MQ=255)
                                                                                                         tggttTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                             <  1:1703734/150‑1 (MQ=255)
                                                                                                         tggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACAACCAAAACCACGACCACCAGGGGGAGGGGACCGCTCGAGCCCCCAGTGTGAGGGATATCTTAGTGTGTGTTGTGACAGAAGTAc                                             >  2:1275698/1‑150 (MQ=12)
                                                                                                          ggttTTCTCCCTCTCCGGGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                            <  2:3667663/150‑1 (MQ=255)
                                                                                                          ggttTTCCCCCACTCCCTGTGGGCGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTGTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                            <  1:4018247/150‑1 (MQ=255)
                                                                                                            tggtagcccgctcgctGTGTGTTAGGGCGGGTGTGAGGCCATCAGCGCGATAGTTCACCCTCACCCTAACCCTCTCCCTCATGGGAGAGGGGAACGATAGAGCACAAATTTTTTGAGATATCTTCTTCTGTCTTGTAACATAATAACAGa                                          <  1:2176256/144‑1 (MQ=255)
                                                                                                                tccgatctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAt                                      <  1:70675/145‑1 (MQ=255)
                                                                                                                  ccccccccctGTGGGTGTGGGCCGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                    <  1:488405/144‑1 (MQ=255)
                                                                                                                   cctctccctGGGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                   <  1:2150835/150‑1 (MQ=255)
                                                                                                                    ctctccctGTGGGTGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                  <  1:2997873/150‑1 (MQ=255)
                                                                                                                            gggggACAGGCCCGGTTTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcg                          <  2:494591/148‑1 (MQ=255)
                                                                                                                             tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                         <  2:4115496/150‑1 (MQ=255)
                                                                                                                               ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  2:2362583/150‑1 (MQ=255)
                                                                                                                                gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACTCCCACAAGGGCGAGGGGCCCGATCGAGCCCAAAATTTGAGAGATAAAATATTATGTCTTGTTAAACAAGAAAAGAAAAAAGGGCTTTTGCGccgg                      >  2:2080890/1‑148 (MQ=255)
                                                                                                                                       ccGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg               <  2:3088393/150‑1 (MQ=255)
                                                                                                                                        gggggTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc              <  2:3215342/149‑1 (MQ=255)
                                                                                                                                         ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  1:2547289/150‑1 (MQ=255)
                                                                                                                                           ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga           <  1:4120618/150‑1 (MQ=255)
                                                                                                                                                gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa      <  1:355881/150‑1 (MQ=255)
                                                                                                                                                gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa      <  1:355863/150‑1 (MQ=255)
                                                                                                                                                gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa      <  1:3257588/150‑1 (MQ=255)
                                                                                                                                                gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa      <  2:3553659/150‑1 (MQ=255)
                                                                                                                                                 ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa     <  2:511831/150‑1 (MQ=255)
                                                                                                                                                  gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg    <  2:1723538/150‑1 (MQ=255)
                                                                                                                                                    aTCCGCGCGCACGTGCAACCACACCCCAAACCAATCCCTCAAGGGAGAGGGGACCCAAACAGCACAAATTTTTAGAGATATAATCATAATTATTTTAAAAGAAGAAAACAAAAAACGGCTTTTTCGCCCGAATAATTAGGAGAAAAAGCa  >  2:3754665/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164433‑164730

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: