Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,146,761 T→G 5.1% intergenic (‑24/+95) yghW ← / ← yghX DUF2623 family protein/pseudogene, dienlactone hydrolase family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,146,7610TG 5.1% 413.9 / 9.7 138intergenic (‑24/+95)yghW/yghXDUF2623 family protein/pseudogene, dienlactone hydrolase family
Reads supporting (aligned to +/- strand):  ref base T (55/76);  new base G (5/2);  total (60/78)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

Only 100 of 141 total aligned reads displayed.

CGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGC  >  NC_000913/3146613‑3146909
                                                                                                                                                    |                                                                                                                                                    
cGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTg                                                                                                                                                     <  1:447392/150‑1 (MQ=255)
cGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTg                                                                                                                                                     <  1:1180919/150‑1 (MQ=255)
 ggTGTGAGTATCCTAATACAAAACCTCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGt                                                                                                                                                    <  2:1951682/150‑1 (MQ=255)
 ggTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCCCCTCATTCGCTATATATTGt                                                                                                                                                    >  1:1564875/1‑150 (MQ=255)
 ggTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACAATAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGt                                                                                                                                                    <  1:580409/150‑1 (MQ=255)
     tgAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGAACTGTCGGCATGCGTTGCTATTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCata                                                                                                                                                >  1:306441/1‑150 (MQ=255)
           gCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTTGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATGACAATATATAGc                                                                                                                                          >  1:451005/1‑150 (MQ=255)
            ccTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCg                                                                                                                                         >  1:1867246/1‑150 (MQ=255)
             cTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAgaga                                                                                                                                        >  1:1144735/1‑150 (MQ=255)
             cTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGa                                                                                                                                        >  1:1143419/1‑150 (MQ=255)
              tAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGat                                                                                                                                       >  2:1497759/1‑150 (MQ=255)
                  aCAAAACCGCGTTAATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATAttt                                                                                                                                   <  2:259469/150‑1 (MQ=255)
                   cAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATAtttt                                                                                                                                  <  2:1173523/150‑1 (MQ=255)
                     aaaCCGCGTTTATAATCGGCAGAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAg                                                                                                                                <  2:1879467/150‑1 (MQ=255)
                     aaaCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAg                                                                                                                                <  1:1646905/150‑1 (MQ=255)
                      aaCCGCGTTTATAATCGGCACAGAATTTTGTAACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAATCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAgg                                                                                                                               <  1:1334480/150‑1 (MQ=255)
                      aaCCGCGTTTATAATAGGCACAGAATTTTGTCACATTAACCGCACGGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAgg                                                                                                                               <  2:295046/150‑1 (MQ=255)
                       aCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGc                                                                                                                              >  1:1952759/1‑150 (MQ=255)
                       aCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGc                                                                                                                              <  2:1558555/150‑1 (MQ=255)
                       aCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGc                                                                                                                              <  2:1031304/150‑1 (MQ=255)
                         cgcgTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGACaa                                                                                                                            >  1:812332/1‑150 (MQ=255)
                          gcgTTTATAATCGGCACAGAATTTTGTCACATTAACCTCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCaaa                                                                                                                           <  1:597551/150‑1 (MQ=255)
                           cgTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAg                                                                                                                          <  2:1542323/150‑1 (MQ=255)
                             tttATAATCGGCACAGAATATTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTAATATATAGCGATATTTTAGGCAAAGaa                                                                                                                        >  2:1905964/1‑150 (MQ=255)
                              ttATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAAc                                                                                                                       >  2:832288/1‑150 (MQ=255)
                               tataATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCAGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACg                                                                                                                      <  2:1017843/150‑1 (MQ=255)
                                  aaTCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACt                                                                                                                   >  2:1343134/1‑150 (MQ=255)
                                  aaTCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACt                                                                                                                   <  1:1068040/150‑1 (MQ=255)
                                    tCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACttt                                                                                                                 <  1:693689/150‑1 (MQ=255)
                                      ggCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTg                                                                                                               >  2:1115167/1‑150 (MQ=255)
                                       gCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCGATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTgc                                                                                                              >  2:1541542/1‑150 (MQ=255)
                                        cacaGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTgcg                                                                                                             >  2:469466/1‑150 (MQ=255)
                                        cacaGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTgcg                                                                                                             >  1:99522/1‑150 (MQ=255)
                                        cacaGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTgcg                                                                                                             >  2:1168678/1‑150 (MQ=255)
                                           aGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCAc                                                                                                          <  1:738214/150‑1 (MQ=255)
                                                    tCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTAttttt                                                                                                 <  2:1952758/150‑1 (MQ=255)
                                                     cacaTTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTAtttttt                                                                                                >  2:951424/1‑150 (MQ=255)
                                                      acaTTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTc                                                                                               <  2:450908/150‑1 (MQ=255)
                                                       caTTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATAATGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCa                                                                                              >  2:147259/1‑150 (MQ=255)
                                                       caTTAACCGCACTGTCGGCATGAGTTGCTTTTAATCCCGCCATTAATCCTTTACCAAAATGGGTATTCATCGTTCACCTCATTCGCTATATATTGTCATCTATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCa                                                                                              <  1:949840/150‑1 (MQ=255)
                                                            aCCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATAAATTTTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATTCtt                                                                                         >  2:1041882/1‑150 (MQ=255)
                                                             ccGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTg                                                                                        <  1:831634/150‑1 (MQ=255)
                                                             ccGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTg                                                                                        <  2:1937320/150‑1 (MQ=255)
                                                              cGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGt                                                                                       <  2:808705/150‑1 (MQ=255)
                                                               ccACTGTCGGCATGCGGTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGATCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGtt                                                                                      <  2:1972429/149‑1 (MQ=255)
                                                                cACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGttt                                                                                     >  1:1043448/1‑150 (MQ=255)
                                                                 aCTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTa                                                                                    <  2:1867245/150‑1 (MQ=255)
                                                                  cTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTaa                                                                                   >  1:239082/1‑150 (MQ=255)
                                                                  cTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTaa                                                                                   >  2:416171/1‑150 (MQ=255)
                                                                   tGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAg                                                                                  <  2:637624/150‑1 (MQ=255)
                                                                    gTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGAGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGa                                                                                 <  1:149860/150‑1 (MQ=255)
                                                                      cGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACg                                                                               <  1:1532299/150‑1 (MQ=255)
                                                                          aTGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAAt                                                                           >  2:1679620/1‑150 (MQ=255)
                                                                               ttGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGt                                                                      <  2:162777/150‑1 (MQ=255)
                                                                                tGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTa                                                                     >  2:1519477/1‑150 (MQ=255)
                                                                                 gcgtttAATCCCGCCATTAAGACTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAAAGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTTAGACGTAATCCCGaaa                                                                    >  1:1096406/4‑150 (MQ=255)
                                                                                  cTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAag                                                                   >  2:1182766/1‑150 (MQ=255)
                                                                                  cTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGAGCACGCTTTTTTTTCATGCTTGTTTAAGACGTAGTCCCGTatg                                                                   >  1:175380/1‑148 (MQ=255)
                                                                                  cTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATAATTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAag                                                                   >  1:1153359/1‑150 (MQ=255)
                                                                                     ttAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAgaga                                                                <  2:230423/150‑1 (MQ=255)
                                                                                      tAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGAt                                                               <  1:1526259/150‑1 (MQ=255)
                                                                                       aaTCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATAACGTAAGAGATa                                                              >  1:1327755/1‑150 (MQ=255)
                                                                                       aaTCCAGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATa                                                              <  2:1849253/150‑1 (MQ=255)
                                                                                          cccGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGCACGCCTGTTGCGCACGCGATTTTTTCAGGCGTGTGGATGACGTATTCCGGTATGAGTTaaag                                                           >  2:1915520/1‑149 (MQ=255)
                                                                                           ccGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAAcc                                                          >  2:484383/1‑150 (MQ=255)
                                                                                           ccGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAAcc                                                          <  1:425876/150‑1 (MQ=255)
                                                                                                aTTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCa                                                     <  1:835635/150‑1 (MQ=255)
                                                                                                  tAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAgg                                                   <  1:1168679/150‑1 (MQ=255)
                                                                                                    aGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGgag                                                 <  2:1360389/150‑1 (MQ=255)
                                                                                                      ccattACCAAAATGGTTATTCATCGTTCACCTCATTCGCTTTATATTGTCATATCTAGCGTTATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGATTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGgagag                                               <  2:390537/147‑1 (MQ=255)
                                                                                                      ccTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGgagag                                               <  2:1043447/150‑1 (MQ=255)
                                                                                                       cTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGgagaga                                              <  1:778880/150‑1 (MQ=255)
                                                                                                       cTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGgagaga                                              >  2:482904/1‑150 (MQ=255)
                                                                                                           aCCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATAtt                                          <  2:174960/150‑1 (MQ=255)
                                                                                                           aCCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGAGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATAtt                                          <  2:751116/150‑1 (MQ=255)
                                                                                                            ccAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATAttt                                         >  2:926293/1‑150 (MQ=255)
                                                                                                             cAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTa                                        <  2:1591900/150‑1 (MQ=255)
                                                                                                              aaaaTGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTAt                                       <  1:1802912/150‑1 (MQ=255)
                                                                                                               aaaTGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATc                                      <  1:1686239/150‑1 (MQ=255)
                                                                                                                   ggTTATTCATCGTTCCCCTCATTCGCTATAGATTGTCATATATAGCGATATTTTAGGCAAAGAACGAGTTTAGCGCACGCTATTTTTTCATCCTTGTGTAAGACGTAATCCCGTAAGAGAAAAACCTATCAGGAGAGATATTTATCGAAc                                  >  2:332806/1‑150 (MQ=255)
                                                                                                                    gTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATAAAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAAACCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAAcc                                 <  1:1109632/150‑1 (MQ=255)
                                                                                                                       aTTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGCATTTTGCTCTCGCTATTTTTTCATACTTGTTTAAGACGTAATACCTTAAGAGTTAAACCTATCAGGAAAAATATTTATAGAACAAtt                              >  1:352960/1‑150 (MQ=255)
                                                                                                                           aTCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAg                          <  1:1764010/150‑1 (MQ=255)
                                                                                                                             cGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCTTGCTTGTTTAAGACGTAATACCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAgtg                        >  1:1282105/1‑150 (MQ=255)
                                                                                                                                   ccTCATTCGCTATATATTGTCAGATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTTAGACGTAATCCCGTAAGAGATAACCCTATCAGGTGAGATATTTATAGAACCATTTCAGTGTCCttt                  >  1:1431051/1‑150 (MQ=255)
                                                                                                                                    cTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTg                 <  1:962377/150‑1 (MQ=255)
                                                                                                                                    cTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTg                 <  2:1988129/150‑1 (MQ=255)
                                                                                                                                    cTCATTCGCTATATATTGTCATATATAGCGATATTTTAGACAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTg                 <  2:812331/150‑1 (MQ=255)
                                                                                                                                      cATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTAATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGcc               >  1:896850/1‑150 (MQ=255)
                                                                                                                                        atCGCTATATATGACAATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAg             >  2:610577/2‑150 (MQ=255)
                                                                                                                                        atCGCTATATATGACAATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAg             >  1:1646489/2‑150 (MQ=255)
                                                                                                                                          cGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGc           >  2:1241155/1‑150 (MQ=255)
                                                                                                                                          cGCTATATATGACAATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGc           >  1:1588072/1‑150 (MQ=255)
                                                                                                                                           gCTATATATGACAATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGTCa          >  1:409947/1‑150 (MQ=255)
                                                                                                                                            caatatatTGTCATATAAAGCGATATTTTAGGCAAAGAACTACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACTTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCaa         <  2:115947/148‑1 (MQ=255)
                                                                                                                                              atatatTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGa       >  2:794083/1‑150 (MQ=255)
                                                                                                                                                 tatTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTACGACGTAATCCCGTAAGAGATAAACCTATCAGGATAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATcg    <  2:1118426/150‑1 (MQ=255)
                                                                                                                                                   ttGTCATTTATAGCGATATTTTAGGCAAAGACCGACTTTTGCGCACTCTATTTTTTCATGATTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTACTTTTCCAGGCAAGATcggc  <  1:1041883/150‑1 (MQ=255)
                                                                                                                                                   ttGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATcggc  <  1:1519478/150‑1 (MQ=255)
                                                                                                                                                   ttGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATcggc  <  1:712399/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGCTATATATTGTCATATATAGCGATATTTTAGGCAAAGAACGACTTTTGCGCACGCTATTTTTTCATGCTTGTTTAAGACGTAATCCCGTAAGAGATAAACCTATCAGGAGAGATATTTATCGAACCATTTCAGTGTCCTTTGCCAGGCAAGATCGGC  >  NC_000913/3146613‑3146909

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