Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 312,308 A→T 6.3% N66I (AAC→ATC)  ykgL → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913312,3080AT 6.3% 253.0 / 6.3 95N66I (AAC→ATC) ykgLuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (51/38);  new base T (3/3);  total (54/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

TAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACATATAGCCGTCCTTTTCGCT  >  NC_000913/312159‑312453
                                                                                                                                                     |                                                                                                                                                 
tAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATaaa                                                                                                                                                   >  2:953329/1‑150 (MQ=255)
tAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATaaa                                                                                                                                                   >  1:666873/1‑150 (MQ=255)
 aaaaGAGGGTTGATTCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGATTATTCTGAGAGAAACGTCATGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAATATTAATGTTGGTGTTATAAac                                                                                                                                                  <  1:213731/150‑1 (MQ=255)
     gagGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAacacac                                                                                                                                              <  1:579527/150‑1 (MQ=255)
       gggTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCa                                                                                                                                            <  2:946219/150‑1 (MQ=255)
       gggTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCa                                                                                                                                            >  2:385654/1‑150 (MQ=255)
        ggTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCaa                                                                                                                                           >  2:83538/1‑150 (MQ=255)
            gATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTGTTTATAACACCAACAtt                                                                                                                                       <  2:114111/150‑1 (MQ=255)
             aTGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTGTTTATAACACCAACATTa                                                                                                                                      <  1:135031/150‑1 (MQ=255)
               gCCAGCATACAGGGGCGAAGCTGGGCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTGTTTATAACACCAACATTAAt                                                                                                                                    <  2:780507/150‑1 (MQ=255)
                   gCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGa                                                                                                                                <  1:751674/150‑1 (MQ=255)
                     aTACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATc                                                                                                                              <  2:153977/150‑1 (MQ=255)
                      tACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCt                                                                                                                             >  2:1263649/1‑150 (MQ=255)
                        cAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATAAATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGAACTGa                                                                                                                           >  1:122971/1‑150 (MQ=255)
                          ggggCGAAGCTGGTCAACAAGTCAACATTAAGATAATGCAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAAACCACCAACATTAATCCGATCTGatt                                                                                                                         >  2:365457/1‑150 (MQ=255)
                           gggCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGatta                                                                                                                        <  2:864865/150‑1 (MQ=255)
                           gggCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGatta                                                                                                                        <  2:614934/150‑1 (MQ=255)
                            ggCGAAGCTGGTCAACAAGGCAACATTAAGATAATGGAATATTCGGAGAGACACGTCAGGCAGTTGGCCTCAAATTAACAGGAAGAATATAACCCCATGAAAATTAATGTTGGTGGTATAAACACACCAACATTAATCCGAGCTGattat                                                                                                                       <  1:184770/150‑1 (MQ=255)
                              cGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGAttatta                                                                                                                     <  1:605096/150‑1 (MQ=255)
                                 aGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGTAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCACATGAAAATTAATGTTTGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATc                                                                                                                  <  2:685021/150‑1 (MQ=255)
                                 aGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAACATTAATGTTGGTGTTATTAACCCACCATCATGAGTCCGATGTGATGATTAATc                                                                                                                  >  2:1146953/1‑150 (MQ=255)
                                 aGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTTATTATTAAta                                                                                                                  >  2:637085/1‑149 (MQ=255)
                                 aGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATc                                                                                                                  <  1:133471/150‑1 (MQ=255)
                                  gCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCa                                                                                                                 >  2:403914/1‑150 (MQ=255)
                                     ggTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGAc                                                                                                              >  1:641055/1‑150 (MQ=255)
                                     ggTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGGTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATACGATCTGATTATTAATCataa                                                                                                              >  1:287877/1‑147 (MQ=255)
                                     ggTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGACAGGCAGTTGGCCTCAAATGAACAGGAAGAAAATAACCCCAGGAAAATTAATGTTGGTGTTAAAAACACACCAACATTAATACGATATGATTATTAATCAGaa                                                                                                              >  1:980098/1‑149 (MQ=255)
                                      gTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAAAGGAAGAATATATGCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAAACCGATCTGATTATTAATAAGACa                                                                                                             >  2:592498/1‑150 (MQ=255)
                                        caacaaGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGTCATTTTGCCTCAAATGCTCAGGAAGAATATATCACCAGGAAAATTAATGTTTTTGTTATAAACACACCAACATTAATCCGATCTGATTAATAATCAGatcgg                                                                                                           <  2:196326/150‑5 (MQ=255)
                                        caacaaGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATAAGAaaaa                                                                                                           >  1:345767/1‑150 (MQ=255)
                                              gTCAACATTAAGATAATGGAATATTCGGCGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCTCCAGGAAAATTAATGTTGGTGTTATAAAAACACCAACATTAATACGATCTGATTATTAATCAGACAAATATCGt                                                                                                     >  1:368472/1‑150 (MQ=255)
                                              gTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGt                                                                                                     >  1:287099/1‑150 (MQ=255)
                                                cAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTa                                                                                                   >  1:878178/1‑150 (MQ=255)
                                                 aaCATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTaa                                                                                                  <  2:723873/150‑1 (MQ=255)
                                                   cATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACGCCAACATTAATACGATCTGATTATTAATCAGACAAATATCGTTAATc                                                                                                >  2:1209979/1‑150 (MQ=255)
                                                     ttAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCg                                                                                              >  1:430288/1‑150 (MQ=255)
                                                     ttAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCg                                                                                              <  2:499457/150‑1 (MQ=255)
                                                        aGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCg                                                                                           >  1:468510/1‑150 (MQ=255)
                                                           tAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACCCCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAAt                                                                                        >  2:502666/1‑150 (MQ=255)
                                                              tGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAAtttt                                                                                     >  2:351134/1‑150 (MQ=255)
                                                               ggAATATTCGGAGAGAACCGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATGATTAATCAGACAAATATCGTTAATCCGGCGAATTTTg                                                                                    >  1:1293136/1‑150 (MQ=255)
                                                               ggAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATTTATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTg                                                                                    >  2:435162/1‑150 (MQ=255)
                                                                  atatTCGGAGAGAACCGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAAt                                                                                 <  2:584857/150‑1 (MQ=255)
                                                                   tatTCGGAGAGAAAGGTCAGGCAGTTGGCCTGAAATGAACAGGAAGAATATATTCCCAGGAAATTTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGTCAAATATCGTTAATCCGTCGAATTTTGAATc                                                                                <  2:596834/150‑1 (MQ=255)
                                                                   tatTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATc                                                                                >  2:1281783/1‑150 (MQ=255)
                                                                       cGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACg                                                                            <  2:266967/150‑1 (MQ=255)
                                                                       cGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACg                                                                            >  1:119504/1‑150 (MQ=255)
                                                                        ggAGAGAAAAGTCACGCAGTTGGCCTAAAATGATCAAGAAGAATTTATAACAAGGAAAATTAATGTTGGTGTTTTAAAAACAAAAAAAATAATCAGATATGATTTATAATCATACAAAAAAATGTAAACCGTCATAGTTTGAATCAACGt                                                                           >  2:1187972/1‑150 (MQ=255)
                                                                            agaAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCACAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGttttt                                                                       <  1:83538/150‑1 (MQ=255)
                                                                             gaAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGtttttt                                                                      >  1:529345/1‑150 (MQ=255)
                                                                             gaAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGtttttt                                                                      >  2:423716/1‑150 (MQ=255)
                                                                               aaCGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCt                                                                    >  1:215540/1‑150 (MQ=255)
                                                                               aaCGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCt                                                                    <  1:1263649/150‑1 (MQ=255)
                                                                                 cGTCAGGCAGGTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTa                                                                  <  1:1092096/150‑1 (MQ=255)
                                                                                  gTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGATAATCCGTCGAATTTTGAATAAAATTTTTTTCTTAc                                                                 >  1:90879/1‑150 (MQ=255)
                                                                                       gCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACTAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGAccc                                                            <  2:125674/150‑1 (MQ=255)
                                                                                            tGGCCTCAAATGAAAATTAAGAATATTTCCCAAGGTAAATTAATTGTGGTGATATAAACACAGCAACATTAATCCGATCTGATTATTAATCAGACTAATATCTTTAATCCGGCTAATTTTGAATCAACGTTTTTTCTTACGACCCTGAAc                                                       <  1:190479/150‑1 (MQ=255)
                                                                                              gCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACgg                                                     <  2:311024/150‑1 (MQ=255)
                                                                                               ccTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGc                                                    >  2:396958/1‑150 (MQ=255)
                                                                                               ccTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATAAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGc                                                    >  2:1064228/1‑150 (MQ=255)
                                                                                                cTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGcc                                                   >  2:746511/1‑150 (MQ=255)
                                                                                                 tCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCTTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCt                                                  <  2:510691/150‑1 (MQ=255)
                                                                                                 tCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCt                                                  >  1:1292162/1‑150 (MQ=255)
                                                                                                 tCAAATGAACAGGAACAATATATCCACAGGAAAATTAAATTTGGTGTTTTAAACACACCAACATTAATCCGATCTGATTAATAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCt                                                  <  2:121050/150‑1 (MQ=255)
                                                                                                  cAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCtt                                                 >  2:1161702/1‑150 (MQ=255)
                                                                                                  cAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATACGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCtt                                                 <  2:1136635/150‑1 (MQ=255)
                                                                                                  aaaaTGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCtt                                                 <  1:704919/149‑1 (MQ=255)
                                                                                                   aaaTGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTa                                                <  2:716642/150‑1 (MQ=255)
                                                                                                    attGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACTACCCTGAACTGCCTTaa                                               >  2:567629/3‑150 (MQ=255)
                                                                                                      tGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTaaaa                                             >  2:303598/1‑150 (MQ=255)
                                                                                                        aaCAGGAAGAATATATCCACAGGAAAATTAATGTTGGTGTTATAAACACACAAAGATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTGTTTTCTTACGACCCTGAACGGCCTTAAAACg                                           <  1:646342/150‑1 (MQ=255)
                                                                                                          cAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCATTAAAACGTg                                         >  2:116429/1‑150 (MQ=255)
                                                                                                           aGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTgg                                        >  1:1056646/1‑150 (MQ=255)
                                                                                                            ggAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTTATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGa                                       >  2:446883/1‑150 (MQ=255)
                                                                                                            ggAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGa                                       <  1:403914/150‑1 (MQ=255)
                                                                                                             gaagaaTATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGAAAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTAAGACACTGAACGGCATTAAAACGTGGAt                                      >  1:160951/1‑150 (MQ=255)
                                                                                                              aagaaTATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGAtt                                     <  1:1161702/150‑1 (MQ=255)
                                                                                                                           cAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATc                        >  1:1058741/1‑150 (MQ=255)
                                                                                                                             ggAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATcac                      <  1:1260661/150‑1 (MQ=255)
                                                                                                                              gAAAATTAATGTTGGTGTTATAAACACACCATCATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATcaca                     <  1:351134/150‑1 (MQ=255)
                                                                                                                              gAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATcaca                     <  2:929705/150‑1 (MQ=255)
                                                                                                                               aaaaTTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACat                    >  1:50428/1‑150 (MQ=255)
                                                                                                                                aaaTTAATGTTGGTGTTATAAACACACCAACCTTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGCCCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACata                   >  1:940817/1‑150 (MQ=255)
                                                                                                                                 gaTTAATGTTGGTGTGTTTATAACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACataa                  >  2:386214/2‑149 (MQ=39)
                                                                                                                                 aaTTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACatat                  <  1:746511/150‑1 (MQ=255)
                                                                                                                                    tAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAACACGGGGATTAGATTTACAAATCACATATAGc               <  1:562625/150‑1 (MQ=255)
                                                                                                                                    tAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACATATAGc               <  2:902847/150‑1 (MQ=255)
                                                                                                                                     aaTGTTGGTGTTATAAACACACCAACATTAATCCGATGTGATTATTAATCAGACAAATATCGTTAATCCGTCGTATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAATACGTGGATTAGATATACAAATCACATATAGcc              >  2:643720/1‑150 (MQ=255)
                                                                                                                                     aaTGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACATATAGcc              >  2:643718/1‑150 (MQ=255)
                                                                                                                                      aTGTTGGTGTGTTTATAACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACATATAGCCg             >  1:1024215/1‑150 (MQ=255)
                                                                                                                                       tgttgttgtTATAAACACAACAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACATATAGCCGt            <  1:231280/150‑1 (MQ=255)
                                                                                                                                       tGTTGGTGTGTTTATAACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACATATAGCCGt            >  1:1125898/1‑150 (MQ=255)
                                                                                                                                             tgtTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTAAATCCGTCGAATATTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACATATAGCCGTcggttt      >  1:727917/1‑145 (MQ=255)
                                                                                                                                              gtTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACATATAGCCGTCCTTTTc     >  2:1241616/1‑150 (MQ=255)
                                                                                                                                              gtTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACATATAGCCGTCCTTTTc     >  1:544476/1‑150 (MQ=255)
                                                                                                                                              gtTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACATATAGCCGTCCTTTTc     >  2:1250453/1‑150 (MQ=255)
                                                                                                                                               ttATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACATATAGCCGTCCTTTTCg    <  1:1027702/150‑1 (MQ=255)
                                                                                                                                                tataAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACATATAGCCGTCCTTTTCGc   >  1:451866/1‑150 (MQ=255)
                                                                                                                                                 ataAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACATATAGCCGTCCTTTTCGCt  <  2:1000307/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAATCCGATCTGATTATTAATCAGACAAATATCGTTAATCCGTCGAATTTTGAATCAACGTTTTTTCTTACGACCCTGAACGGCCTTAAAACGTGGATTAGATTTACAAATCACATATAGCCGTCCTTTTCGCT  >  NC_000913/312159‑312453

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