Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,743,155 C→G 9.0% D30H (GAC→CAC)  ribC ← riboflavin synthase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,743,1550CG 9.0% 285.5 / 15.9 89D30H (GAC→CAC) ribCriboflavin synthase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (42/39);  new base G (4/4);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

CCCCACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAA  >  NC_000913/1743006‑1743299
                                                                                                                                                     |                                                                                                                                                 
ccccACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTc                                                                                                                                                   >  1:463303/1‑150 (MQ=255)
ccccACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTc                                                                                                                                                   >  1:373448/1‑150 (MQ=255)
ccccACTTTTAAATCGCCAAGATTGGTAATGCGTAACGATTCTTTCATCAGGTCAAAACTGACATGGTTGCAGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTc                                                                                                                                                   >  1:884897/1‑150 (MQ=255)
 cccACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCg                                                                                                                                                  <  2:758302/150‑1 (MQ=255)
    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                                                                                                                          >  2:883146/1‑150 (MQ=255)
                        ggTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACg                                                                                                                          >  1:948964/1‑150 (MQ=255)
                        ggTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGA‑GTACg                                                                                                                          >  2:381029/1‑150 (MQ=255)
                          tAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGaa                                                                                                                        >  2:1664335/1‑150 (MQ=255)
                                gTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTg                                                                                                                  <  2:576531/150‑1 (MQ=255)
                                gTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTg                                                                                                                  <  1:654696/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                                                                                                          <  2:516877/150‑1 (MQ=255)
                                        tctttcATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCg                                                                                                          <  1:601381/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                                                                                                >  1:789524/1‑150 (MQ=255)
                                                  ggTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGaca                                                                                                >  1:671520/1‑150 (MQ=255)
                                                  ggTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGaca                                                                                                >  2:727011/1‑150 (MQ=255)
                                                   gTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGacac                                                                                               >  2:484489/1‑150 (MQ=255)
                                                       aaaCTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGt                                                                                           >  1:1667524/1‑150 (MQ=255)
                                                       aaaCTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGt                                                                                           <  2:1453091/150‑1 (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)
                                                                                                                                                     |                                                                                                                                                 
CCCCACTTTTAAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGC‑GTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAA  >  NC_000913/1743006‑1743299

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: ‑