Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,331,987 G→A 6.8% A250V (GCG→GTG)  yhbE ← EamA family inner membrane putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,331,9870GA 6.8% 319.8 / 6.7 104A250V (GCG→GTG) yhbEEamA family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base G (51/44);  new base A (3/4);  total (55/49)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 107 total aligned reads displayed.

AATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGGAAGAGC  >  NC_000913/3331839‑3332133
                                                                                                                                                      |                                                                                                                                                  
aaTGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATcgc                                                                                                                                                   >  1:835002/1‑150 (MQ=255)
 aTGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATcgcg                                                                                                                                                  >  1:1637841/1‑150 (MQ=255)
  tGGCGGAATACAACGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCgcgc                                                                                                                                                 <  2:1394968/150‑1 (MQ=255)
   ggCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCt                                                                                                                                                >  1:1310521/1‑150 (MQ=255)
    gCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTc                                                                                                                                               <  2:860118/150‑1 (MQ=255)
      ggAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAATCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTctc                                                                                                                                             >  2:578106/1‑148 (MQ=255)
      ggAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCT‑AAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGCGATCATCACGCTCAcc                                                                                                                                            <  1:458606/150‑1 (MQ=255)
       gAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCAcc                                                                                                                                            >  1:1009312/1‑150 (MQ=255)
       gAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGAGGTGAGCGCGATCATCACGCTCAcc                                                                                                                                            <  1:102971/150‑1 (MQ=255)
        aaTACATCGCTCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGCGATCATCACGCTCAccc                                                                                                                                           <  2:52463/150‑2 (MQ=255)
        aaTACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGCGATCATCACGCTCAccc                                                                                                                                           <  1:773500/150‑2 (MQ=255)
         aTACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTg                                                                                                                                          <  1:627350/150‑1 (MQ=255)
            cATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTgcgc                                                                                                                                       <  2:345281/150‑1 (MQ=255)
             aTCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGct                                                                                                                                      >  1:1461801/1‑150 (MQ=255)
             aTCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGct                                                                                                                                      >  1:519007/1‑150 (MQ=255)
                gcgcCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcc                                                                                                                                   <  2:889436/150‑1 (MQ=255)
                 cgcCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcct                                                                                                                                  >  1:812051/1‑150 (MQ=255)
                  gcCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcctg                                                                                                                                 >  2:942878/1‑150 (MQ=255)
                  gcCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcctg                                                                                                                                 <  1:1608287/150‑1 (MQ=255)
                   cccGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑GAAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcctg                                                                                                                                 <  1:801205/150‑1 (MQ=255)
                   cccGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcctgc                                                                                                                                >  1:1337959/1‑150 (MQ=255)
                         aCCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGa                                                                                                                          >  2:302385/1‑150 (MQ=255)
                          ccACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCTGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAg                                                                                                                         <  2:860038/150‑1 (MQ=255)
                          ccACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAg                                                                                                                         <  2:156897/150‑1 (MQ=255)
                          ccACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAg                                                                                                                         <  2:229857/150‑1 (MQ=255)
                           cACGACAAACGCACCTAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATATAAGATCT‑G‑AAAAAAACAGCGTAAACCGTGGGGTGAGCTTGATGATCCCGCTCCCCTCCGCTTCCTGCCACCGAGc                                                                                                                        >  2:1055131/1‑150 (MQ=255)
                            aCGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGcc                                                                                                                       <  1:964603/150‑1 (MQ=255)
                              gACAAACGCACCGAGATAACCTAAAAGGTTTAACAGCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCAt                                                                                                                     <  2:48438/150‑1 (MQ=255)
                                          gagaTAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAg                                                                                                         >  2:449611/1‑150 (MQ=255)
                                          gagaTAACCT‑AAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAgg                                                                                                        >  1:568950/1‑150 (MQ=255)
                                                  cTAAACGGATTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCat                                                                                                 <  2:462808/150‑1 (MQ=255)
                                                  cTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGCGCCATCGCTTCCGCCAGGGCGCCat                                                                                                 >  1:431498/1‑150 (MQ=255)
                                                     aaaGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATc                                                                                              <  1:1422967/150‑1 (MQ=255)
                                                     aaaGGTTTAACATAGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTCA‑AAAAAAACAGAGTAAACAGTGGGGTGAGCGTGATGATCGAGCTCACCTGAGCTGACTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCataa                                                                                               >  1:369689/1‑149 (MQ=255)
                                                       agttttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATGT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCt                                                                                            <  2:886861/147‑1 (MQ=255)
                                                        ggTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTa                                                                                           >  2:1599519/1‑150 (MQ=255)
                                                        ggTTTAACATCGGTCTCGCGAAGAAATCGGGCCAGGCCAGTGATAAAATATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTa                                                                                           <  1:1349977/150‑1 (MQ=255)
                                                          tttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                                                                         >  1:1590138/1‑150 (MQ=255)
                                                          tttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                                                                         >  1:441806/1‑150 (MQ=255)
                                                           ttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAAGAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCa                                                                                        >  2:295521/1‑150 (MQ=255)
                                                            tAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCaa                                                                                       >  1:873059/1‑150 (MQ=255)
                                                             aaCCTCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAg                                                                                      <  2:1505586/150‑1 (MQ=255)
                                                              aCATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAgg                                                                                     <  1:1605383/150‑1 (MQ=255)
                                                               cATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTCCCAAGGt                                                                                    >  2:734577/1‑150 (MQ=255)
                                                                   ggTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTc                                                                                <  2:1483456/150‑1 (MQ=255)
                                                                   ggTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTc                                                                                <  1:21471/150‑1 (MQ=255)
                                                                       tGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTc                                                                            <  1:646741/150‑1 (MQ=255)
                                                                       tGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTc                                                                            <  2:1576387/150‑1 (MQ=255)
                                                                        ggCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTcc                                                                           >  2:704688/1‑150 (MQ=255)
                                                                          cGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGc                                                                         >  2:1019457/1‑150 (MQ=255)
                                                                          cGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGc                                                                         >  1:638770/1‑150 (MQ=255)
                                                                           gaagaaATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCa                                                                        >  2:1304385/1‑150 (MQ=255)
                                                                            aagaaATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaa                                                                       >  1:1548057/1‑150 (MQ=255)
                                                                              gaaATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaa                                                                     >  2:781172/1‑150 (MQ=255)
                                                                                aaTCGGGCCAGGCCAGTGCTAAAAGATCT‑G‑AAAAAAACCGCGGAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGCGCCCTCGCTTCCGCCAGGGCGCCATATCCTACCACGGTATTCAGTCCGCaaaaaa                                                                   >  1:1577217/1‑150 (MQ=255)
                                                                                    gggCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTaa                                                               >  1:1207627/1‑150 (MQ=255)
                                                                                     ggccaggccaGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGGATTCAGTCCGCAAAAAATTaaa                                                              >  2:136414/1‑150 (MQ=255)
                                                                                       ccaggccagTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACa                                                            <  1:1057280/150‑1 (MQ=255)
                                                                                       ccaggccagTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACa                                                            <  1:362562/150‑1 (MQ=255)
                                                                                          ggccagTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGc                                                         >  1:362108/1‑150 (MQ=255)
                                                                                          ggccagTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGc                                                         >  2:129455/1‑150 (MQ=255)
                                                                                          ggccagTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGc                                                         <  2:1548057/150‑1 (MQ=255)
                                                                                          ggccagTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCGGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGc                                                         <  2:835002/150‑1 (MQ=255)
                                                                                           gccagTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCg                                                        <  2:411240/150‑1 (MQ=255)
                                                                                             cagTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAg                                                      <  1:1299741/150‑1 (MQ=255)
                                                                                             cagTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAg                                                      <  1:658794/150‑1 (MQ=255)
                                                                                             cagTGATAAAAGAGCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAg                                                      <  2:294928/150‑1 (MQ=255)
                                                                                                  aTAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGcc                                                 >  1:1363472/1‑150 (MQ=255)
                                                                                                    aaaaGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGCGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCTAGCTGCAAg                                               >  2:57914/1‑150 (MQ=255)
                                                                                                     aaaGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGt                                              <  2:315734/150‑1 (MQ=255)
                                                                                                       aGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCTCGCTCACCTGCGCTGCCTGCCAGCGAGACATCGCTTCCGCCAGGGCGCCATATCCTACCAATGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTgg                                            <  2:1169700/150‑1 (MQ=255)
                                                                                                        gATCT‑G‑AAACACACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGc                                           <  1:1012091/150‑1 (MQ=255)
                                                                                                        gATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGc                                           <  2:885485/150‑1 (MQ=255)
                                                                                                          tCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCAGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTa                                         >  1:431688/1‑150 (MQ=255)
                                                                                                            t‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCTCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAg                                       <  1:1025979/150‑1 (MQ=255)
                                                                                                            t‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGCGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTAGTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGTTaac                                       >  1:268301/1‑149 (MQ=255)
                                                                                                            t‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCCTCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAg                                       >  2:1208017/1‑150 (MQ=255)
                                                                                                            t‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAg                                       >  1:247472/1‑150 (MQ=255)
                                                                                                            cccccacaaacCAGCGTAAACACCCGGCTGCGCGCGAGGCGCGCGCGCACCGGCGCGGCCTGCCCGCGAGCCATCGCGTCCGCCAGGGCGCCATATCCTACCCAGGTATTCCGCCCGCAAACAATTAAACATGCGAGCTGCCAGTGGCTa                                         <  2:473278/143‑1 (MQ=255)
                                                                                                              g‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGc                                      <  1:302385/150‑1 (MQ=255)
                                                                                                                 aaaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTg                                    >  1:159575/1‑150 (MQ=255)
                                                                                                                 aaaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCACACCTGCCCTGACTCCCAGCGAGCCATCGCTTCCGCCAGAGCGCCATATCCTACCAAGGAATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTg                                    >  1:1152335/1‑150 (MQ=255)
                                                                                                                   aaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGTTATTCAGTCCGCAAAAAATTAAACATGCTATCTGCCAGTGGCTAAGCTgcg                                  >  2:1047915/1‑150 (MQ=255)
                                                                                                                        aGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCACATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTacc                             <  2:187878/150‑3 (MQ=255)
                                                                                                                          cGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCAc                           >  1:811376/1‑150 (MQ=255)
                                                                                                                          cGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCAc                           >  1:1511995/1‑150 (MQ=255)
                                                                                                                             aaaCAGTGGGGTGAGCTGGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGATTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTcc                        <  2:1707292/150‑1 (MQ=255)
                                                                                                                              aaCAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCa                       <  1:734577/150‑1 (MQ=255)
                                                                                                                                  cgcaggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCt                   >  1:778045/5‑150 (MQ=255)
                                                                                                                                    ggggTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCAAGTGGCTAAGCTGCGCTATCACTCCACGAttt                 >  2:1728452/1‑150 (MQ=255)
                                                                                                                                      ggTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGc               <  1:520524/150‑1 (MQ=255)
                                                                                                                                      ggTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGc               >  2:1155698/1‑150 (MQ=255)
                                                                                                                                      cgTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTAGTCAGTACGCAAAAAATTAAACATGCCAGCTGCCAGAGGCTAAGCTCCGCTATCACTCCAGGCTGTGc               >  1:1410143/2‑150 (MQ=255)
                                                                                                                                        tGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCa             >  2:1008136/1‑150 (MQ=255)
                                                                                                                                        tGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCa             >  2:564437/1‑150 (MQ=255)
                                                                                                                                        ggAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCa             <  1:258927/149‑1 (MQ=255)
                                                                                                                                             gTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAggg        >  2:1421384/1‑150 (MQ=255)
                                                                                                                                                 tgatCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTACACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGGAagc    >  1:98964/1‑149 (MQ=255)
                                                                                                                                                  gatCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCGAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGGAtgag   >  1:1548647/1‑150 (MQ=255)
                                                                                                                                                   atCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGGAAGAgc  <  1:449611/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
AATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGGAAGAGC  >  NC_000913/3331839‑3332133

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