Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,331,979 T→C 8.3% T253A (ACG→GCG)  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,9790TC 8.3% 306.1 / 13.3 109T253A (ACG→GCG) yhbEEamA family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base T (53/46);  new base C (3/6);  total (56/53)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

Only 100 of 118 total aligned reads displayed.

CGATGACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGG  >  NC_000913/3331831‑3332127
                                                                                                                                                      |                                                                                                                                                    
cGATGACCAATGGCGGAATACATCGCGCCCGCGACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGtg                                                                                                                                                     >  2:362562/1‑150 (MQ=255)
  aTGACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGtgat                                                                                                                                                   <  2:819172/150‑1 (MQ=255)
  aTGACCAATGGCGGAATACATAGCGCACGCAACCACGACAAACGCACCGAGATAACCTAAAATGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGTGTGAGCGtgat                                                                                                                                                   <  1:94220/150‑1 (MQ=255)
   tGACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatg                                                                                                                                                  >  1:1483456/1‑150 (MQ=255)
   tGACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatg                                                                                                                                                  >  2:501294/1‑150 (MQ=255)
   tGACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatg                                                                                                                                                  <  1:1727692/150‑1 (MQ=255)
    gACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatga                                                                                                                                                 >  2:1299741/1‑150 (MQ=255)
      ccAATGTCGGAATTCATCGCGCCCGCTACCACGCCAAACGCACCTAGATATCCTATATGTTTTAACATCGGTCTGGCGAAGAATTCGGCCCAGGCCAGTGAAAAAAGAGCT‑G‑AATAAAACCGCGTAAACATTGGGGTGAGCTTGATTATc                                                                                                                                               <  2:469783/150‑1 (MQ=255)
      ccAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATc                                                                                                                                               <  1:1233310/150‑1 (MQ=255)
      ccAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATc                                                                                                                                               >  1:187878/1‑150 (MQ=255)
      ccAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGATTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGTGTGATGATc                                                                                                                                               <  1:1504593/150‑1 (MQ=255)
        aaTGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATcgc                                                                                                                                             >  1:835002/1‑150 (MQ=255)
        aaTGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCAAAAAGGCTTAACATCGGGCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGCGatcatcac                                                                                                                                             <  1:311262/150‑3 (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)
              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)
                 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)
                                  ccACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAg                                                                                                                   <  2:229857/150‑1 (MQ=255)
                                  ccACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAg                                                                                                                   <  2:156897/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)
                                            cGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTc                                                                                                         >  1:1482145/1‑150 (MQ=255)
                                                  gagaTAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAg                                                                                                   >  2:449611/1‑150 (MQ=255)
                                                          cTAAACGGATTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCat                                                                                           <  2:462808/150‑1 (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)
                                                                     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)
                                                                        aTCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTa                                                                             >  1:315734/1‑150 (MQ=255)
                                                                           ggTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTc                                                                          <  1:21471/150‑1 (MQ=255)
                                                                           ggTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTc                                                                          <  2:1483456/150‑1 (MQ=255)
                                                                               tGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTc                                                                      <  2:1576387/150‑1 (MQ=255)
                                                                               tGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTc                                                                      <  1:646741/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)
                                                                                    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)
                                                                                             ggccaggccaGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGCGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTaaa                                                        >  2:1099472/1‑150 (MQ=255)
                                                                                             ggccaggccaGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGGATTCAGTCCGCAAAAAATTaaa                                                        >  2:136414/1‑150 (MQ=255)
                                                                                               ccaggccagTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACa                                                      <  1:362562/150‑1 (MQ=255)
                                                                                               ccaggccagTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACa                                                      <  1:1057280/150‑1 (MQ=255)
                                                                                                  ggccagTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGc                                                   >  2:129455/1‑150 (MQ=255)
                                                                                                  ggccagTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGc                                                   <  2:1548057/150‑1 (MQ=255)
                                                                                                  ggccagTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGc                                                   >  1:362108/1‑150 (MQ=255)
                                                                                                  ggccagTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCGGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGc                                                   <  2:835002/150‑1 (MQ=255)
                                                                                                   gccagTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCg                                                  <  2:519007/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)
                                                                                                               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‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCCTCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAg                                 >  2:1208017/1‑150 (MQ=255)
                                                                                                                    t‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAg                                 >  1:247472/1‑150 (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)
                                                                                                                                     aaaCAGTGGGGTGAGCTGGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGATTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTcc                  <  2:1707292/150‑1 (MQ=255)
                                                                                                                                      aaCAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCa                 <  1:734577/150‑1 (MQ=255)
                                                                                                                                          cgcaggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCt             >  1:778045/5‑150 (MQ=255)
                                                                                                                                              ggTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGc         <  1:520524/150‑1 (MQ=255)
                                                                                                                                              ggTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGc         >  2:1155698/1‑150 (MQ=255)
                                                                                                                                               gTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTTTCACTCCAGGCTTTGcc        >  2:168681/1‑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)
                                                                                                                                                     gTGATGAGCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAggg  <  1:652747/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CGATGACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCT‑G‑AAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGG  >  NC_000913/3331831‑3332127

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: ‑