Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,561,969 A→T 6.6% intergenic (‑146/‑44) glpE ← / → glpD thiosulfate sulfurtransferase GlpE/aerobic glycerol 3‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,561,9690AT 6.6% 337.3 / 10.4 106intergenic (‑146/‑44)glpE/glpDthiosulfate sulfurtransferase GlpE/aerobic glycerol 3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (56/43);  new base T (5/2);  total (61/45)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 106 total aligned reads displayed.

CATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGC  >  NC_000913/3561821‑3562115
                                                                                                                                                    |                                                                                                                                                  
cATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGAt                                                                                                                                                   >  1:1278200/1‑150 (MQ=255)
   gCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAAc                                                                                                                                                >  1:1166726/1‑150 (MQ=255)
       ctctctTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAAc                                                                                                                                            >  2:416412/1‑150 (MQ=255)
       ctctctTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAAc                                                                                                                                            >  2:1678442/1‑150 (MQ=255)
         ctctTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAt                                                                                                                                          >  1:462007/1‑150 (MQ=255)
          tctttcttTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGTTATCGAACAta                                                                                                                                         <  2:441552/150‑2 (MQ=255)
          tctttcttTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAtt                                                                                                                                         >  1:772403/1‑150 (MQ=255)
            tttctttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTa                                                                                                                                       <  1:955883/150‑1 (MQ=255)
               ctttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGTTATCGAACAtattagc                                                                                                                                    <  1:521277/150‑7 (MQ=255)
               ctttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGa                                                                                                                                    >  1:813772/1‑150 (MQ=255)
                 ttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGc                                                                                                                                  <  2:335471/150‑1 (MQ=255)
                 ttACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGc                                                                                                                                  >  2:1290102/1‑150 (MQ=255)
                   acaaacaaGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCtt                                                                                                                                >  2:557760/1‑150 (MQ=255)
                    caaacaaGTGGGCAAATTTACCGCACAGTTTACGTAGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCttt                                                                                                                               >  1:331495/1‑150 (MQ=255)
                     aaacaaGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTa                                                                                                                              <  2:171615/150‑1 (MQ=255)
                     aaacaaGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTa                                                                                                                              <  2:1278200/150‑1 (MQ=255)
                        caaGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACg                                                                                                                           <  1:1258742/150‑1 (MQ=255)
                              ggCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTg                                                                                                                     <  1:1922302/150‑1 (MQ=255)
                                cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGag                                                                                                                   >  2:306746/1‑149 (MQ=255)
                                cAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGaa                                                                                                                   >  1:1347505/1‑150 (MQ=255)
                                 aaaTTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAAt                                                                                                                  <  1:1158340/150‑1 (MQ=255)
                                  aaTTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATg                                                                                                                                     >  1:114996/1‑130 (MQ=255)
                                  aaTTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATg                                                                                                                                     <  2:114996/130‑1 (MQ=255)
                                    tttACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAg                                                                                                               <  1:828290/150‑1 (MQ=255)
                                         cGCACAGTTTACGTCGAAGCTGCAGATAAACGCCATAATGTTATACATATCACCCTAAAATTTTTTTTCAATGTTACCTAAAGCGCTATTCTTTGGTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGgcag                                                                                                          <  2:1101825/150‑1 (MQ=255)
                                           cacaGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGgcagca                                                                                                        <  2:1307525/150‑1 (MQ=255)
                                                 ttACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGaaa                                                                                                  <  1:1678442/150‑1 (MQ=255)
                                                  tACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAAc                                                                                                 >  1:1561842/1‑150 (MQ=255)
                                                    cGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACTAAAGTGAATGAGGGCAGCATGGAAACCa                                                                                               >  1:1260843/1‑150 (MQ=255)
                                                     gTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCaa                                                                                              >  1:172040/1‑150 (MQ=255)
                                                     gTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCaa                                                                                              >  2:1559190/1‑150 (MQ=255)
                                                      tCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCaaa                                                                                             >  2:671102/1‑150 (MQ=255)
                                                       cGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAg                                                                                            <  2:1166726/150‑1 (MQ=255)
                                                        gAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGa                                                                                           >  1:1331363/1‑150 (MQ=255)
                                                            cGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTg                                                                                       >  1:625910/1‑150 (MQ=255)
                                                             ggCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGa                                                                                      >  2:1716336/1‑150 (MQ=255)
                                                              gCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGAt                                                                                     >  1:1517826/1‑150 (MQ=255)
                                                               cAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGAtt                                                                                    >  1:191504/1‑150 (MQ=255)
                                                                aGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATtg                                                                                   <  2:245935/150‑1 (MQ=255)
                                                                aGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATtg                                                                                   <  1:1299004/150‑1 (MQ=255)
                                                                  aTAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATtgtg                                                                                 <  2:625910/150‑1 (MQ=255)
                                                                  aTAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATtgtg                                                                                 >  1:1795647/1‑150 (MQ=255)
                                                                      aCGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAg                                                                             >  1:376544/1‑150 (MQ=255)
                                                                       cGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAgg                                                                            >  2:812715/1‑150 (MQ=255)
                                                                           aTAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGggc                                                                        <  2:479510/150‑1 (MQ=255)
                                                                           aTAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGggc                                                                        <  1:928153/150‑1 (MQ=255)
                                                                            tAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGggcg                                                                       >  1:163783/1‑150 (MQ=255)
                                                                             aaTGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGggcgg                                                                      <  1:1683164/150‑1 (MQ=255)
                                                                              aTGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGggcggc                                                                     >  1:1376976/1‑150 (MQ=255)
                                                                              aTGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGggcggc                                                                     <  2:1435444/150‑1 (MQ=255)
                                                                                   ataCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCaa                                                                >  2:337933/1‑150 (MQ=255)
                                                                                     aCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATg                                                              >  2:1861353/1‑150 (MQ=255)
                                                                                     aCATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATg                                                              <  1:1619193/150‑1 (MQ=255)
                                                                                             ctctAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGt                                                      >  1:1160534/1‑150 (MQ=255)
                                                                                              tctAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTa                                                     <  2:772403/150‑1 (MQ=255)
                                                                                                tAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATc                                                   >  1:1400338/1‑150 (MQ=255)
                                                                                                 aaaaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATcg                                                  >  2:1811339/1‑150 (MQ=255)
                                                                                                  aaaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATcgc                                                 >  1:99903/1‑150 (MQ=255)
                                                                                                   aaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATgg                                                             <  2:844304/137‑1 (MQ=255)
                                                                                                   aaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATgg                                                             >  1:844304/1‑137 (MQ=255)
                                                                                                   aaTGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATcgcg                                                <  1:1559190/150‑1 (MQ=255)
                                                                                                     tGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGc                                              <  2:1347505/150‑1 (MQ=255)
                                                                                                      gTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATtgtg                                                                                 >  2:462331/1‑114 (MQ=255)
                                                                                                      gTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATtgtg                                                                                 <  1:462331/114‑1 (MQ=255)
                                                                                                      gTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCa                                             >  1:921198/1‑150 (MQ=255)
                                                                                                      gTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCa                                             >  1:1304345/1‑150 (MQ=255)
                                                                                                       ttttttCAATGTTACCTAAAGCTCGATTATTTGCTAATATTTTCGATAACGAACATTTATGAGATTTAACGAAAGTGAATGAGGGCAGAATGGAAACCAAAGATCTGATTGTGATATGGGGCGGCATCAATGGTGCTGGTATAGCGGCAg                                            >  2:504876/1‑150 (MQ=255)
                                                                                                       ttttttCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAg                                            <  1:1716336/150‑1 (MQ=255)
                                                                                                       ttttttCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGATGGTATCGCGGCAg                                            >  2:512146/1‑150 (MQ=255)
                                                                                                          tttCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGAcg                                         >  2:777939/1‑150 (MQ=255)
                                                                                                           ttCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGAcgc                                        <  1:591718/150‑1 (MQ=255)
                                                                                                            tCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAATATCTGATTGTGATAGGGGTCGTCATCAATGGTGCTGGTATCGCTGGAGAcgcc                                       <  2:607723/150‑1 (MQ=255)
                                                                                                                tGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTg                                   <  2:331495/150‑1 (MQ=255)
                                                                                                                  ttACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGa                                 <  1:1363719/150‑1 (MQ=255)
                                                                                                                   tACCTAAAGCGCGATTCTTTTCTAATATCTTCGATAACCAACATTTATGAGCTTTAACGAAAGAGAATGAGGGCAGCATGGAAACCAAAGATATGATTTTGATAGGTCTCCGCATCAATGGTGCTGGTATCGCGGCAGAAGCCGCTGGAc                                >  1:1863138/1‑150 (MQ=255)
                                                                                                                    aCCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGAcg                               <  2:813772/150‑1 (MQ=255)
                                                                                                                     ccTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGAcgc                              >  1:1699043/1‑150 (MQ=255)
                                                                                                                     ccTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGAcgc                              <  2:462007/150‑1 (MQ=255)
                                                                                                                     ccTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGAcgc                              >  1:1183199/1‑150 (MQ=255)
                                                                                                                       tAAAGCTCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATg                                                              <  2:1616815/116‑1 (MQ=255)
                                                                                                                       tAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATg                                                              >  1:1616815/1‑116 (MQ=255)
                                                                                                                        aaaGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGt                           <  2:1376976/150‑1 (MQ=255)
                                                                                                                        aaaGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGt                           >  2:1811566/1‑150 (MQ=255)
                                                                                                                          aGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGttt                         >  1:977525/1‑150 (MQ=255)
                                                                                                                           gcgcGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTa                        <  2:191504/150‑1 (MQ=255)
                                                                                                                           gcgcGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTa                        <  2:1517826/150‑1 (MQ=255)
                                                                                                                              cgATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATcc                     >  1:217028/1‑150 (MQ=255)
                                                                                                                                    tttGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTg               >  1:593629/1‑150 (MQ=255)
                                                                                                                                      ctcataaatGTTCGTTATCGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGAt             >  2:1380771/8‑150 (MQ=255)
                                                                                                                                      ctcataaatGTTCGTTATCGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGAt             >  2:1465801/8‑150 (MQ=255)
                                                                                                                                      ctcataaatGTTCGTTATCGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGAt             >  2:1472757/8‑150 (MQ=255)
                                                                                                                                        cataaatGTTCGTTATCGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGc           >  2:1779139/6‑150 (MQ=255)
                                                                                                                                        cTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGc           >  1:302452/1‑150 (MQ=255)
                                                                                                                                           atatGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTgg        <  2:1331363/150‑1 (MQ=255)
                                                                                                                                            tatGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGa       >  2:355904/1‑150 (MQ=255)
                                                                                                                                            aatGTTCGTTATCGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGa       >  1:1169957/2‑150 (MQ=255)
                                                                                                                                             atGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAg      >  2:1889176/1‑150 (MQ=255)
                                                                                                                                             atGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAg      <  2:376544/150‑1 (MQ=255)
                                                                                                                                              tGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAgg     >  2:811300/1‑150 (MQ=255)
                                                                                                                                                 tCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGgcgc  <  2:1282698/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CATGCGTCTCTCTTTCTTTACAAACAAGTGGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAACGAAAGTGAATGAGGGCAGCATGGAAACCAAAGATCTGATTGTGATAGGGGGCGGCATCAATGGTGCTGGTATCGCGGCAGACGCCGCTGGACGCGGTTTATCCGTGCTGATGCTGGAGGCGC  >  NC_000913/3561821‑3562115

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