Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,322,049 T→A 8.7% intergenic (‑436/+23) 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,0490TA 8.7% 308.6 / 15.9 104intergenic (‑436/+23)yhbX/leuUputative hydrolase, inner membrane/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base T (51/44);  new base A (7/2);  total (58/46)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.91e-01

Only 100 of 107 total aligned reads displayed.

CAATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCTGGT  >  NC_000913/3321901‑3322193
                                                                                                                                                    |                                                                                                                                                
cAATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCtttt                                                                                                                                                 >  1:334514/1‑150 (MQ=255)
cAATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCtttt                                                                                                                                                 >  1:1341976/1‑150 (MQ=255)
cAATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCtttt                                                                                                                                                 >  2:341469/1‑150 (MQ=255)
 aaTTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTc                                                                                                                                                <  2:1214168/150‑1 (MQ=255)
 aaTTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTc                                                                                                                                                <  1:1538834/150‑1 (MQ=255)
     aTGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCg                                                                                                                                            <  1:1112594/150‑1 (MQ=255)
     aTGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCg                                                                                                                                            <  2:399313/150‑1 (MQ=255)
     aTGTCAGTAGAAGGGAAAAATTAATAACTAAGGCGTACCGCCACCATCGTTTCAAGGTACCTGCTACGAGTAATGCAACTGCACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCg                                                                                                                                            <  2:1618727/150‑1 (MQ=255)
      tGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGt                                                                                                                                           >  2:320181/1‑150 (MQ=255)
       gTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGtc                                                                                                                                          <  2:1391517/150‑1 (MQ=255)
        tCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGtct                                                                                                                                         >  1:1186183/1‑150 (MQ=255)
          aGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTctct                                                                                                                                       >  2:360058/1‑150 (MQ=255)
          aGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTctct                                                                                                                                       >  2:92256/1‑150 (MQ=255)
             aGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCtttt                                                                                                                                    >  2:89193/1‑150 (MQ=255)
             aGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCtttt                                                                                                                                    >  1:385876/1‑150 (MQ=255)
                aGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTg                                                                                                                                 <  2:614882/150‑1 (MQ=255)
                 gggAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTgg                                                                                                                                >  1:1013676/1‑150 (MQ=255)
                  ggAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGa                                                                                                                               >  2:27608/1‑150 (MQ=255)
                   gAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGaa                                                                                                                              >  1:816712/1‑150 (MQ=255)
                   gAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGaa                                                                                                                              >  1:671852/1‑150 (MQ=255)
                    aaaaaTTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAAt                                                                                                                             <  1:41666/150‑1 (MQ=255)
                     aaaaTTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAGAGACGCTGAAAAGCGTCTTTTTTcgttttgg                                                                                                                            <  1:429567/150‑8 (MQ=38)
                     aaaaTTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAGAGACGCTGAAAAGCGTCTTTTTTcgttttgg                                                                                                                            <  1:384744/150‑8 (MQ=38)
                        aTTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTgg                                                                                                                         <  2:390786/150‑1 (MQ=255)
                         tttATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGt                                                                                                                        <  2:1024732/150‑1 (MQ=255)
                         tttATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGt                                                                                                                        >  1:195031/1‑150 (MQ=255)
                           tataACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTAc                                                                                                                      <  1:92256/150‑1 (MQ=255)
                                    ggCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACgg                                                                                                             <  2:671852/150‑1 (MQ=255)
                                    ggCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACgg                                                                                                             >  1:1289447/1‑150 (MQ=255)
                                          ccGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTg                                                                                                       <  2:1441435/150‑1 (MQ=255)
                                            ggCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGac                                                                                                     >  2:379975/1‑149 (MQ=255)
                                            ggCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGaa                                                                                                     >  1:269318/1‑150 (MQ=255)
                                            ggCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGaa                                                                                                     >  1:1369068/1‑150 (MQ=255)
                                             gCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGATACCGAGGACGGGACTTGAAc                                                                                                    <  2:498885/150‑1 (MQ=255)
                                             gCACCATCGTTTCAAGGTACAAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAAc                                                                                                    >  2:526300/1‑150 (MQ=255)
                                                 cATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGt                                                                                                >  2:1229433/1‑150 (MQ=255)
                                                  aTCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTa                                                                                               >  1:504557/1‑150 (MQ=255)
                                                    cGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAg                                                                                             >  2:278665/1‑150 (MQ=255)
                                                    cGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAg                                                                                             <  2:515053/150‑1 (MQ=255)
                                                     gTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGc                                                                                            >  1:1608105/1‑150 (MQ=255)
                                                      tttCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGcc                                                                                           <  1:1229433/150‑1 (MQ=255)
                                                          aaGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTAt                                                                                       >  2:809237/1‑150 (MQ=255)
                                                           aGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTAtt                                                                                      <  1:1173974/150‑1 (MQ=255)
                                                              tACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTggg                                                                                   >  1:439732/1‑150 (MQ=255)
                                                               aCCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGc                                                                                  <  1:360058/150‑1 (MQ=255)
                                                               aCCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGc                                                                                  >  2:1004412/1‑150 (MQ=255)
                                                                ccAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCa                                                                                 <  2:1398609/150‑1 (MQ=255)
                                                                 cAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCAc                                                                                >  1:17667/1‑150 (MQ=255)
                                                                    cTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTac                                                                             >  1:1545445/1‑150 (MQ=255)
                                                                     tACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTacc                                                                            <  2:504557/150‑1 (MQ=255)
                                                                       cGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTaccac                                                                          <  1:187389/150‑1 (MQ=255)
                                                                           tAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCa                                                                      >  1:1527671/1‑150 (MQ=255)
                                                                           tAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCa                                                                      >  1:1197716/1‑150 (MQ=255)
                                                                            aaaGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAAcc                                                                                                   >  1:1388385/1‑120 (MQ=255)
                                                                            aaaGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAAcc                                                                                                   <  2:1388385/120‑1 (MQ=255)
                                                                            aaaGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCaa                                                                     >  1:105461/1‑150 (MQ=255)
                                                                                     ggACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCgtg                                                            <  1:27608/150‑1 (MQ=255)
                                                                                     ggACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCgtg                                                            <  2:1472793/150‑1 (MQ=255)
                                                                                        cGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGAAGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACCGGACTTGAACCCGTAACCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTca                                                         >  1:299710/1‑149 (MQ=255)
                                                                                         gagaTACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTa                                                        <  2:1186183/150‑1 (MQ=255)
                                                                                           gatacagatacCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTAcc                                                      <  2:816712/150‑1 (MQ=255)
                                                                                             tacagatacCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCaa                                                    <  2:1424589/150‑1 (MQ=255)
                                                                                                     ccTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTccacca                                            >  1:1370019/1‑150 (MQ=255)
                                                                                                      cTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCcaccac                                           >  2:1169/1‑150 (MQ=255)
                                                                                                      cTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCcaccac                                           >  2:1334589/1‑150 (MQ=255)
                                                                                                      cTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCcaccac                                           <  1:89193/150‑1 (MQ=255)
                                                                                                           aaCCATTCCTCAGACCAGGACCAAAACGAAAAAAGATGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCgg                                      >  1:865456/1‑150 (MQ=255)
                                                                                                           aaCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCgg                                      <  1:526300/150‑1 (MQ=255)
                                                                                                            aCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCTGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGc                                     <  1:1099663/150‑1 (MQ=255)
                                                                                                            aCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTc                                        <  2:814118/147‑1 (MQ=255)
                                                                                                            aCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTc                                        >  1:814118/1‑147 (MQ=255)
                                                                                                            aCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGc                                     >  2:474289/1‑150 (MQ=255)
                                                                                                               aTTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACg                                  <  2:1290535/150‑1 (MQ=255)
                                                                                                               aTTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACg                                  >  1:518013/1‑150 (MQ=255)
                                                                                                                  ccTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAt                               <  2:1289447/150‑1 (MQ=255)
                                                                                                                  ccTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAt                               >  1:898375/1‑150 (MQ=255)
                                                                                                                  ccTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAt                               >  1:518319/1‑150 (MQ=255)
                                                                                                                    tCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGTCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAtac                             <  2:334514/150‑1 (MQ=255)
                                                                                                                       gACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCCCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAtactac                          >  1:90090/1‑150 (MQ=255)
                                                                                                                              tccAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACttttttg                   <  1:477298/149‑3 (MQ=255)
                                                                                                                              aCCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAg                   <  1:847071/150‑1 (MQ=255)
                                                                                                                                cAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGtt                 >  1:1624406/1‑150 (MQ=255)
                                                                                                                                attccagAAAAGAGACGCTGAAAAGCGTCTTTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGtt                 >  1:1591641/7‑150 (MQ=38)
                                                                                                                                 aaaaCGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTc                >  2:130269/1‑150 (MQ=255)
                                                                                                                                   aaCGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCgg              <  2:1013676/150‑1 (MQ=255)
                                                                                                                                    cagAAAAGAGACGCTGAAAAGCGTCTTTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCggg             >  1:441674/3‑150 (MQ=255)
                                                                                                                                    cagAAAAGAGACGCTGAAAAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCggg             >  1:842240/3‑150 (MQ=255)
                                                                                                                                    aCGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCggg             <  2:105461/150‑1 (MQ=255)
                                                                                                                                     cGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGa            >  2:515036/1‑150 (MQ=255)
                                                                                                                                     cGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGa            >  2:772024/1‑150 (MQ=255)
                                                                                                                                     agAAAAGAGACGCTGAAAAGCGTCTTTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGa            >  1:1078507/2‑150 (MQ=255)
                                                                                                                                      taaaaaaGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGat           <  2:1534709/149‑1 (MQ=255)
                                                                                                                                      gAAAAGAGACGCTGAAAAGCGTCTTTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGat           >  1:36816/1‑150 (MQ=255)
                                                                                                                                       aaaagaGACGCTGAAAAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGata          >  1:147814/1‑150 (MQ=255)
                                                                                                                                       aaaagaGACGCTGAAAAGCGTCTCTATTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGata          >  2:906383/1‑150 (MQ=255)
                                                                                                                                         aaaaGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATc        <  2:209793/150‑1 (MQ=255)
                                                                                                                                           aaGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGc      <  2:1608105/150‑1 (MQ=255)
                                                                                                                                            aGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCt     <  1:348357/150‑1 (MQ=255)
                                                                                                                                              ccGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCTgg   <  1:824471/149‑1 (MQ=255)
                                                                                                                                               cGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCTGGt  >  1:1395504/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CAATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCTGGT  >  NC_000913/3321901‑3322193

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