Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,377,770 A→C 5.0% intergenic (‑350/+45) sspA ← / ← rpsI stringent starvation protein A/30S ribosomal subunit protein S9

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,377,7700AC 5.0% 375.4 / 9.7 120intergenic (‑350/+45)sspA/rpsIstringent starvation protein A/30S ribosomal subunit protein S9
Reads supporting (aligned to +/- strand):  ref base A (58/56);  new base C (4/2);  total (62/58)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 122 total aligned reads displayed.

AAAAAGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTC  >  NC_000913/3377623‑3377919
                                                                                                                                                   |                                                                                                                                                     
aaaaaGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAg                                                                                                                                                     >  1:1058295/1‑150 (MQ=255)
   aaGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCgg                                                                                                                                                  >  1:916964/1‑150 (MQ=255)
    agagATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCggg                                                                                                                                                 <  1:566569/150‑1 (MQ=255)
    agagATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCggg                                                                                                                                                 <  1:551550/150‑1 (MQ=255)
      agaTGCGCTTTAGTCTGAAATAGTTGAGTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCATCTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGtt                                                                                                                                               <  2:1907100/150‑1 (MQ=255)
       gaTGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGttt                                                                                                                                              >  2:868576/1‑150 (MQ=255)
         tGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGttttt                                                                                                                                            >  1:1090287/1‑150 (MQ=255)
          gcgcTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGtttttt                                                                                                                                           <  1:343110/150‑1 (MQ=255)
           cgcTTTAGTCTGAAATAGTTGACTTTGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTc                                                                                                                                          <  1:1542238/150‑1 (MQ=255)
           cgcTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCTTCGGCGGGttttttt                                                                                                                                          <  2:718230/150‑2 (MQ=255)
            gcTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCTTCGGCGGGttttttta                                                                                                                                         <  1:316774/150‑3 (MQ=255)
                tAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGaaaa                                                                                                                                     <  2:916964/150‑1 (MQ=255)
                  gTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAAtt                                                                                                                                   >  1:1850669/1‑150 (MQ=255)
                    cTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGt                                                                                                                                 >  2:1625263/1‑150 (MQ=255)
                    cTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGt                                                                                                                                 >  1:952617/1‑150 (MQ=255)
                    cTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGt                                                                                                                                 >  2:722993/1‑150 (MQ=255)
                     tGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGtt                                                                                                                                >  2:425779/1‑150 (MQ=255)
                      gAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGACGAAGCGGGTTTTTTCGAAAATTGttt                                                                                                                               >  2:813258/1‑150 (MQ=255)
                      gAAATAGTTGACTTAGTCCCTTATTGGCCATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGttt                                                                                                                               >  2:1090275/1‑150 (MQ=255)
                        aaTAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTc                                                                                                                             <  1:907569/150‑1 (MQ=255)
                             ttGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCg                                                                                                                        <  2:1090287/150‑1 (MQ=255)
                               gACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAATAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGa                                                                                                                      <  2:924043/150‑1 (MQ=255)
                                aCTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAg                                                                                                                     <  2:1782970/150‑1 (MQ=255)
                                       cccTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGc                                                                                                              >  1:1544527/1‑150 (MQ=255)
                                       cccTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGc                                                                                                              <  2:1278757/150‑1 (MQ=255)
                                       cccTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGc                                                                                                              >  1:274612/1‑150 (MQ=255)
                                             ttGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTa                                                                                                        >  2:1924942/1‑150 (MQ=255)
                                               ggCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTACTGCCGGAGCAGAAGCCAATTAAc                                                                                                      >  2:1004431/1‑150 (MQ=255)
                                                gCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACg                                                                                                     <  2:274612/150‑1 (MQ=255)
                                                    tgtgGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTg                                                                                                 <  1:813258/150‑1 (MQ=255)
                                                     gtgGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTgg                                                                                                >  2:1174234/1‑150 (MQ=255)
                                                     gtgGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTgg                                                                                                <  2:1806301/150‑1 (MQ=255)
                                                     gtgGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTgg                                                                                                <  1:167406/150‑1 (MQ=255)
                                                      tgGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGga                                                                                               <  2:770680/150‑1 (MQ=255)
                                                       ggTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGgag                                                                                              >  2:314016/1‑150 (MQ=255)
                                                       ggTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGgag                                                                                              >  1:893573/1‑150 (MQ=255)
                                                       ggTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACATTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGgag                                                                                              >  1:210272/1‑150 (MQ=255)
                                                            ttGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTg                                                                                         >  2:1995089/1‑150 (MQ=255)
                                                             tGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGc                                                                                        >  1:1619822/1‑150 (MQ=255)
                                                                 ttACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgac                                                                                    <  1:1174234/150‑1 (MQ=255)
                                                                   aCCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacga                                                                                  >  2:711548/1‑150 (MQ=255)
                                                                    ccTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacgac                                                                                 <  1:722804/150‑1 (MQ=255)
                                                                    ccTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacgac                                                                                 >  1:1893331/1‑150 (MQ=255)
                                                                          tgtcAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTg                                                                           <  2:1058295/150‑1 (MQ=255)
                                                                          tgtcAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTg                                                                           >  1:163292/1‑150 (MQ=255)
                                                                          tgtcAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTg                                                                           >  2:1011586/1‑150 (MQ=255)
                                                                           gtcAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGc                                                                          >  2:132629/1‑150 (MQ=255)
                                                                            tcAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCt                                                                         <  2:424713/150‑1 (MQ=255)
                                                                             cAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCtt                                                                        <  1:722993/150‑1 (MQ=255)
                                                                               ggTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTa                                                                      >  2:1295658/1‑150 (MQ=255)
                                                                                  ggCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGc                                                                   <  1:132056/150‑1 (MQ=255)
                                                                                  ggCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGc                                                                   <  2:773977/150‑1 (MQ=255)
                                                                                   gcagcaAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCa                                                                  <  2:893573/150‑1 (MQ=255)
                                                                                    cagcaAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCAGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAg                                                                 >  2:163814/1‑150 (MQ=255)
                                                                                       caAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGAcc                                                              >  1:1822931/1‑150 (MQ=255)
                                                                                           aaGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGCCt                                                          >  2:1935897/1‑150 (MQ=255)
                                                                                           aaGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACt                                                          <  2:1639285/150‑1 (MQ=255)
                                                                                           aaGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTACGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGa                                                            <  1:1967686/148‑1 (MQ=255)
                                                                                             gCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACttt                                                        >  1:1559811/1‑150 (MQ=255)
                                                                                                aCTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCtt                                                     <  2:218671/150‑1 (MQ=255)
                                                                                                 ttttCCAGTTTTTAGTCTGATTGAGATTTTAGGTATAAAAAAACCCGCCGAATCGGGTTTTTTCGAAAATTGTTTTCTGCGGGAGCAGAATCCAATTAACGTTTGGAGAACTGGGGACGACGACGTGCTTTACGCTGACCGACTTTCTTa                                                    <  2:888722/149‑1 (MQ=255)
                                                                                                     ccAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGtt                                                <  1:868576/150‑1 (MQ=255)
                                                                                                      cAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTc                                               >  1:1510228/1‑150 (MQ=255)
                                                                                                        gTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCaa                                             <  1:314016/150‑1 (MQ=255)
                                                                                                          ttttACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAAcc                                           <  1:1011586/150‑1 (MQ=255)
                                                                                                          ttttACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAAcc                                           <  1:1295658/150‑1 (MQ=255)
                                                                                                          ttttACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAAcc                                           <  1:711548/150‑1 (MQ=255)
                                                                                                              aCGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGAc                                       >  1:477620/1‑150 (MQ=255)
                                                                                                               cGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACg                                      >  1:1524954/1‑150 (MQ=255)
                                                                                                                gCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGa                                     >  2:1087602/1‑150 (MQ=255)
                                                                                                                     ttCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTc                                >  1:75069/1‑150 (MQ=255)
                                                                                                                      tCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCa                               <  2:1544527/150‑1 (MQ=255)
                                                                                                                      tCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCa                               >  2:707729/1‑150 (MQ=255)
                                                                                                                       cAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCAc                              <  2:457676/150‑1 (MQ=255)
                                                                                                                        aGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAAc                                            <  1:1405093/135‑1 (MQ=255)
                                                                                                                        aGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAAc                                            >  2:1405093/1‑135 (MQ=255)
                                                                                                                        aGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACg                             <  2:163292/150‑1 (MQ=255)
                                                                                                                        aGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACg                             <  1:1040099/150‑1 (MQ=255)
                                                                                                                        aGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACg                             <  1:221385/150‑1 (MQ=255)
                                                                                                                         gATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCTGACCGACTTTCTTACGTTCAACCTGACGAGCGTCTCGa                            >  1:1300986/1‑150 (MQ=255)
                                                                                                                         gATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGGTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACTTTCAACCTGACGAGCGTCACGa                            <  2:1968573/150‑1 (MQ=255)
                                                                                                                          gttttAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTc                                >  2:1431614/2‑145 (MQ=255)
                                                                                                                          aTTTTAGCTATAAAAAAACCCGCCTAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAg                           <  2:952617/150‑1 (MQ=255)
                                                                                                                                 tttcgaaaaaaCCCGCTTCGGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGaa                    >  2:1776236/6‑150 (MQ=255)
                                                                                                                                    taAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTc                                <  1:983573/135‑1 (MQ=255)
                                                                                                                                    cgaaaaaaCCCGCTTCGGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGAACGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGcc                 >  1:347441/3‑150 (MQ=255)
                                                                                                                                     aaaaaaaCCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCa                >  2:1168320/1‑150 (MQ=255)
                                                                                                                                      aaaaaaCCCGCTTCGGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTTCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAg               >  2:1935973/1‑150 (MQ=255)
                                                                                                                                      aaaaaaCCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGcc                 <  1:1591950/148‑1 (MQ=255)
                                                                                                                                      aaaaaaCCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAg               <  2:159513/150‑1 (MQ=255)
                                                                                                                                       aaaaaCCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGc              <  2:1822931/150‑1 (MQ=255)
                                                                                                                                       aaaaaCCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTCCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGc              <  1:1004431/150‑1 (MQ=255)
                                                                                                                                        aaaaCCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCt             <  2:1383129/150‑1 (MQ=255)
                                                                                                                                        aaaaCCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCt             >  1:1216217/1‑150 (MQ=255)
                                                                                                                                           aCCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACTCAGACCAACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTa          >  2:1193321/1‑150 (MQ=255)
                                                                                                                                            cccGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTAc         >  2:583132/1‑150 (MQ=255)
                                                                                                                                            cccGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTAc         >  1:1033041/1‑150 (MQ=255)
                                                                                                                                              cgccgAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGc       <  1:1090275/150‑1 (MQ=255)
                                                                                                                                                 cgAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGt    <  2:1898902/150‑1 (MQ=255)
                                                                                                                                                   aaGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTc  <  1:1087602/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AAAAAGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTC  >  NC_000913/3377623‑3377919

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: