Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 805,959 A→G 5.2% intergenic (+195/+39) ybhJ → / ← ybhC aconitase family protein/acyl‑CoA thioesterase, lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913805,9590AG 5.2% 341.7 / 7.1 119intergenic (+195/+39)ybhJ/ybhCaconitase family protein/acyl‑CoA thioesterase, lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (59/51);  new base G (2/4);  total (64/55)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01

Only 100 of 120 total aligned reads displayed.

GATGCTTGCGCATCTTATCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGTTATCATC  >  NC_000913/805813‑806108
                                                                                                                                                  |                                                                                                                                                     
gatgcTTGCGCATCTTATCCGACCTACGAATCGCATCGATTCTGTAGGCCAGATAAGGCATTTTCGCAACATCCGGCACTTATTGTCGGATGCGATGCTAGCGCATCTTATACGACATACAAATCGCATAGAACCGTAGGCCGGATAAgg                                                                                                                                                    >  2:76415/1‑150 (MQ=38)
    cTTGCGCATCTTATCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAAACGCATCGAACCGTAGGCCGGATAAGGAGtt                                                                                                                                                >  1:329155/1‑150 (MQ=255)
     ttGAGCATCTTATCCGACCTACGAATCGCCTCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGttt                                                                                                                                               <  2:224097/150‑1 (MQ=255)
       gcgcATCTTATCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTAc                                                                                                                                             >  1:1015372/1‑150 (MQ=255)
           aTCTTATCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                         >  1:426690/1‑150 (MQ=255)
            tCTTATCCGCCCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCCAACCGTAGGCCGGATAAGGCGTTTAcgccgc                                                                                                                                        >  2:1798481/1‑150 (MQ=255)
              ttATCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCAt                                                                                                                                      >  1:1536024/1‑150 (MQ=255)
               tATCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATGCGGCGTAAAcgccttatc                                                                                                                                     <  1:927339/150‑6 (MQ=37)
                aTCCGACCTACGAATCGCATCGAATCTGTAGGCGAGATAAGGCATTTTCGCAGCCTCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAACTAGCATCGAACCGTAGGCCGGATTAGGCGTTTTCGCCGCATcc                                                                                                                                    >  2:1186638/1‑150 (MQ=25)
                aTCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                    <  2:1313936/150‑1 (MQ=37)
                  ccGACCTACGAATCGCATCTAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                  <  1:617273/150‑1 (MQ=25)
                   cGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGCTGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCGTCGAACCGTCGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                 >  2:346527/1‑150 (MQ=255)
                   cGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                 >  2:1393398/1‑150 (MQ=255)
                    gACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                <  1:1215558/150‑1 (MQ=255)
                    gACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCATATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                <  2:519421/150‑2 (MQ=37)
                          cGAATCGCATCGAATCTGTATGCCAGATAAGGCATTTTCGCATCATGCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGt                                                                                                                          <  2:1386144/150‑1 (MQ=255)
                           taaTCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGtt                                                                                                                         <  2:773102/149‑1 (MQ=255)
                           gAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGtt                                                                                                                         >  1:1639343/1‑150 (MQ=255)
                           gAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGtt                                                                                                                         <  2:1070523/150‑1 (MQ=255)
                             aTCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCCTCCGGCAAATAGTTaa                                                                                                                       >  1:1194106/1‑150 (MQ=255)
                             aTCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTaa                                                                                                                       >  2:394858/1‑150 (MQ=255)
                               cGCATCGAATCTGTATGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGAATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAAtt                                                                                                                     <  2:159907/150‑1 (MQ=255)
                               cGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAAtt                                                                                                                     >  2:636002/1‑150 (MQ=255)
                                gCATCGTATGTTTAGGGCAGATTATGCATTTTCGCAGCATCCGTCACTTATTGTCGGATGCCATGCTTGCTCATCTTATCCGACCTACAAAACGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTg                                                                                                                    <  2:267052/150‑1 (MQ=37)
                                gCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGtt                                                                                                                         <  1:1845053/145‑1 (MQ=255)
                                 cATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGc                                                                                                                   <  1:1028461/150‑1 (MQ=255)
                                   tctaaTCTGTAGGCCAGATAAGGCATTTTCGCCGCATCCGGCACATATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGctc                                                                                                                 <  1:121123/147‑1 (MQ=255)
                                     gAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCtt                                                                                                               >  2:6960/1‑150 (MQ=255)
                                         cTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTActt                                                                                                           >  1:1626948/1‑150 (MQ=255)
                                         cTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTActt                                                                                                           >  1:1619736/1‑150 (MQ=255)
                                           gTAGGCGAGAGATGGCCTAATCGCAGCATCCGGCACTTATTTTCGGATGCGATGATTGCGCAGCTTATGCGAGCTACAAATGGCATCGAACGGTAGGCCGGTTTGGGGGTTTACGCCGCATCCGGCAAATAGTTAAGTTCTCTTActtag                                                                                                         >  2:324836/1‑148 (MQ=255)
                                           gTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTActtct                                                                                                         >  2:151720/1‑150 (MQ=255)
                                              ggCCAGATAATGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCg                                                                                                      <  2:939853/150‑1 (MQ=255)
                                              ggCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCg                                                                                                      <  1:1238674/150‑1 (MQ=255)
                                              ggCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCg                                                                                                      >  1:1025680/1‑150 (MQ=255)
                                              ggCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCAGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCg                                                                                                      <  2:20074/150‑1 (MQ=255)
                                                ccAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGCTGCTTGCGCATCTTATCCGCCCGACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTGAATTGCTCTTACTTCTTCGcc                                                                                                    >  2:1739643/1‑150 (MQ=255)
                                                ccAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGcc                                                                                                    >  2:205857/1‑150 (MQ=255)
                                                    aTAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTg                                                                                                <  1:896586/150‑1 (MQ=255)
                                                      aaGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCa                                                                                              >  2:1360290/1‑150 (MQ=255)
                                                        ggCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAAc                                                                                            >  1:44971/1‑150 (MQ=255)
                                                        ggCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAAc                                                                                            >  1:889945/1‑150 (MQ=255)
                                                         gCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAAcg                                                                                           <  2:1092903/150‑2 (MQ=255)
                                                          cATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCa                                                                                          >  1:779688/1‑150 (MQ=255)
                                                          cATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCa                                                                                          >  1:865437/1‑150 (MQ=255)
                                                           aTTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCAc                                                                                         <  2:997707/150‑1 (MQ=255)
                                                              ttCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACttt                                                                                      >  1:469245/1‑150 (MQ=255)
                                                               tCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTta                                                                                     >  1:1400468/1‑150 (MQ=255)
                                                                cGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTtac                                                                                    <  2:1760471/150‑1 (MQ=255)
                                                                cGCAGCATACGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCAGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCATCGCCTCTCCAACCACtgtcc                                                                                    >  2:1637132/1‑146 (MQ=25)
                                                                     caTCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTTTCCGACCTACAAATCGCATCGCTCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATTGTTAATTGCGGTTACTTCTTCGCCTCTGCAACCGCTGTACTAccc                                                                               >  2:1170407/1‑150 (MQ=25)
                                                                        ccGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTTCTACCCAcg                                                                            <  1:1045898/150‑1 (MQ=255)
                                                                        ccGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCAcg                                                                            >  2:990629/1‑150 (MQ=255)
                                                                        ccGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTAGCCTCTGCAACCACTTTACTACCCAcg                                                                            >  2:1174907/1‑150 (MQ=255)
                                                                           gCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAACTCGCATCGAACCGTAGGCCGGAGAAGGCGTGTACGCCGCATGCCGGAAATAGGTAATTGCTCGTACGTCTTCGCCTCTGCAGCCACTATCCGACCCAcgccg                                                                         >  1:289539/1‑150 (MQ=21)
                                                                             aCTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCcgcg                                                                       <  1:1434340/150‑1 (MQ=255)
                                                                              cTTATTGTCGGATGCGATGCTTTCGCATCTTATCAGACCTACAAATCGCATCGAACCGTAGGCCCGATAAGGCGTTTACGCCGCCTGCGGCAAATAGATATTTGCTATTACTTCTTCGAATATGCAACCACTTTACTACCCACGCCGCgg                                                                      >  2:1335087/1‑150 (MQ=25)
                                                                                tATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCAACGCCGCgggt                                                                    >  2:61312/1‑148 (MQ=255)
                                                                                       ggATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTAt                                                             <  1:19385/150‑1 (MQ=255)
                                                                                          tgcgatgcTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGTCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTcc                                                          <  1:1037663/150‑1 (MQ=255)
                                                                                            ggatgcTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGTATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCca                                                        <  2:24832/149‑1 (MQ=255)
                                                                                            cgatgcTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCACCCACTTTACTACCCACGCCGCGGTTATTGTATTCCca                                                        >  1:684536/1‑150 (MQ=255)
                                                                                              atgcTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATTAGGCGTTTACGCCGGATCCGGCAAATAGTTAATTGCTCTTACGTCTTCGCCTCTGCAACCACTGTCCTACCCACGCCGCGGTTATTGTATTcacaca                                                      >  1:321544/1‑150 (MQ=38)
                                                                                               tgcTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACAt                                                     >  1:1428984/1‑150 (MQ=255)
                                                                                                 cTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTGATTGTATTCCCACATGc                                                   >  1:36449/1‑150 (MQ=255)
                                                                                                   tGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCgg                                                 <  1:373333/150‑1 (MQ=255)
                                                                                                   tGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCgg                                                 <  1:1081255/150‑1 (MQ=255)
                                                                                                   tGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTAGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCgg                                                 >  1:687289/1‑150 (MQ=255)
                                                                                                     cgcATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGtt                                               >  2:755786/1‑150 (MQ=255)
                                                                                                     cgcATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGtt                                               <  2:54200/150‑1 (MQ=255)
                                                                                                     cgcATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGtt                                               >  2:1120601/1‑150 (MQ=255)
                                                                                                      gcATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCTGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTg                                              <  2:718509/150‑1 (MQ=255)
                                                                                                      gcATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTg                                              <  2:1639343/150‑1 (MQ=255)
                                                                                                      gcATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTg                                              <  2:1432852/150‑1 (MQ=255)
                                                                                                         tCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTg                                                                >  2:1257636/1‑129 (MQ=255)
                                                                                                         tCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTg                                                                <  1:1257636/129‑1 (MQ=255)
                                                                                                           ttATTCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGtagtt                                         >  1:809931/1‑150 (MQ=255)
                                                                                                              tCCGACCTACAAATCGCATCGTACCGTAGGCCGGATAAGGCGTTTACGCCGCATCAGTAAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGtagttagt                                      <  1:1798481/150‑1 (MQ=255)
                                                                                                              tCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGtagttagt                                      <  1:999630/150‑1 (MQ=255)
                                                                                                                cGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAAACACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAgtgt                                    >  2:1068480/1‑150 (MQ=255)
                                                                                                                cGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCGTACTTCTTCGCCTCTGCACCCACTTTACTACCCCCGCCGCGGTTATTGTATTCCCACATCCGGTTGTAGGTAgtgt                                    >  1:1636766/1‑150 (MQ=39)
                                                                                                                     tACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTc                               <  2:195991/150‑1 (MQ=255)
                                                                                                                     tACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTc                               <  1:1787355/150‑1 (MQ=255)
                                                                                                                     tACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTc                               <  2:426690/150‑1 (MQ=255)
                                                                                                                      aCAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTGCCACATGCGGTTGTAGTTAGTGTCATTCa                              <  2:881006/150‑1 (MQ=255)
                                                                                                                      aCAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGATATTGTATTCCCACATGCGGTTGTATTTAGTGTCATTCa                              >  1:180270/1‑150 (MQ=255)
                                                                                                                        aaaTCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGa                            >  1:71451/1‑150 (MQ=255)
                                                                                                                          aTCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGGTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGAtt                          >  1:1137043/1‑150 (MQ=255)
                                                                                                                           tCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTg                         <  2:44971/150‑1 (MQ=255)
                                                                                                                             gCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTgcg                       >  1:466599/1‑150 (MQ=255)
                                                                                                                               aTCGAACCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCt                     >  2:540147/1‑150 (MQ=255)
                                                                                                                                    actatttgCCGGATGCGGCGTAAACGCCTTATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTAtt                >  2:570305/8‑150 (MQ=38)
                                                                                                                                       gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTtcg             <  1:959374/150‑1 (MQ=255)
                                                                                                                                        tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTtcgt            <  2:976774/150‑1 (MQ=255)
                                                                                                                                          ggCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTtcgtcg          <  2:329155/150‑1 (MQ=255)
                                                                                                                                          ggCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTAGTAGTTAGTGTCATTCAGATTGACCTGTATTtcgtcg          >  2:107629/1‑150 (MQ=255)
                                                                                                                                          ggCCGGATAAGGCGTTTACGCCGCATCCGGCAAAGAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGGGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTACGATTGCGCTGTATTtcgtcg          >  2:335000/1‑150 (MQ=255)
                                                                                                                                            ccGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGtt        <  2:469245/150‑1 (MQ=255)
                                                                                                                                              ggATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGTTat      >  2:1662549/1‑150 (MQ=255)
                                                                                                                                                  aaGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGTTatcatc  >  2:1763746/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GATGCTTGCGCATCTTATCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGTTATCATC  >  NC_000913/805813‑806108

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: