Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,322,051 C→A 6.0% intergenic (‑438/+21) yhbX ← / ← leuU putative hydrolase, inner membrane/tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,322,0510CA 6.0% 392.0 / 8.0 115intergenic (‑438/+21)yhbX/leuUputative hydrolase, inner membrane/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base C (48/60);  new base A (4/3);  total (52/63)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 116 total aligned reads displayed.

AATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCTGGTCGGCTT  >  NC_000913/3321902‑3322199
                                                                                                                                                     |                                                                                                                                                    
aaTTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTc                                                                                                                                                      <  2:584242/150‑1 (MQ=255)
    aTGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCg                                                                                                                                                  >  2:824185/1‑150 (MQ=255)
      gTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGtc                                                                                                                                                <  1:294366/150‑1 (MQ=255)
       tCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTAAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGtct                                                                                                                                               >  2:281020/1‑150 (MQ=255)
          gTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCtt                                                                                                                                            <  2:1333063/150‑1 (MQ=255)
              aaGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCt                                                                                                                                        <  2:45089/150‑1 (MQ=255)
               aGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTg                                                                                                                                       <  1:1453303/150‑1 (MQ=255)
                  gAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGaa                                                                                                                                    >  1:1064068/1‑150 (MQ=255)
                  gAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGaa                                                                                                                                    >  1:1089186/1‑150 (MQ=255)
                  gAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGaa                                                                                                                                    >  1:798862/1‑150 (MQ=255)
                   aaaaaTTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAGAGACGCTGAAAAGCGTCTTTTTTcgttttg                                                                                                                                   <  1:1117692/150‑7 (MQ=38)
                       aTTTATAACTAAGGCGTACCGGCACCATATTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTgg                                                                                                                               <  1:507782/150‑1 (MQ=255)
                        tttATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGt                                                                                                                              >  1:1343951/1‑150 (MQ=255)
                        tttATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGt                                                                                                                              <  2:1499117/150‑1 (MQ=255)
                          tataACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTAc                                                                                                                            >  2:1291878/1‑150 (MQ=255)
                              aCTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAg                                                                                                                        <  2:1372482/150‑1 (MQ=255)
                              aCTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAg                                                                                                                        <  1:83117/150‑1 (MQ=255)
                               cTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAgg                                                                                                                       <  2:54452/150‑1 (MQ=255)
                               cTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAgg                                                                                                                       <  2:1446369/150‑1 (MQ=255)
                               cTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAgg                                                                                                                       >  1:710659/1‑150 (MQ=255)
                                 aaGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGAc                                                                                                                     >  1:493854/1‑150 (MQ=255)
                                   ggCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACgg                                                                                                                   <  2:556430/150‑1 (MQ=255)
                                    gCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACggg                                                                                                                  <  2:798862/150‑1 (MQ=255)
                                        aCCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACtt                                                                                                              >  1:93486/1‑150 (MQ=255)
                                         ccGGAACCATCTTTTCAAGTTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTg                                                                                                             <  2:653344/150‑1 (MQ=255)
                                          cGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGa                                                                                                            <  2:224704/150‑1 (MQ=255)
                                            gCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAAc                                                                                                          >  2:933764/1‑150 (MQ=255)
                                             cACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAAcc                                                                                                         <  1:1459915/150‑1 (MQ=255)
                                                    gTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGc                                                                                                  <  2:112713/150‑1 (MQ=255)
                                                    gTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGc                                                                                                  <  2:1089186/150‑1 (MQ=255)
                                                     tttCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGcc                                                                                                 <  1:1010202/150‑1 (MQ=255)
                                                      ttCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGccc                                                                                                <  2:458150/150‑1 (MQ=255)
                                                         aaGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTAt                                                                                             <  2:365342/150‑1 (MQ=255)
                                                           ggTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTg                                                                                           >  1:586598/1‑150 (MQ=255)
                                                           ggTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTg                                                                                           <  1:1143999/150‑1 (MQ=255)
                                                             tACCAGCTACGAGTAAAGCAACGGGACGAGATCCCGATACCTGACACCCATTCCTCCGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGGCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTggg                                                                                         <  2:1303569/150‑1 (MQ=255)
                                                              aCCAGCTACGAGTAAAGCAACTTGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCAATTGGGc                                                                                        <  2:1206977/150‑1 (MQ=255)
                                                              aCCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGc                                                                                        <  2:372105/150‑1 (MQ=255)
                                                              aCCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGc                                                                                        <  2:767700/150‑1 (MQ=255)
                                                              aCCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGc                                                                                        <  2:1081311/150‑1 (MQ=255)
                                                               ccAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGTGACTTGAACCCGTAAGCCCTATTGGGCa                                                                                       <  1:62472/150‑1 (MQ=255)
                                                               ccAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCa                                                                                       >  2:718879/1‑150 (MQ=255)
                                                                 aGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACt                                                                                     <  2:419086/150‑1 (MQ=255)
                                                                  gCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTa                                                                                    <  1:708485/150‑1 (MQ=255)
                                                                   cTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTac                                                                                   >  2:157905/1‑150 (MQ=255)
                                                                    tACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTCTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTacc                                                                                  <  2:1239053/150‑1 (MQ=255)
                                                                     aCGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTacca                                                                                 >  1:1228521/1‑150 (MQ=255)
                                                                      cGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTaccac                                                                                >  2:76901/1‑150 (MQ=255)
                                                                       gAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGGACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTaccacc                                                                               >  1:1186505/1‑150 (MQ=255)
                                                                        aGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCt                                                                              >  1:330399/1‑150 (MQ=255)
                                                                          tAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCa                                                                            <  2:345151/150‑1 (MQ=255)
                                                                          tAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCa                                                                            <  2:74873/150‑1 (MQ=255)
                                                                           aaaGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCaa                                                                           >  2:1172102/1‑150 (MQ=255)
                                                                             aGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAgg                                                                         >  2:1024613/1‑150 (MQ=255)
                                                                             aGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAgg                                                                         <  2:1290377/150‑1 (MQ=255)
                                                                              gCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGt                                                                        <  1:953285/150‑1 (MQ=255)
                                                                                    ggACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCgtg                                                                  <  2:1364743/150‑1 (MQ=255)
                                                                                       cGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCt                                                               <  1:132101/150‑1 (MQ=255)
                                                                                              cagatacCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAAtt                                                        >  2:738327/1‑150 (MQ=255)
                                                                                              cagatacCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAAtt                                                        <  2:52983/150‑1 (MQ=255)
                                                                                               agatacCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTGTTTTCTGGAATTTGGTACCGAGGACGGGAGTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTc                                                       <  1:1406710/150‑1 (MQ=255)
                                                                                               agatacCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTc                                                       <  1:933764/150‑1 (MQ=255)
                                                                                                  tacCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTccac                                                    <  2:510268/150‑1 (MQ=255)
                                                                                                   acCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTccacc                                                   <  1:76901/150‑1 (MQ=255)
                                                                                                    ccTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTccacca                                                  >  2:157925/1‑150 (MQ=255)
                                                                                                      tGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCaccacc                                                <  2:722436/150‑1 (MQ=255)
                                                                                                          aaCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCgg                                            <  1:1286237/150‑1 (MQ=255)
                                                                                                           aCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGc                                           >  1:1288765/1‑150 (MQ=255)
                                                                                                           aCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGc                                           >  1:671135/1‑150 (MQ=255)
                                                                                                            ccATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCCCCTCAAGGTAGCGTGTCTACCAATTCCACCACCACGGCa                                          >  1:1405568/1‑150 (MQ=255)
                                                                                                             cATTCCTCAGACAAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCGCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCAc                                         >  2:1111529/1‑150 (MQ=255)
                                                                                                             aaTTCCTCAGACCAGTTCCAAAACTAAAAAAGACGGTTTTCAGAGTCTCTTTTCTGGAATTTGGTACCGAGGAGGGGAGTTGATCCCGTAAGGCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCAGCTCGGCag                                         <  2:508389/149‑2 (MQ=255)
                                                                                                               ttCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACgg                                       <  2:713629/150‑1 (MQ=255)
                                                                                                                 ccTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAt                                     <  2:12873/150‑1 (MQ=255)
                                                                                                                 ccTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAt                                     >  1:348427/1‑150 (MQ=255)
                                                                                                                  cTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAta                                    <  2:1421533/150‑1 (MQ=255)
                                                                                                                  cTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAta                                    >  2:252089/1‑150 (MQ=255)
                                                                                                                    cAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAtact                                  >  2:155102/1‑150 (MQ=255)
                                                                                                                      gACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGACCCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAtactac                                >  1:160411/1‑150 (MQ=255)
                                                                                                                      gACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCCCCTCAAGGTAGCGTGTCTACCAAGTCCACCCCCTCGGCACGGAtactac                                >  1:154975/1‑150 (MQ=255)
                                                                                                                      gACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAtactac                                <  1:153535/150‑1 (MQ=255)
                                                                                                                        ccAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACtt                              >  1:63533/1‑150 (MQ=255)
                                                                                                                              ccAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGt                        >  1:1324361/1‑150 (MQ=255)
                                                                                                                              aattccagAAAAGAGACGCTGAAAAGCGTCTTTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCCCCACCTCGGCACGGATACTACTTTTTAGt                        >  1:304993/8‑150 (MQ=38)
                                                                                                                                ttccagAAAAGAGACGCTGAAAAGCGTCTTTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTc                      >  1:1694/6‑150 (MQ=255)
                                                                                                                                aaaaCGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTc                      >  1:683382/1‑150 (MQ=255)
                                                                                                                                aaaaCGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTc                      >  1:64991/1‑150 (MQ=255)
                                                                                                                                aaaaCGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAAACTACTTTTTAGTTc                      >  1:600757/1‑150 (MQ=255)
                                                                                                                                  ccagAAAAGAGACGCTGAAAAGCGTCTTTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCgg                    >  1:94015/4‑150 (MQ=255)
                                                                                                                                  ccagAAAAGAGACGCTGAAAAGCGTCTTTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCgg                    >  1:113853/4‑150 (MQ=255)
                                                                                                                                  aaCGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCgg                    >  1:829397/1‑150 (MQ=255)
                                                                                                                                   cagAAAAGAGACGCTGAAAAGCGTCTTTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCggg                   <  2:1694/148‑1 (MQ=255)
                                                                                                                                     gAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGCACTTGAACCCGTAAGCCCTAATGGGCACTACCACCTCAAGGTAGCCTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTAAGTTCGGGat                 >  1:993923/1‑150 (MQ=255)
                                                                                                                                        aaaaGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATc              >  1:764478/1‑150 (MQ=255)
                                                                                                                                        aaaaGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATc              >  2:92819/1‑150 (MQ=255)
                                                                                                                                          aaGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGc            >  1:559725/1‑150 (MQ=255)
                                                                                                                                                cTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACCGATACTACATTTTAGTTCGGGATATCGCTGGTCg      >  2:1157338/1‑150 (MQ=255)
                                                                                                                                                   ttCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCTGGTCGGCt   >  1:586884/1‑150 (MQ=255)
                                                                                                                                                    tCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCTGGTCGGCtt  <  1:934890/150‑1 (MQ=255)
                                                                                                                                                    tCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCTGGTCGGCtt  <  1:824185/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCTGGTCGGCTT  >  NC_000913/3321902‑3322199

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: