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

ACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCT  >  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)
                taTTGGTAATGCGTAACGTTTCTTTCATCAGGGCAAAACTGACCTGGTGGCCGTTAATTTCCGTCACGGGCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑g                                                                                                                                     <  1:347426/149‑1 (MQ=255)
                gATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑g                                                                                                                                     <  2:560430/150‑1 (MQ=255)
                gATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCACCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑g                                                                                                                                     <  1:905430/150‑1 (MQ=255)
                caTTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACATGct                                                                                                                                      <  1:131633/149‑2 (MQ=255)
                   tGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTAc                                                                                                                                  >  1:1651322/1‑150 (MQ=255)
                    ggTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACg                                                                                                                                 >  1:948964/1‑150 (MQ=255)
                    ggTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACg                                                                                                                                 >  2:883146/1‑150 (MQ=255)
                    ggTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGA‑GTACg                                                                                                                                 >  2:381029/1‑150 (MQ=255)
                      tAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGaa                                                                                                                               >  2:1664335/1‑150 (MQ=255)
                            gTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTg                                                                                                                         <  1:654696/150‑1 (MQ=255)
                            gTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTg                                                                                                                         <  2:576531/150‑1 (MQ=255)
                              aaCGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAAATTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGt                                                                                                                       <  1:30108/150‑1 (MQ=255)
                              aaCGTTTCTTTCATCAGGTCAAAACTCACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGt                                                                                                                       <  1:1515933/150‑1 (MQ=255)
                                 gTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTtc                                                                                                                    <  1:426380/150‑1 (MQ=255)
                                    tctttcATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCg                                                                                                                 <  1:601381/150‑1 (MQ=255)
                                    tctttcATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCg                                                                                                                 <  2:516877/150‑1 (MQ=255)
                                     ctttcATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGt                                                                                                                <  1:1104076/150‑1 (MQ=255)
                                     ctttcATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGt                                                                                                                >  1:596226/1‑150 (MQ=255)
                                      tttcATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTc                                                                                                               <  2:1036276/150‑1 (MQ=255)
                                         catcaGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAAt                                                                                                            >  1:606439/1‑150 (MQ=255)
                                              ggTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGaca                                                                                                       >  2:727011/1‑150 (MQ=255)
                                              ggTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGaca                                                                                                       >  1:671520/1‑150 (MQ=255)
                                              ggTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGaca                                                                                                       >  1:789524/1‑150 (MQ=255)
                                               gTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGacac                                                                                                      >  2:484489/1‑150 (MQ=255)
                                                   aaaCTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGt                                                                                                  <  2:1453091/150‑1 (MQ=255)
                                                   aaaCTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGt                                                                                                  >  1:1667524/1‑150 (MQ=255)
                                                    aaCTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCGACCGTTATGCGCAACGGAAGGAACGTTTTCCAGGGCGTCCAGCATGTGGTCGGGTAACTACGCCACATGC‑GTACGAAAATTTGGTTTCTCTTCAATCGACGCCTGtt                                                                                                 >  2:262041/1‑150 (MQ=255)
                                                      cTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGtttt                                                                                               <  2:162543/150‑1 (MQ=255)
                                                          cATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCgg                                                                                           <  1:883146/150‑1 (MQ=255)
                                                            tGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTg                                                                                         <  1:370739/150‑1 (MQ=255)
                                                              gTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGcc                                                                                       >  1:898397/1‑150 (MQ=255)
                                                                tGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGCCACCAGTTTTGCGGTGCCCt                                                                                     >  2:1173523/1‑150 (MQ=255)
                                                                   cGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACTGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCTGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTa                                                                                  <  2:583926/150‑1 (MQ=255)
                                                                   cGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTa                                                                                  >  1:434157/1‑150 (MQ=255)
                                                                     ttAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACa                                                                                >  2:1242186/1‑150 (MQ=255)
                                                                          ttCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATccc                                                                           <  2:1363392/150‑1 (MQ=255)
                                                                           tCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATcccc                                                                          >  2:1553446/1‑150 (MQ=255)
                                                                            ccGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCg                                                                         <  2:1437592/150‑1 (MQ=255)
                                                                             cGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGt                                                                        >  2:1176882/1‑150 (MQ=255)
                                                                               ccACGGTCATGCAGCAACCGTTAAGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCAGTACAATCCCCGTaa                                                                      <  1:1470520/149‑1 (MQ=255)
                                                                                     tCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACTTAc                                                                <  1:727011/150‑1 (MQ=255)
                                                                                        ggCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTa                                                             >  2:1691147/1‑150 (MQ=255)
                                                                                           agcaACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACaa                                                          >  2:312339/1‑150 (MQ=255)
                                                                                             caACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACt                                                        >  1:1520612/1‑150 (MQ=255)
                                                                                              aaCCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACACTCCCCGTAAACATACTTACAACTc                                                       >  1:1283496/1‑150 (MQ=255)
                                                                                              aaCCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTc                                                       >  2:674073/1‑150 (MQ=255)
                                                                                                 cGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTg                                                    <  2:463303/150‑1 (MQ=255)
                                                                                                      tGCGCCACGGAAGCACCTGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGA‑GTACGAAAATTTGGTTTCTCGTCAATAGACACCAGGTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATc                                               <  1:393116/150‑1 (MQ=255)
                                                                                                      tGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATc                                               <  2:746736/150‑1 (MQ=255)
                                                                                                        cgcCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAg                                             <  2:948964/150‑1 (MQ=255)
                                                                                                             cGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAg                                        <  2:977532/150‑1 (MQ=255)
                                                                                                                     ccGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTg                                <  2:898397/150‑1 (MQ=255)
                                                                                                                     ccGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGGACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTg                                >  1:455513/1‑150 (MQ=255)
                                                                                                                       ggTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCACGTAAACATACATACAACTCCTGAAATCAGTTATGAAATTCTGtt                              >  2:647485/1‑150 (MQ=255)
                                                                                                                          ttCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAg                           <  2:671520/150‑1 (MQ=255)
                                                                                                                            ccAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGca                         >  1:582212/1‑150 (MQ=255)
                                                                                                                             cAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGCCATTCTGTTCAGctc                        >  2:246965/1‑148 (MQ=255)
                                                                                                                              aGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGcaca                       <  2:1320443/150‑1 (MQ=255)
                                                                                                                                 ttcgtacgCATGTGGTGGAGTTACCCGACCACATGCTGGACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAAt                     >  1:397359/8‑150 (MQ=37)
                                                                                                                                  cGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACTAAAATTTGGTGTCTCGTCAATCGCCACCAGTTGTGCTGTGCCCTGTACAATCCCCGTAAACCTACTTCCACCGCCTGAAATCAGGTAAGACATTCTTTTCAGCACACTat                   >  1:174336/1‑149 (MQ=255)
                                                                                                                                  cGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAg                   >  1:1395373/1‑150 (MQ=255)
                                                                                                                                    gtacgCATGTGGTGGAGTTACCCGACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCa                 >  2:365828/5‑150 (MQ=255)
                                                                                                                                     ccAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGCCATTCTGTTCAGCACAATAGCAg                >  2:1188275/1‑150 (MQ=255)
                                                                                                                                     ccAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAg                <  2:1313548/150‑1 (MQ=255)
                                                                                                                                     ccAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAg                <  1:312339/150‑1 (MQ=255)
                                                                                                                                      cAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAgg               >  1:1367161/1‑150 (MQ=255)
                                                                                                                                       cgCATGTGGTGGAGTTACCCGACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACTCCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAAGTCCTTAAATCAGTTAAGACATTCTGTACAGCACAATAGCAGGt              >  1:1215610/2‑150 (MQ=255)
                                                                                                                                         cATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGGTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGGTCAGCACAATAGCAGGTgg            >  2:900153/1‑150 (MQ=255)
                                                                                                                                          aTGTGGTGGAGTTACCCGACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGCAATCAGTTAAGCCATTCTGTTCAGCACAATAGCAGGTGGa           >  1:1384146/1‑150 (MQ=255)
                                                                                                                                            gtgGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAAGATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGaaa         >  1:1646608/1‑150 (MQ=255)
                                                                                                                                                  gggTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGccc   >  1:498634/1‑150 (MQ=255)
                                                                                                                                                   ggTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGAAAAACGCCCt  >  1:396309/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
ACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCT  >  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: ‑