Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,743,157 G→A 9.1% P29L (CCC→CTC)  ribC ← riboflavin synthase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,743,1570GA 9.1% 258.4 / 11.6 88P29L (CCC→CTC) ribCriboflavin synthase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (42/38);  new base A (4/4);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

ACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCT  >  NC_000913/1743010‑1743306
                                                                                                                                                   |                                                                                                                                                      
aCTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGt                                                                                                                                                      <  1:1001411/150‑1 (MQ=255)
 cTTTTAATTCTCCAAGTTTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAATCTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCATCCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTa                                                                                                                                                     <  2:676940/150‑1 (MQ=255)
 cTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTa                                                                                                                                                     >  1:1036276/1‑150 (MQ=255)
   tttAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCACCCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCCGCATGTGGTCGGGTAAc                                                                                                                                                   >  2:732236/1‑150 (MQ=255)
    ttAAATCGCCAAGATTGGTAATGCGTATCGTTTCTTTCATCAGGTCACAACTGACATGGGTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGGCGGGTAACt                                                                                                                                                  <  1:46616/150‑1 (MQ=255)
    ttAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACt                                                                                                                                                  >  1:59625/1‑150 (MQ=255)
    ttAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACt                                                                                                                                                  <  2:566906/150‑1 (MQ=255)
     tAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACTGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACta                                                                                                                                                 >  1:1721078/1‑149 (MQ=255)
      aaaTCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTcc                                                                                                                                                >  1:1151611/1‑150 (MQ=255)
        aTCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGGCCAGCATGTGGTCGGGTAACTccac                                                                                                                                              <  2:1192320/150‑1 (MQ=255)
           gCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCccccac                                                                                                                                           >  1:283201/1‑150 (MQ=255)
            ccAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACcaca                                                                                                                                          >  1:1320443/1‑150 (MQ=255)
            ccAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGGCAAAACTGACATGGTTGCCGTTAATTTCCGTCCCGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACcaca                                                                                                                                          <  2:780671/150‑1 (MQ=255)
             cAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACAt                                                                                                                                         <  1:607228/150‑1 (MQ=255)
             aaaGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACAt                                                                                                                                         <  1:581305/149‑1 (MQ=255)
              aaGATTGGTACTGCGTAACGTTTCTTTCATCCTGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTGTCCTGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATg                                                                                                                                        <  1:1084315/150‑1 (MQ=255)
               aGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCAACCACATGc                                                                                                                                       >  1:977532/1‑150 (MQ=255)
                taTTGGTAATGCGTAACGTTTCTTTGATCAGGTCAAAACTGACATGGATGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGGGGAGTTACCCGACCACATGct                                                                                                                                      <  1:43731/149‑2 (MQ=255)
                taTTGGTAATGCGTAACGTTTCTTTCATCAGGGCAAAACTGACCTGGTGGCCGTTAATTTCCGTCACGGGCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCg                                                                                                                                     <  1:347426/149‑1 (MQ=255)
                gATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCg                                                                                                                                     <  2:560430/150‑1 (MQ=255)
                gATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCACCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCg                                                                                                                                     <  1:905430/150‑1 (MQ=255)
                caTTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACATGct                                                                                                                                      <  1:131633/149‑2 (MQ=255)
                   tGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTAc                                                                                                                                  >  1:1651322/1‑150 (MQ=255)
                    ggTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACg                                                                                                                                 >  1:948964/1‑150 (MQ=255)
                    ggTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACg                                                                                                                                 >  2:883146/1‑150 (MQ=255)
                    ggTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGAGTACg                                                                                                                                 >  2:381029/1‑150 (MQ=255)
                      tAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGaa                                                                                                                               >  2:1664335/1‑150 (MQ=255)
                            gTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTg                                                                                                                         <  1:654696/150‑1 (MQ=255)
                            gTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTg                                                                                                                         <  2:576531/150‑1 (MQ=255)
                              aaCGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAAATTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGt                                                                                                                       <  1:30108/150‑1 (MQ=255)
                              aaCGTTTCTTTCATCAGGTCAAAACTCACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGt                                                                                                                       <  1:1515933/150‑1 (MQ=255)
                                 gTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTtc                                                                                                                    <  1:426380/150‑1 (MQ=255)
                                    tctttcATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCg                                                                                                                 <  1:601381/150‑1 (MQ=255)
                                    tctttcATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCg                                                                                                                 <  2:516877/150‑1 (MQ=255)
                                     ctttcATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGt                                                                                                                <  1:1104076/150‑1 (MQ=255)
                                     ctttcATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGt                                                                                                                >  1:596226/1‑150 (MQ=255)
                                      tttcATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTc                                                                                                               <  2:1036276/150‑1 (MQ=255)
                                         catcaGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAAt                                                                                                            >  1:606439/1‑150 (MQ=255)
                                              ggTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGaca                                                                                                       >  2:727011/1‑150 (MQ=255)
                                              ggTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGaca                                                                                                       >  1:671520/1‑150 (MQ=255)
                                              ggTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGaca                                                                                                       >  1:789524/1‑150 (MQ=255)
                                               gTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGacac                                                                                                      >  2:484489/1‑150 (MQ=255)
                                                   aaaCTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGt                                                                                                  <  2:1453091/150‑1 (MQ=255)
                                                   aaaCTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGt                                                                                                  >  1:1667524/1‑150 (MQ=255)
                                                    aaCTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCGACCGTTATGCGCAACGGAAGGAACGTTTTCCAGGGCGTCCAGCATGTGGTCGGGTAACTACGCCACATGCGTACGAAAATTTGGTTTCTCTTCAATCGACGCCTGtt                                                                                                 >  2:262041/1‑150 (MQ=255)
                                                      cTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGtttt                                                                                               <  2:162543/150‑1 (MQ=255)
                                                          cATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCgg                                                                                           <  1:883146/150‑1 (MQ=255)
                                                            tGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTg                                                                                         <  1:370739/150‑1 (MQ=255)
                                                              gTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGcc                                                                                       >  1:898397/1‑150 (MQ=255)
                                                                tGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGCCACCAGTTTTGCGGTGCCCt                                                                                     >  2:1173523/1‑150 (MQ=255)
                                                                   cGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACTGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCTGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTa                                                                                  <  2:583926/150‑1 (MQ=255)
                                                                   cGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTa                                                                                  >  1:434157/1‑150 (MQ=255)
                                                                     ttAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACa                                                                                >  2:1242186/1‑150 (MQ=255)
                                                                          ttCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATccc                                                                           <  2:1363392/150‑1 (MQ=255)
                                                                           tCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATcccc                                                                          >  2:1553446/1‑150 (MQ=255)
                                                                            ccGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCg                                                                         <  2:1437592/150‑1 (MQ=255)
                                                                             cGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGt                                                                        >  2:1176882/1‑150 (MQ=255)
                                                                               ccACGGTCATGCAGCAACCGTTAAGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCAGTACAATCCCCGTaa                                                                      <  1:1470520/149‑1 (MQ=255)
                                                                                     tCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACTTAc                                                                <  1:727011/150‑1 (MQ=255)
                                                                                        ggCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTa                                                             >  2:1691147/1‑150 (MQ=255)
                                                                                           agcaACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACaa                                                          >  2:312339/1‑150 (MQ=255)
                                                                                             caACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACt                                                        >  1:1520612/1‑150 (MQ=255)
                                                                                              aaCCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACACTCCCCGTAAACATACTTACAACTc                                                       >  1:1283496/1‑150 (MQ=255)
                                                                                              aaCCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTc                                                       >  2:674073/1‑150 (MQ=255)
                                                                                                 cGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTg                                                    <  2:463303/150‑1 (MQ=255)
                                                                                                      tGCGCCACGGAAGCACCTGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGAGTACGAAAATTTGGTTTCTCGTCAATAGACACCAGGTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATc                                               <  1:393116/150‑1 (MQ=255)
                                                                                                      tGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATc                                               <  2:746736/150‑1 (MQ=255)
                                                                                                        cgcCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAg                                             <  2:948964/150‑1 (MQ=255)
                                                                                                             cGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAg                                        <  2:977532/150‑1 (MQ=255)
                                                                                                                     ccGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTg                                <  2:898397/150‑1 (MQ=255)
                                                                                                                     ccGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGGACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTg                                >  1:455513/1‑150 (MQ=255)
                                                                                                                       ggTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCACGTAAACATACATACAACTCCTGAAATCAGTTATGAAATTCTGtt                              >  2:647485/1‑150 (MQ=255)
                                                                                                                          ttCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAg                           <  2:671520/150‑1 (MQ=255)
                                                                                                                            ccAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGca                         >  1:582212/1‑150 (MQ=255)
                                                                                                                             cAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGCCATTCTGTTCAGctc                        >  2:246965/1‑148 (MQ=255)
                                                                                                                              aGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGcaca                       <  2:1320443/150‑1 (MQ=255)
                                                                                                                                 ttcgtacgCATGTGGTGGAGTTACCCGACCACATGCTGGACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAAt                     >  1:397359/8‑150 (MQ=37)
                                                                                                                                  cGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACTAAAATTTGGTGTCTCGTCAATCGCCACCAGTTGTGCTGTGCCCTGTACAATCCCCGTAAACCTACTTCCACCGCCTGAAATCAGGTAAGACATTCTTTTCAGCACACTat                   >  1:174336/1‑149 (MQ=255)
                                                                                                                                  cGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAg                   >  1:1395373/1‑150 (MQ=255)
                                                                                                                                    gtacgCATGTGGTGGAGTTACCCGACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCa                 >  2:365828/5‑150 (MQ=255)
                                                                                                                                     ccAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGCCATTCTGTTCAGCACAATAGCAg                >  2:1188275/1‑150 (MQ=255)
                                                                                                                                     ccAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAg                <  2:1313548/150‑1 (MQ=255)
                                                                                                                                     ccAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAg                <  1:312339/150‑1 (MQ=255)
                                                                                                                                      cAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAgg               >  1:1367161/1‑150 (MQ=255)
                                                                                                                                       cgCATGTGGTGGAGTTACCCGACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACTCCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAAGTCCTTAAATCAGTTAAGACATTCTGTACAGCACAATAGCAGGt              >  1:1215610/2‑150 (MQ=255)
                                                                                                                                         cATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGGTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGGTCAGCACAATAGCAGGTgg            >  2:900153/1‑150 (MQ=255)
                                                                                                                                          aTGTGGTGGAGTTACCCGACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGCAATCAGTTAAGCCATTCTGTTCAGCACAATAGCAGGTGGa           >  1:1384146/1‑150 (MQ=255)
                                                                                                                                            gtgGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAAGATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGaaa         >  1:1646608/1‑150 (MQ=255)
                                                                                                                                                  gggTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGccc   >  1:498634/1‑150 (MQ=255)
                                                                                                                                                   ggTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGAAAAACGCCCt  >  1:396309/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
ACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCT  >  NC_000913/1743010‑1743306

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