Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,597,121 C→T 6.4% intergenic (+109/+29) lgoD → / ← opgB L‑galactonate oxidoreductase/phosphoglycerol transferase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,597,1210CT 6.4% 430.4 / 12.2 140intergenic (+109/+29)lgoD/opgBL‑galactonate oxidoreductase/phosphoglycerol transferase I
Reads supporting (aligned to +/- strand):  ref base C (55/75);  new base T (6/3);  total (61/79)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-01

Only 100 of 144 total aligned reads displayed.

TAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCA  >  NC_000913/4596973‑4597265
                                                                                                                                                    |                                                                                                                                                
tAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCAcc                                                                                                                                                 >  1:761460/1‑150 (MQ=255)
 aaCAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCAccc                                                                                                                                                >  1:1969299/1‑150 (MQ=255)
 aaCAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCAccc                                                                                                                                                >  1:1811738/1‑150 (MQ=255)
  aCAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCt                                                                                                                                               >  1:519852/1‑150 (MQ=255)
  aCAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCt                                                                                                                                               >  1:1961689/1‑150 (MQ=255)
      tCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGc                                                                                                                                           >  1:872399/1‑150 (MQ=255)
       cGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCg                                                                                                                                          <  2:1216791/150‑1 (MQ=255)
       cGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGGTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCg                                                                                                                                          <  1:2478388/150‑1 (MQ=255)
         tGAGTTAATTACAGGCGGAATTAGTTTCTGATAAGGCGCGATTATTAAAGTTGCCATGCAGCGGCCGGGGAAGGGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTTGAAGGCGTGTAGGTCGCACCCTATGCGAc                                                                                                                                        <  2:2251709/150‑1 (MQ=255)
         ggAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGTTCGCACCCTATGCGAc                                                                                                                                        <  1:406904/149‑1 (MQ=255)
          gAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGAcc                                                                                                                                       <  2:1275071/150‑1 (MQ=255)
            gTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTTCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACACTATGCGACCTa                                                                                                                                     >  1:1154781/1‑150 (MQ=255)
            gTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTa                                                                                                                                     >  2:431683/1‑150 (MQ=255)
             ttaattaaAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTAGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTac                                                                                                                                    <  1:1710345/150‑1 (MQ=255)
             ttaattaaAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGGTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTac                                                                                                                                    <  1:1701757/150‑1 (MQ=255)
               aattaaAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTacac                                                                                                                                  >  1:27717/1‑150 (MQ=255)
                attaaAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCATAGGGTGCGACCTacacg                                                                                                                                 <  1:1828089/150‑2 (MQ=255)
                 ttaaAGGCGTAGTTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCATAGGGTGCGACCTacacgc                                                                                                                                <  2:35810/150‑3 (MQ=255)
                  taaAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACAtc                                                                                                                               <  2:774697/150‑1 (MQ=255)
                    aaGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACAtcag                                                                                                                             >  2:1322359/1‑150 (MQ=255)
                     aGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACAtcagt                                                                                                                            >  1:2300999/1‑150 (MQ=255)
                      ggCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACAtcagtc                                                                                                                           >  1:638863/1‑150 (MQ=255)
                      ggCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACAtcagtca                                                                                                                          >  2:1755117/1‑150 (MQ=255)
                         gTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAgg                                                                                                                        <  2:2168025/150‑1 (MQ=255)
                               aCTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTAc                                                                                                                  >  2:279261/1‑150 (MQ=255)
                               aCTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTAc                                                                                                                  <  1:992630/150‑1 (MQ=255)
                               aCTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTAc                                                                                                                  >  1:1849450/1‑150 (MQ=255)
                                  ttCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCt                                                                                                               <  1:299485/150‑1 (MQ=255)
                                       aTAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACg                                                                                                          >  1:367593/1‑150 (MQ=255)
                                        tAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGt                                                                                                         >  1:417612/1‑150 (MQ=255)
                                           ggCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCt                                                                                                      <  1:1126938/150‑1 (MQ=255)
                                             cGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTAc                                                                                                    >  1:652638/1‑150 (MQ=255)
                                              gagaTTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTAcc                                                                                                   <  2:183740/150‑1 (MQ=255)
                                                  ttattaAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACtt                                                                                               <  2:1171281/150‑1 (MQ=255)
                                                   tattaAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACttt                                                                                              >  1:346121/1‑150 (MQ=255)
                                                   tattaAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACttt                                                                                              <  2:977507/150‑1 (MQ=255)
                                                     ttaAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTaa                                                                                            <  2:1961689/150‑1 (MQ=255)
                                                          gTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTcc                                                                                       >  2:970368/1‑150 (MQ=255)
                                                          gTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTcc                                                                                       >  2:929742/1‑150 (MQ=255)
                                                          gTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTcc                                                                                       <  2:1423506/150‑1 (MQ=255)
                                                             gCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTccacc                                                                                    <  1:1546105/150‑1 (MQ=255)
                                                             gCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTccacc                                                                                    <  1:2041567/150‑1 (MQ=255)
                                                               cATGCAGCGTCCGGCGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCAt                                                                                  <  1:2474748/150‑1 (MQ=255)
                                                                  gCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGcc                                                                               <  1:327596/150‑1 (MQ=255)
                                                                  gCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGcc                                                                               <  2:2195401/150‑1 (MQ=255)
                                                                   cAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGCGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCg                                                                              <  2:2029971/150‑1 (MQ=255)
                                                                    aGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGa                                                                             <  1:1863345/150‑1 (MQ=255)
                                                                       gTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTATAAGTCTTGTATGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATcc                                                                          <  2:213294/150‑1 (MQ=255)
                                                                            gggAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGc                                                                     <  1:1615372/150‑1 (MQ=255)
                                                                            gggAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGc                                                                     <  2:161255/150‑1 (MQ=255)
                                                                                 ttgtTGGGCGCTTTTTTTTTTGTTTCTTAATAATGTGTTGTATGGCGTAGAAGGCTTTTAGGTCTTACCCTTTTTGACCTAGACATCAGTCAGGATTTACGCTTGACTTTCTAACACTTTAATTTCCACCATTCCGATCCGGAGCTTAcg                                                                <  2:441928/149‑1 (MQ=255)
                                                                                  tgtTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTAcgc                                                               <  1:1757142/150‑1 (MQ=255)
                                                                                   gtTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTAcgcg                                                              >  2:7312/1‑150 (MQ=255)
                                                                                   gtTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTAcgcg                                                              <  1:2081252/150‑1 (MQ=255)
                                                                                     tGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACgcggc                                                            <  2:1816640/150‑1 (MQ=255)
                                                                                      gggCGCTGTTTTTTTTGTTTCTTAATAATGTTTTGTAAGCCGTAGTAGGCGTGTAGGTCGCACCCTATGTGACCTACACATCAGTCAGTAATTACCCTTCACGTTCTACCTCTTTAATTTCCACCATTCCGATCCCTAGCTTACgctgcg                                                           <  2:2266829/150‑1 (MQ=255)
                                                                                       ggCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGa                                                          <  1:1573037/150‑1 (MQ=255)
                                                                                       ggCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGa                                                          >  1:651868/1‑150 (MQ=255)
                                                                                        gcgcTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAg                                                         <  1:733054/150‑1 (MQ=255)
                                                                                         cgcTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAgt                                                        >  1:877091/1‑150 (MQ=255)
                                                                                             tgtttttttGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTcc                                                    <  1:1755117/148‑1 (MQ=255)
                                                                                                 tttttGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAgg                                                >  2:1197401/1‑150 (MQ=255)
                                                                                                 tttttGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAgg                                                >  2:2391427/1‑150 (MQ=255)
                                                                                                 tttttGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAgg                                                <  2:27717/150‑1 (MQ=255)
                                                                                                  ttttGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGa                                               <  1:1666452/150‑1 (MQ=255)
                                                                                                  ttttGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGa                                               <  1:2409545/150‑1 (MQ=255)
                                                                                                  ttttGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGa                                               <  2:1716271/150‑1 (MQ=255)
                                                                                                      gTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTAACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATAtt                                           >  1:1814002/1‑150 (MQ=255)
                                                                                                        ttCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGGCCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTcc                                         >  1:1287249/1‑150 (MQ=255)
                                                                                                        ttCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTcc                                         >  1:1944270/1‑150 (MQ=255)
                                                                                                         tCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTccc                                        >  2:870060/1‑150 (MQ=255)
                                                                                                             aataatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCg                                    <  1:931090/150‑1 (MQ=255)
                                                                                                             aataatGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCg                                    <  2:2241645/150‑1 (MQ=255)
                                                                                                                 atGTGTTGTAAGCCGTAGAATGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGATTGTCAGAGGATATTCCCTTCgttgg                                <  2:1969299/150‑1 (MQ=255)
                                                                                                                 atGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttgg                                <  2:367593/150‑1 (MQ=255)
                                                                                                                   gtgtTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttggtt                              >  2:678387/1‑150 (MQ=255)
                                                                                                                    tgttgtAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttggttg                             <  2:519852/150‑1 (MQ=255)
                                                                                                                    tgttgtAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCgttggttg                             <  2:328238/150‑1 (MQ=255)
                                                                                                                      ttgtAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTgag                           <  2:652638/150‑1 (MQ=255)
                                                                                                                       tgtAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTgaga                          <  2:2075616/150‑1 (MQ=255)
                                                                                                                         tAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACa                        <  2:31234/150‑1 (MQ=255)
                                                                                                                          aaGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAg                       <  2:289104/150‑1 (MQ=255)
                                                                                                                          aaGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAg                       >  2:1888352/1‑150 (MQ=255)
                                                                                                                             ccGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGtt                    >  1:2204325/1‑150 (MQ=255)
                                                                                                                             ccGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGtt                    <  2:1894815/150‑1 (MQ=255)
                                                                                                                                  taaGGCGTGTAGGTCGCACGCTATGCGACCTACACATCAGTCAGTAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCGGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCcggc               <  2:2215184/149‑1 (MQ=255)
                                                                                                                                  gAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCcggc               >  1:1023672/1‑150 (MQ=255)
                                                                                                                                  gAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCcggc               <  1:1964912/150‑1 (MQ=255)
                                                                                                                                  gAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCcggc               >  2:1830809/1‑150 (MQ=255)
                                                                                                                                   aaGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCCGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGAGCCCGAGCTTCCGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCcggcg              <  2:2415076/150‑1 (MQ=255)
                                                                                                                                        gtgtAGGTGGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGcaca         <  2:1969897/150‑1 (MQ=255)
                                                                                                                                        gtgtAGGTCGCATAGGGTGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGcaca         >  1:2418488/1‑150 (MQ=255)
                                                                                                                                         tgtAGGTCGCATAGGGTGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACaa        >  2:494505/1‑150 (MQ=255)
                                                                                                                                         tgtAGGTCGCATAGGGTGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACaa        >  2:1650749/1‑150 (MQ=255)
                                                                                                                                           tAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAAtt      >  1:1664744/1‑150 (MQ=255)
                                                                                                                                           tAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAAtt      >  1:2343748/1‑150 (MQ=255)
                                                                                                                                            aGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTTATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTa     <  1:1322359/150‑1 (MQ=255)
                                                                                                                                             ggTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTAc    >  2:769402/1‑150 (MQ=255)
                                                                                                                                              gTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTAcc   <  1:472098/150‑1 (MQ=255)
                                                                                                                                               tCGCACCCTATGCGGCCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCa  <  2:1287249/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCGCTGTTTTTTTTGTTTCTTAATAATGTGTTGTAAGCCGTAGAAGGCGTGTAGGTCGCACCCTATGCGACCTACACATCAGTCAGGAATTACCCTTCACGTTCTACCACTTTAATTTCCACCATGCCGATCCCGAGCTTACGCGGCGAGTGTCCGAGGATATTCCCTTCGTTGGTTGAGACAGGTTCCGGCGGCACAATTACCA  >  NC_000913/4596973‑4597265

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