Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,331,983 G→C 6.5% I251M (ATC→ATG 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,9830GC 6.5% 339.1 / 12.0 108I251M (ATC→ATGyhbEEamA family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base G (53/47);  new base C (3/4);  total (56/52)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

Only 100 of 116 total aligned reads displayed.

TGACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGGAAGAG  >  NC_000913/3331834‑3332132
                                                                                                                                                       |                                                                                                                                                     
tGACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatg                                                                                                                                                       >  1:1483456/1‑150 (MQ=255)
tGACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatg                                                                                                                                                       >  2:501294/1‑150 (MQ=255)
tGACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatg                                                                                                                                                       >  2:1512054/1‑150 (MQ=255)
 gACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatga                                                                                                                                                      >  2:1299741/1‑150 (MQ=255)
   ccAATGTCGGAATTCATCGCGCCCGCTACCACGCCAAACGCACCTAGATATCCTATATGTTTTAACATCGGTCTGGCGAAGAATTCGGCCCAGGCCAGTGAAAAAAGAGCTGAATAAAACCGCGTAAACATTGGGGTGAGCTTGATTATc                                                                                                                                                    <  2:469783/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)
      aTGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATcgcg                                                                                                                                                 >  1:1637841/1‑150 (MQ=255)
       tGGCGGAATACAACGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCgcgc                                                                                                                                                <  2:1394968/150‑1 (MQ=255)
        ggCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCt                                                                                                                                               >  1:1310521/1‑150 (MQ=255)
         gCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTc                                                                                                                                              <  2:860118/150‑1 (MQ=255)
           ggAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAATCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTctc                                                                                                                                            >  2:578106/1‑148 (MQ=255)
            gAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCAcc                                                                                                                                           >  1:1009312/1‑150 (MQ=255)
            gAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGAGGTGAGCGCGATCATCACGCTCAcc                                                                                                                                           <  1:102971/150‑1 (MQ=255)
             aaTACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGCGATCATCACGCTCAccc                                                                                                                                          <  1:773500/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)
                       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)
                               ccACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCTGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAg                                                                                                                        <  2:860038/150‑1 (MQ=255)
                               ccACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAg                                                                                                                        <  2:229857/150‑1 (MQ=255)
                               ccACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAg                                                                                                                        <  2:156897/150‑1 (MQ=255)
                                 aCGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGcc                                                                                                                      <  1:964603/150‑1 (MQ=255)
                                   gACAAACGCACCGAGATAACCTAAAAGGTTTAACAGCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCAt                                                                                                                    <  2:48438/150‑1 (MQ=255)
                                         cGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTc                                                                                                              >  1:1482145/1‑150 (MQ=255)
                                               gagaTAACCTAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAgg                                                                                                       >  1:568950/1‑150 (MQ=255)
                                                       cTAAACGGATTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCat                                                                                                <  2:462808/150‑1 (MQ=255)
                                                       cTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGCGCCATCGCTTCCGCCAGGGCGCCat                                                                                                >  1:431498/1‑150 (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:441806/1‑150 (MQ=255)
                                                               tttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                                                                        >  1:1590138/1‑150 (MQ=255)
                                                                ttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAAGAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCa                                                                                       >  2:295521/1‑150 (MQ=255)
                                                                 tAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCaa                                                                                      >  1:873059/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)
                                                                        ggTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTc                                                                               <  2:1483456/150‑1 (MQ=255)
                                                                        ggTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTc                                                                               <  1:21471/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                                                                        >  1:638770/1‑150 (MQ=255)
                                                                               cGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGc                                                                        >  2:1019457/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)
                                                                                         gggCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTaa                                                              >  1:1207627/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                                                        >  1:362108/1‑150 (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)
                                                                                                          aaaGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGt                                             <  2:315734/150‑1 (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)
                                                                                                                 tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGCGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTAGTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGTTaac                                      >  1:268301/1‑149 (MQ=255)
                                                                                                                 tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCCTCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAg                                      >  2:1208017/1‑150 (MQ=255)
                                                                                                                 tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAg                                      >  1:247472/1‑150 (MQ=255)
                                                                                                                 cccccacaaacCAGCGTAAACACCCGGCTGCGCGCGAGGCGCGCGCGCACCGGCGCGGCCTGCCCGCGAGCCATCGCGTCCGCCAGGGCGCCATATCCTACCCAGGTATTCCGCCCGCAAACAATTAAACATGCGAGCTGCCAGTGGCTa                                        <  2:473278/143‑1 (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)
                                                                                                                               cGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCAc                          >  1:1511995/1‑150 (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)
                                                                                                                                         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)
                                                                                                                                            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)
                                                                                                                                                      tgatCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTACACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGGAag  >  1:98964/1‑149 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
TGACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGGAAGAG  >  NC_000913/3331834‑3332132

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: