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.

CGATGACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGG  >  NC_000913/3331831‑3332127
                                                                                                                                                      |                                                                                                                                                    
cGATGACCAATGGCGGAATACATCGCGCCCGCGACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtg                                                                                                                                                     >  2:362562/1‑150 (MQ=255)
  aTGACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgat                                                                                                                                                   <  2:819172/150‑1 (MQ=255)
  aTGACCAATGGCGGAATACATAGCGCACGCAACCACGACAAACGCACCGAGATAACCTAAAATGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGTGTGAGCGtgat                                                                                                                                                   <  1:94220/150‑1 (MQ=255)
   tGACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatg                                                                                                                                                  >  1:1483456/1‑150 (MQ=255)
   tGACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatg                                                                                                                                                  >  2:501294/1‑150 (MQ=255)
   tGACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatg                                                                                                                                                  <  1:1727692/150‑1 (MQ=255)
    gACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatga                                                                                                                                                 >  2:1299741/1‑150 (MQ=255)
      ccAATGTCGGAATTCATCGCGCCCGCTACCACGCCAAACGCACCTAGATATCCTATATGTTTTAACATCGGTCTGGCGAAGAATTCGGCCCAGGCCAGTGAAAAAAGAGCTGAATAAAACCGCGTAAACATTGGGGTGAGCTTGATTATc                                                                                                                                               <  2:469783/150‑1 (MQ=255)
      ccAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATc                                                                                                                                               <  1:1233310/150‑1 (MQ=255)
      ccAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATc                                                                                                                                               >  1:187878/1‑150 (MQ=255)
      ccAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGATTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGTGTGATGATc                                                                                                                                               <  1:1504593/150‑1 (MQ=255)
        aaTGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATcgc                                                                                                                                             >  1:835002/1‑150 (MQ=255)
        aaTGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCAAAAAGGCTTAACATCGGGCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGCGatcatcac                                                                                                                                             <  1:311262/150‑3 (MQ=255)
          tGGCGGAATACAACGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCgcgc                                                                                                                                           <  2:1394968/150‑1 (MQ=255)
           ggCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCt                                                                                                                                          >  1:1310521/1‑150 (MQ=255)
            gCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTc                                                                                                                                         <  2:860118/150‑1 (MQ=255)
              ggAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGCGATCATCACGCTCAcc                                                                                                                                      <  1:458606/150‑1 (MQ=255)
               gAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCAcc                                                                                                                                      >  1:1009312/1‑150 (MQ=255)
               gAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGAGGTGAGCGCGATCATCACGCTCAcc                                                                                                                                      <  1:102971/150‑1 (MQ=255)
                aaTACATCGCTCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGCGATCATCACGCTCAccc                                                                                                                                     <  2:52463/150‑2 (MQ=255)
                 aTACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTg                                                                                                                                    <  1:627350/150‑1 (MQ=255)
                    cATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTgcgc                                                                                                                                 <  2:345281/150‑1 (MQ=255)
                     aTCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGct                                                                                                                                >  1:1461801/1‑150 (MQ=255)
                     aTCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGct                                                                                                                                >  1:519007/1‑150 (MQ=255)
                        gcgcCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcc                                                                                                                             <  2:889436/150‑1 (MQ=255)
                         cgcCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcct                                                                                                                            >  1:812051/1‑150 (MQ=255)
                          gcCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcctg                                                                                                                           >  2:942878/1‑150 (MQ=255)
                          gcCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcctg                                                                                                                           <  1:1608287/150‑1 (MQ=255)
                           cccGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcctg                                                                                                                           <  1:801205/150‑1 (MQ=255)
                           cccGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcctgc                                                                                                                          >  1:1337959/1‑150 (MQ=255)
                                 aCCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGa                                                                                                                    >  2:302385/1‑150 (MQ=255)
                                  ccACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAg                                                                                                                   <  2:229857/150‑1 (MQ=255)
                                  ccACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAg                                                                                                                   <  2:156897/150‑1 (MQ=255)
                                   cACGACAAACGCACCTAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATATAAGATCTGAAAAAAACAGCGTAAACCGTGGGGTGAGCTTGATGATCCCGCTCCCCTCCGCTTCCTGCCACCGAGc                                                                                                                  >  2:1055131/1‑150 (MQ=255)
                                    aCGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGcc                                                                                                                 <  1:964603/150‑1 (MQ=255)
                                      gACAAACGCACCGAGATAACCTAAAAGGTTTAACAGCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCAt                                                                                                               <  2:48438/150‑1 (MQ=255)
                                            cGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTc                                                                                                         >  1:1482145/1‑150 (MQ=255)
                                                  gagaTAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAg                                                                                                   >  2:449611/1‑150 (MQ=255)
                                                          cTAAACGGATTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCat                                                                                           <  2:462808/150‑1 (MQ=255)
                                                             aaaGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATc                                                                                        <  1:1422967/150‑1 (MQ=255)
                                                             aaaGGTTTAACATAGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTCAAAAAAAACAGAGTAAACAGTGGGGTGAGCGTGATGATCGAGCTCACCTGAGCTGACTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCataa                                                                                         >  1:369689/1‑149 (MQ=255)
                                                               agttttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATGTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCt                                                                                      <  2:886861/147‑1 (MQ=255)
                                                                ggTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTa                                                                                     >  2:1599519/1‑150 (MQ=255)
                                                                ggTTTAACATCGGTCTCGCGAAGAAATCGGGCCAGGCCAGTGATAAAATATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTa                                                                                     <  1:1349977/150‑1 (MQ=255)
                                                                  tttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                                                                   >  1:1590138/1‑150 (MQ=255)
                                                                  tttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                                                                   >  1:441806/1‑150 (MQ=255)
                                                                   ttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAAGAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCa                                                                                  >  2:295521/1‑150 (MQ=255)
                                                                     aaCCTCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAg                                                                                <  2:1505586/150‑1 (MQ=255)
                                                                      aCATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAgg                                                                               <  1:1605383/150‑1 (MQ=255)
                                                                       cATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTCCCAAGGt                                                                              >  2:734577/1‑150 (MQ=255)
                                                                        aTCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTa                                                                             >  1:315734/1‑150 (MQ=255)
                                                                           ggTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTc                                                                          <  1:21471/150‑1 (MQ=255)
                                                                           ggTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTc                                                                          <  2:1483456/150‑1 (MQ=255)
                                                                               tGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTc                                                                      <  2:1576387/150‑1 (MQ=255)
                                                                               tGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTc                                                                      <  1:646741/150‑1 (MQ=255)
                                                                                ggCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTcc                                                                     >  2:704688/1‑150 (MQ=255)
                                                                                  cGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGc                                                                   >  2:1019457/1‑150 (MQ=255)
                                                                                  cGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGc                                                                   >  1:638770/1‑150 (MQ=255)
                                                                                    aagaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaa                                                                 >  1:1548057/1‑150 (MQ=255)
                                                                                      gaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaa                                                               >  2:781172/1‑150 (MQ=255)
                                                                                        aaTCGGGCCAGGCCAGTGCTAAAAGATCTGAAAAAAACCGCGGAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGCGCCCTCGCTTCCGCCAGGGCGCCATATCCTACCACGGTATTCAGTCCGCaaaaaa                                                             >  1:1577217/1‑150 (MQ=255)
                                                                                             ggccaggccaGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGCGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTaaa                                                        >  2:1099472/1‑150 (MQ=255)
                                                                                             ggccaggccaGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGGATTCAGTCCGCAAAAAATTaaa                                                        >  2:136414/1‑150 (MQ=255)
                                                                                               ccaggccagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACa                                                      <  1:362562/150‑1 (MQ=255)
                                                                                               ccaggccagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACa                                                      <  1:1057280/150‑1 (MQ=255)
                                                                                                  ggccagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGc                                                   >  2:129455/1‑150 (MQ=255)
                                                                                                  ggccagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGc                                                   <  2:1548057/150‑1 (MQ=255)
                                                                                                  ggccagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGc                                                   >  1:362108/1‑150 (MQ=255)
                                                                                                  ggccagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCGGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGc                                                   <  2:835002/150‑1 (MQ=255)
                                                                                                   gccagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCg                                                  <  2:519007/150‑1 (MQ=255)
                                                                                                   gccagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCg                                                  <  2:411240/150‑1 (MQ=255)
                                                                                                     cagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAg                                                <  1:1299741/150‑1 (MQ=255)
                                                                                                     cagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAg                                                <  1:658794/150‑1 (MQ=255)
                                                                                                     cagTGATAAAAGAGCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAg                                                <  2:294928/150‑1 (MQ=255)
                                                                                                          aTAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGcc                                           >  1:1363472/1‑150 (MQ=255)
                                                                                                            aaaaGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGCGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCTAGCTGCAAg                                         >  2:57914/1‑150 (MQ=255)
                                                                                                               aGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCTCGCTCACCTGCGCTGCCTGCCAGCGAGACATCGCTTCCGCCAGGGCGCCATATCCTACCAATGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTgg                                      <  2:1169700/150‑1 (MQ=255)
                                                                                                                gATCTGAAACACACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGc                                     <  1:1012091/150‑1 (MQ=255)
                                                                                                                gATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGc                                     <  2:885485/150‑1 (MQ=255)
                                                                                                                  tCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCAGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTa                                   >  1:431688/1‑150 (MQ=255)
                                                                                                                    tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCTCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAg                                 <  1:1025979/150‑1 (MQ=255)
                                                                                                                    tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCCTCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAg                                 >  2:1208017/1‑150 (MQ=255)
                                                                                                                    tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAg                                 >  1:247472/1‑150 (MQ=255)
                                                                                                                      gAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGc                                <  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)
                                                                                                                                                      |                                                                                                                                                    
CGATGACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGG  >  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: