Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,743,163 T→C 5.1% E27G (GAG→GGG)  ribC ← riboflavin synthase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,743,1630TC 5.1% 240.2 / 2.9 79E27G (GAG→GGG) ribCriboflavin synthase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (34/41);  new base C (2/2);  total (36/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 = 9.31e-01

AAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCTTACCA  >  NC_000913/1743016‑1743311
                                                                                                                                                   |                                                                                                                                                    
aaaTCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTcc                                                                                                                                                    >  1:855547/1‑150 (MQ=255)
   tCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTccacc                                                                                                                                                 <  1:317453/150‑1 (MQ=255)
    cGCCAAGATTGGTACTGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTccacca                                                                                                                                                <  1:815168/150‑1 (MQ=255)
     gCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCccccac                                                                                                                                               >  1:1016186/1‑150 (MQ=255)
         aGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACATGc                                                                                                                                           <  1:1124865/150‑1 (MQ=255)
          gATTGGTAATGCGTACCGTTTCTTTCAACAGGTCAAAACTGACATGGTTGCCGTTAAATTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGGTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACATGct                                                                                                                                          <  2:1140596/150‑2 (MQ=255)
             tGGTAATGCGGAACGTTTCTTTCCTCAGCTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTAc                                                                                                                                       <  1:1104971/150‑1 (MQ=255)
              ggTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCATGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACg                                                                                                                                      >  1:728821/1‑150 (MQ=255)
              ggTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACg                                                                                                                                      <  2:544236/150‑1 (MQ=255)
                    gCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAAtt                                                                                                                                >  2:1079359/1‑150 (MQ=255)
                      gTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTg                                                                                                                              <  2:294538/150‑1 (MQ=255)
                         aCGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGtt                                                                                                                           <  2:547859/150‑1 (MQ=255)
                              tctttcATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCg                                                                                                                      <  1:1071613/150‑1 (MQ=255)
                              tctttcATCAGGTCAAAACTGACATGCGTGCCGGAAATTTCCGACACGGTCAGCCAGCAACCGTTATGCGCCACGGAAGCTCCGGGTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCg                                                                                                                      <  1:186615/150‑1 (MQ=255)
                                  tcatcaGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCaa                                                                                                                  <  2:907180/150‑1 (MQ=255)
                                   catcaGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAAt                                                                                                                 >  1:1021115/1‑150 (MQ=255)
                                     tcaGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCg                                                                                                               <  1:542340/150‑1 (MQ=255)
                                         gTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGacac                                                                                                           >  1:19629/1‑150 (MQ=255)
                                         gTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGacac                                                                                                           <  2:614141/150‑1 (MQ=255)
                                         gTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGacac                                                                                                           >  1:648417/1‑150 (MQ=255)
                                             aaaCTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGt                                                                                                       >  2:1220135/1‑150 (MQ=255)
                                             aaaCTGACATGGTTGCCGTTAATTTCCGTCACGGGCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGt                                                                                                       <  1:795006/150‑1 (MQ=255)
                                               aCTGACATGGTTGCCGTTAATTTCCCTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAAGTCCACCACATGCGGACGAAAAGTTGGTTTCTCGTCAATCGACTCCAGtat                                                                                                     >  2:585357/1‑148 (MQ=255)
                                                  gACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGc                                                                                                  >  2:430390/1‑150 (MQ=255)
                                                      tGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTg                                                                                              >  1:700822/1‑150 (MQ=255)
                                                       ggTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTACCTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGc                                                                                             >  1:791332/1‑150 (MQ=255)
                                                              gTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTAc                                                                                      <  1:1247011/150‑1 (MQ=255)
                                                               ttAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGATTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACa                                                                                     <  1:206272/150‑1 (MQ=255)
                                                                   tttCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGGGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAAtgc                                                                                 >  2:253457/1‑148 (MQ=255)
                                                                     tCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACATGcccc                                                                               >  1:509727/1‑150 (MQ=255)
                                                                     tCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATcccc                                                                               <  1:149289/150‑1 (MQ=255)
                                                                      ccGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCg                                                                              >  2:1081755/1‑150 (MQ=255)
                                                                      ccGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCACCg                                                                              >  2:47182/1‑150 (MQ=255)
                                                                       cGTCACGGTCAGGCAGCAACCGTTATGCGCCACTGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAATTTTGGTTTCTCGTCAATCGACACCAGTTTTGCTGTGCCCTGTACAATCCCCGt                                                                             <  1:1209056/150‑1 (MQ=255)
                                                                           aCGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTACCTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCgaacac                                                                         >  2:82158/1‑145 (MQ=255)
                                                                              ttCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATa                                                                      <  1:765585/149‑1 (MQ=255)
                                                                                    cagcaACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACa                                                                <  2:19629/150‑1 (MQ=255)
                                                                                      gcaACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAAc                                                              <  2:1274143/150‑1 (MQ=255)
                                                                                           cGTTCTGCGCCACGGCAGCACCGGTTTCCAGGCCGTCCAGGATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTAAAACTCCTg                                                         <  2:791251/150‑1 (MQ=255)
                                                                                           cGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGGCCAGCATGTGGTCGGGTAACTCCACCAAATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTg                                                         <  2:462885/150‑1 (MQ=255)
                                                                                             ttATGCGCCACGGAAGCACCGGGTTCCATGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCCACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGaa                                                       <  2:665249/150‑1 (MQ=255)
                                                                                              tATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGaaa                                                      <  1:431171/150‑1 (MQ=255)
                                                                                                   gcCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGt                                                 <  1:1267905/150‑1 (MQ=255)
                                                                                                    ccACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGtt                                                <  1:754961/150‑1 (MQ=255)
                                                                                                    ccACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACAAGCGTACGATTATTTGGTTTCTCGTCAATCGCCACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAGCCTACTTCCAGCTCCTGAAATCAGtt                                                >  2:816637/1‑150 (MQ=255)
                                                                                                     cACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTa                                               >  1:538774/1‑150 (MQ=255)
                                                                                                     cACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTa                                               <  2:1293904/150‑1 (MQ=255)
                                                                                                     cACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGAGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTa                                               <  2:777195/150‑1 (MQ=255)
                                                                                                     cACGGAAGCACCGGGTTCCAGGCCGTCCAGCATGGGGGCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTa                                               <  2:126363/150‑1 (MQ=255)
                                                                                                         gAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCCCCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGAc                                           >  2:739660/1‑150 (MQ=255)
                                                                                                         gAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAAGCAGTTAAGAc                                           >  1:550521/1‑150 (MQ=255)
                                                                                                         gAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGGACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGGAAACATACTTACAACACCGGAAATCAGTTAAGAc                                           >  2:865251/1‑150 (MQ=255)
                                                                                                          aaGCACCGGTTTCCAGGCCGGCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACa                                          <  1:1092086/150‑1 (MQ=255)
                                                                                                            gCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAGCTCCTGAAATCAGTTAAGACAtt                                        >  1:68009/1‑150 (MQ=255)
                                                                                                             cACCGGTTTCCAGGCCGTCCAGCATGTCGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACCACTCCTGAAATCAGTTAAGACATTc                                       <  1:63438/150‑1 (MQ=255)
                                                                                                              aCCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCt                                      >  1:764242/1‑150 (MQ=255)
                                                                                                                 ggTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCCCCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACACCTCCTGAAATCAGTTAAGACATTCTGtt                                   >  1:969565/1‑150 (MQ=255)
                                                                                                                  gTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTc                                  <  2:1250396/150‑1 (MQ=255)
                                                                                                                  gTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTc                                  <  2:468540/150‑1 (MQ=255)
                                                                                                                    ttCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAg                                <  2:1315670/150‑1 (MQ=255)
                                                                                                                    ttCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAg                                <  2:1235553/150‑1 (MQ=255)
                                                                                                                    ttCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGATTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAg                                <  2:1181710/150‑1 (MQ=255)
                                                                                                                     tCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTTCAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGc                               >  2:879144/1‑150 (MQ=255)
                                                                                                                      ccAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACCACTCCTGAAATCAGTTAAGACATTCTGTTCAGca                              >  2:905269/1‑150 (MQ=255)
                                                                                                                          gCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAAt                          >  2:68918/1‑150 (MQ=255)
                                                                                                                          gCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAAt                          <  2:677790/150‑1 (MQ=255)
                                                                                                                              tCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCa                      <  1:1025533/150‑1 (MQ=255)
                                                                                                                               tacgCATGTGGTGGAGTTACCCGACCACATGCGTTCGAAAATTTGGTTTCTCGTCAATCGCCACCAGTTTTGCGGTGCCCTCTCCAAACCCCGTAAACCTACTTACCACTCCTGAAATCAGTTAAGACCTTCTGTTCAGCACAATAGCAg                     >  1:380427/4‑150 (MQ=255)
                                                                                                                                cAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCGATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAgg                    >  1:881798/1‑150 (MQ=255)
                                                                                                                                 aGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGt                   <  2:1145966/150‑1 (MQ=255)
                                                                                                                                 aGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGt                   >  2:1130629/1‑150 (MQ=255)
                                                                                                                                  gCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTg                  <  1:216764/150‑1 (MQ=255)
                                                                                                                                   cATGTGGTGGAGTTACCCGACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTgg                 >  2:540798/1‑150 (MQ=255)
                                                                                                                                    aTGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGa                <  2:648417/150‑1 (MQ=255)
                                                                                                                                    aTGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGa                >  1:399722/1‑150 (MQ=255)
                                                                                                                                        ggTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTAAGCACAATAGAAGGTGGAAAAc            >  2:1312857/1‑150 (MQ=255)
                                                                                                                                        ggTCGGGTAACTCCACCACATGCGTACGAAAATATGGTTGCTCGTCAATCGACACCAGTCTAGCGGTGCCCTGTACAAACCCCGTAAACATACTTACAACTCCTGAAATCAGTGAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAAc            <  2:219691/150‑1 (MQ=255)
                                                                                                                                             ggTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAAACGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGATAAGACATTCTGTTCAGCACAATAGCAGGTGGCAAACGCCCt       >  2:1214777/1‑150 (MQ=255)
                                                                                                                                                  cTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAACCGCCCTGACCa  >  2:1284013/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAATCGCCAAGATTGGTAATGCGTAACGTTTCTTTCATCAGGTCAAAACTGACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAGTTAAGACATTCTGTTCAGCACAATAGCAGGTGGAAAACGCCCTTACCA  >  NC_000913/1743016‑1743311

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