Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,322,051 C→A 9.3% 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 9.3% 325.0 / 15.1 102intergenic (‑438/+21)yhbX/leuUputative hydrolase, inner membrane/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base C (50/42);  new base A (8/2);  total (58/44)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

AATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCTGGTCGGCTTA  >  NC_000913/3321902‑3322200
                                                                                                                                                     |                                                                                                                                                     
aaTTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTc                                                                                                                                                       <  1:1538834/150‑1 (MQ=255)
    aTGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCg                                                                                                                                                   <  2:399313/150‑1 (MQ=255)
    aTGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCg                                                                                                                                                   <  1:1112594/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:92256/1‑150 (MQ=255)
         aGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTctct                                                                                                                                              >  2:360058/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:671852/1‑150 (MQ=255)
                  gAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGaa                                                                                                                                     >  1:816712/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)
                      aaTTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTg                                                                                                                                 >  1:1419192/1‑150 (MQ=255)
                       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)
                                     cGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACCAGATACAGATACCTGACAACAATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTAAGCGTCTCTTTTCTGTAAATTGGTACCGAGGACGGGa                                                                                                                  >  2:477298/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:515053/150‑1 (MQ=255)
                                                   cGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAg                                                                                                    >  2:278665/1‑150 (MQ=255)
                                                     tttCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGcc                                                                                                  <  1:1229433/150‑1 (MQ=255)
                                                       tCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCt                                                                                                >  1:624769/1‑150 (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)
                                                               ccAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCa                                                                                        >  2:1141045/1‑150 (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                                                             <  1:341469/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:1334589/1‑150 (MQ=255)
                                                                                                     cTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCcaccac                                                  <  1:89193/150‑1 (MQ=255)
                                                                                                     cTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCcaccac                                                  >  2:1169/1‑150 (MQ=255)
                                                                                                          aaCCATTCCTCAGACCAGGACCAAAACGAAAAAAGATGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCgg                                             >  1:865456/1‑150 (MQ=255)
                                                                                                          aaCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCgg                                             <  1:526300/150‑1 (MQ=255)
                                                                                                           aCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCTGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGc                                            <  1:1099663/150‑1 (MQ=255)
                                                                                                           aCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTc                                               >  1:814118/1‑147 (MQ=255)
                                                                                                           aCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTc                                               <  2:814118/147‑1 (MQ=255)
                                                                                                           aCCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGc                                            >  2:474289/1‑150 (MQ=255)
                                                                                                              aTTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACg                                         >  1:518013/1‑150 (MQ=255)
                                                                                                              aTTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACg                                         <  2:1290535/150‑1 (MQ=255)
                                                                                                                 ccTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAt                                      <  2:1289447/150‑1 (MQ=255)
                                                                                                                 ccTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAt                                      >  1:898375/1‑150 (MQ=255)
                                                                                                                      gACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCCCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGAtactac                                 >  1:90090/1‑150 (MQ=255)
                                                                                                                        ccAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACtt                               >  2:1475171/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)
                                                                                                                                                     cAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCTGGTCGGCTTa  >  2:1314025/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATTATGTCAGTAGAAGGGAAAAATTTATAACTAAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTAAGCCCTATTGGGCACTACCACCTCAAGGTAGCGTGTCTACCAATTCCACCACCTCGGCACGGATACTACTTTTTAGTTCGGGATATCGCTGGTCGGCTTA  >  NC_000913/3321902‑3322200

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: