Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 21,506 G→T A34S (GCG→TCG)  ribF → bifunctional riboflavin kinase/FAD synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091321,5060GT98.7% 231.7 / NA 77A34S (GCG→TCG) ribFbifunctional riboflavin kinase/FAD synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (36/40);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTT  >  NC_000913/21360‑21647
                                                                                                                                                  |                                                                                                                                             
cGCTGTACAAGGTATACTCGGCCGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATCCGCGGCAGACAGAATCTCAGCCAGGCCCCGCAAGAAGGGGGTGGGCTGACTATTGGGAATTTCGACGGCGTGCATCGCGGTCATcgctcg                                                                                                                                             >  1:730121/1‑149 (MQ=255)
 gCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCAtcgctcgc                                                                                                                                            <  1:1438915/149‑1 (MQ=255)
 gCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCAtcgctcgc                                                                                                                                            <  2:355156/149‑1 (MQ=255)
      aCAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATTCGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTa                                                                                                                                       <  1:972177/149‑1 (MQ=255)
      aCAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTa                                                                                                                                       <  1:811943/149‑1 (MQ=255)
      aCAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTa                                                                                                                                       <  1:93661/149‑1 (MQ=255)
        aaGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGTCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGCCGGCGTGCATCGCGGTCATCGCTCGCTGTTACa                                                                                                                                     >  1:1105612/1‑149 (MQ=255)
          ggTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAgg                                                                                                                                   >  1:1185658/1‑149 (MQ=255)
          ggTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAgg                                                                                                                                   >  1:334886/1‑149 (MQ=255)
            tataCTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGc                                                                                                                                 <  2:1282316/149‑1 (MQ=255)
            tataCTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGc                                                                                                                                 <  1:334457/149‑1 (MQ=255)
            tataCTCGGACGATTTTCAATGTTTTGAGCCAGACATGAAGCTGATCCGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTTCTGACTATTTGTAATTTCGACGGCGTGCATCTCGGTCATCGCTCGCTGTTACAGTGc                                                                                                                                 <  1:340159/149‑1 (MQ=255)
                 tCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCa                                                                                                                            >  1:780182/1‑149 (MQ=255)
                     aCGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGGGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGgaa                                                                                                                        <  1:259819/149‑1 (MQ=255)
                         ttttGAGTGTTTTGCGCCAGACATGAAGCGGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGAGTATTGTTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTCTTACCGGGCGTGCAGGaagaag                                                                                                                    <  2:883663/149‑1 (MQ=255)
                          tttCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCTTACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAgg                                                                                                                   >  1:380573/1‑149 (MQ=255)
                          tttCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAgg                                                                                                                   >  1:878087/1‑149 (MQ=255)
                            tCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGgc                                                                                                                 <  1:1649845/149‑1 (MQ=255)
                             cACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGgcg                                                                                                                <  1:1238633/149‑1 (MQ=255)
                               cTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGgcgca                                                                                                              >  1:1158040/1‑149 (MQ=255)
                                ttgtttGAGCCTGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGGGTGTGCTGACTATTGGTAATTTCGACTCCGTGCATCGCGGTCATCGCTCGCTGTTACAGTGCTTGCATGACGAAGGgcgcaa                                                                                                             <  2:730121/146‑1 (MQ=255)
                                tGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGgcgcaa                                                                                                             >  1:1362419/1‑149 (MQ=255)
                                      gAGCGACACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGCCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCGTGCAGGACGAAGGgcgcaagcgcaa                                                                                                       <  1:860729/149‑1 (MQ=255)
                                      gAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGgcgcaagcgcaa                                                                                                       <  1:21868/149‑1 (MQ=255)
                                       aGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAAc                                                                                                      >  1:176677/1‑149 (MQ=255)
                                        gCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACt                                                                                                     >  2:743994/1‑149 (MQ=255)
                                         ccAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACtt                                                                                                    >  1:1481809/1‑149 (MQ=255)
                                          cAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTa                                                                                                   >  2:1407597/1‑149 (MQ=255)
                                          cAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTa                                                                                                   >  2:1328663/1‑149 (MQ=255)
                                           aGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTAc                                                                                                  <  2:906498/149‑1 (MQ=255)
                                           aGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTAc                                                                                                  <  2:1481809/149‑1 (MQ=255)
                                             aCATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCg                                                                                                <  2:994966/149‑1 (MQ=255)
                                             aCATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCg                                                                                                >  1:994966/1‑149 (MQ=255)
                                                              ggCATACATAATCTCAGCCAGGCCCCGCAACAAGGGTGTGTGCTGTCTATTGGTAATTTCGACGGCGTGCATCGCGGTCTTCGCTCGCTGTTCCAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCtttt                                                                               <  2:697/149‑1 (MQ=255)
                                                                  tacataATCTCAGCCAGGTTTTGCCCGCCGGGTGTGTGCTGAATATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGaa                                                                           <  2:1651473/149‑1 (MQ=255)
                                                                   acataATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCGGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAAc                                                                          <  2:269422/149‑1 (MQ=255)
                                                                   acataATCTCAGCCAGGCCCCGCAAGAAGGGGGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAAc                                                                          <  2:1013148/149‑1 (MQ=255)
                                                                      taATCTCATCCAGGCCCCGAAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGCTGGTGATGCTTTTTGAACCTc                                                                       >  2:635562/1‑149 (MQ=255)
                                                                       aaTCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAAc                                                                          <  2:1273044/145‑1 (MQ=255)
                                                                       aaTCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAAc                                                                          >  1:1273044/1‑145 (MQ=255)
                                                                        aTCACAGCCATGCCCCGCAAGAAGGGTGTGTGCTGCCTATTGGGAATTTCGACGGCGTGCAGCGCGGTCAGCTCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCaa                                                                     <  2:563994/149‑1 (MQ=255)
                                                                         tctcAGCCTGGCCCCGCAAGTAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCCTCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAAc                                                                    <  1:445158/149‑1 (MQ=255)
                                                                         tctcAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAAc                                                                    >  2:603253/1‑149 (MQ=255)
                                                                         tctcAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAAc                                                                    <  1:531398/149‑1 (MQ=255)
                                                                         tctcAGCCAGGCCCCGCAAGAAGGGGGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAAc                                                                    <  1:891753/149‑1 (MQ=255)
                                                                           tcAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCACGCGCACCTTACCTGTGATGGTGATGCTTGTTGAACCTCAACCa                                                                  >  1:1232492/1‑149 (MQ=255)
                                                                            cAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCAc                                                                 >  2:119347/1‑149 (MQ=255)
                                                                             aGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACt                                                                <  1:1500526/149‑1 (MQ=255)
                                                                                     cccGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCCACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGt                                                        >  2:1340586/1‑149 (MQ=255)
                                                                                      ccGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGtt                                                       <  2:878087/149‑1 (MQ=255)
                                                                                      ccGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGtt                                                       <  2:380573/149‑1 (MQ=255)
                                                                                       cGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTTCAGGGCATGCAGGAAGAAGGGCGCAAGCGCAACATACCGGTGATGGTGATGCTTTTTGAACCTCGACCACTGGACCTGttt                                                      >  1:1234380/1‑149 (MQ=255)
                                                                                       cGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGttt                                                      <  2:1479195/149‑1 (MQ=255)
                                                                                       cGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGttt                                                      <  2:1158040/149‑1 (MQ=255)
                                                                                       cGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGttt                                                      <  1:1237320/149‑1 (MQ=255)
                                                                                       cGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGttt                                                      <  1:119347/149‑1 (MQ=255)
                                                                                        gCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTg                                                     <  2:1185658/149‑1 (MQ=255)
                                                                                            gtagGGTGTGTGCTGACTATTGGTAATTGCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGAAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTAc                                                 <  2:1105612/147‑1 (MQ=255)
                                                                                               gggTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGa                                              <  1:743994/149‑1 (MQ=255)
                                                                                                ggTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGAt                                             <  1:603253/149‑1 (MQ=255)
                                                                                                 gtgtgtGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATa                                            >  1:1213082/1‑149 (MQ=255)
                                                                                                        cTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGcccc                                     >  1:1333798/1‑149 (MQ=255)
                                                                                                        cTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGcccc                                     >  1:706586/1‑149 (MQ=255)
                                                                                                        cTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGcccc                                     >  1:677867/1‑149 (MQ=255)
                                                                                                          gACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCgg                                   >  2:102362/1‑149 (MQ=255)
                                                                                                          gACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCgg                                   >  2:407213/1‑149 (MQ=255)
                                                                                                           aCTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGc                                  <  2:334886/149‑1 (MQ=255)
                                                                                                             tATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCaa                                >  1:56525/1‑149 (MQ=255)
                                                                                                                 ggTAATTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACt                            >  1:65072/1‑149 (MQ=255)
                                                                                                                    aaTTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTAc                                                                                                  >  1:537061/1‑76 (MQ=255)
                                                                                                                    aaTTTCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTAc                                                                                                  <  2:537061/76‑1 (MQ=255)
                                                                                                                        tCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCgg                     >  1:737295/1‑149 (MQ=255)
                                                                                                                        tCGACGGCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGGGATGGTTATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCgg                     <  1:635562/149‑1 (MQ=255)
                                                                                                                           aCGGCGTGCATCGCGGTCATCGCTCGATGTTACAGGGCTTGCAGGAAGAAGGGCGCCAGCGCCACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGCTAAAGCCCCGGCAAGACTGACCCGGCTg                  >  1:1503421/1‑149 (MQ=255)
                                                                                                                             ggCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCg                <  2:1333798/149‑1 (MQ=255)
                                                                                                                              gCGTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCgg               <  1:1407597/149‑1 (MQ=255)
                                                                                                                                gTGCATCGCGGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCTGCTGCGGGa             >  2:1420620/1‑149 (MQ=255)
                                                                                                                                       gcgGTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGCAGGGCGCACGCGCAACTTACCGGTGATGGTGATGCTTTTTGCACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGCCCCGGCTGCGGGAAAAACTg      >  1:1266545/1‑149 (MQ=255)
                                                                                                                                         ggTCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCg    >  1:1207101/1‑149 (MQ=255)
                                                                                                                                           tCATCGCTCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGCCTGACCCGGCTGCGGGAAACACTGCGtt  >  1:1274549/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTT  >  NC_000913/21360‑21647

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