Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,259,550 A→T 5.1% intergenic (‑254/+169) psuK ← / ← fruA pseudouridine kinase/fused fructose‑specific PTS enzymes: IIBcomponent/IIC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,259,5500AT 5.1% 394.4 / 6.5 129intergenic (‑254/+169)psuK/fruApseudouridine kinase/fused fructose‑specific PTS enzymes: IIBcomponent/IIC components
Reads supporting (aligned to +/- strand):  ref base A (61/61);  new base T (2/5);  total (63/66)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 134 total aligned reads displayed.

TATGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCT  >  NC_000913/2259401‑2259698
                                                                                                                                                     |                                                                                                                                                    
tATGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGaaa                                                                                                                                                      <  1:829090/150‑1 (MQ=255)
         tttACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGc                                                                                                                                             <  1:1251058/150‑1 (MQ=255)
          ttACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGcc                                                                                                                                            <  2:188037/150‑1 (MQ=255)
            aCGTAAAGTAATTTACTGCCAGAGAGGCAAACTGGGCAGATTTCAGCAGGCTCACATAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGTGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTTCCCt                                                                                                                                          <  2:148324/150‑1 (MQ=255)
            aCGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCt                                                                                                                                          >  1:1006568/1‑150 (MQ=255)
             cGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCtt                                                                                                                                         <  2:652782/150‑1 (MQ=255)
                aaagtaaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAg                                                                                                                                      <  2:1181115/150‑1 (MQ=255)
                 aagtaaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGATTATTTTCTGCCCtttttt                                                                                                                                     <  2:183978/150‑5 (MQ=255)
                 aagtaaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAgg                                                                                                                                     >  1:1480742/1‑150 (MQ=255)
                  agtaaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGa                                                                                                                                    <  1:1014634/150‑1 (MQ=255)
                     aaaTTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAAt                                                                                                                                 >  2:168492/1‑150 (MQ=255)
                          aCCGCCAGAGAGGCAAAATGTGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCt                                                                                                                            >  2:1002967/1‑150 (MQ=255)
                           ccGCCAGAGAGGCAAAATGGGTAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTg                                                                                                                           >  2:971518/1‑150 (MQ=255)
                           ccGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTg                                                                                                                           >  1:800146/1‑150 (MQ=255)
                           ccGCCAGAGAGGCAAAATGGGCAGATTTGAGAAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGATGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTg                                                                                                                           <  1:1545597/150‑1 (MQ=255)
                             gCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAg                                                                                                                         >  1:1069176/1‑150 (MQ=255)
                             gCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAg                                                                                                                         >  1:201604/1‑150 (MQ=255)
                              ccAGATAGTCAAAATGGGCAGATTTGAGCATGCTCACAAAATGCGATTAAATGCATGAATTACAAAAGCGAGGAATTTGTGACTGGATGCACAGAACATAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAgt                                                                                                                        <  1:1345367/150‑1 (MQ=255)
                                   gagGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtg                                                                                                                   <  1:955436/150‑1 (MQ=255)
                                   gagGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtg                                                                                                                   <  1:1009551/150‑1 (MQ=255)
                                     gtcAATATGGTCATATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAATCGAGGAATTTGTGACTGGATGCACAGAACAGAACATTAAAAAAGTGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtgat                                                                                                                 <  1:917741/148‑1 (MQ=255)
                                     ggCAAATTGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAATTAACCTGAGTGTTtgat                                                                                                                 <  2:1424345/150‑1 (MQ=255)
                                     ggCAAAATGGGCAGATTTGATCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtgat                                                                                                                 <  2:853833/150‑1 (MQ=255)
                                     ggCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtgat                                                                                                                 >  2:1636924/1‑150 (MQ=255)
                                      gCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGATGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtgatg                                                                                                                <  1:734908/150‑1 (MQ=255)
                                       cAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATTATCTGCCCTTGAGGAAATAACCTGAGTGTTtgatga                                                                                                               >  2:112661/1‑150 (MQ=255)
                                       cAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtgatga                                                                                                               <  1:1608634/150‑1 (MQ=255)
                                       cAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTtgatga                                                                                                               >  2:645890/1‑150 (MQ=255)
                                         aaaTGGGCAGATTTGAGCAGTCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAAc                                                                                                             <  1:1486586/150‑1 (MQ=255)
                                             gggCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGttt                                                                                                         >  1:1046423/1‑150 (MQ=255)
                                              gtcAGATTTGAGCAGGCTCACAAAACGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTg                                                                                                        <  2:293435/148‑1 (MQ=255)
                                              ggCAGATTTGAGCAGGCTCACAAACTGCGCTTACATGCCTGAATTACAAAAGCGAGGAATTTGTGCCGGGATGCCCAGAACAGAACAATAAAAACTGGCAGAAAATAATCTCCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGgttt                                                                                                        >  1:57354/1‑149 (MQ=255)
                                              ggCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATACTCTGCCCTTGAGTAAATAAACTGAGTGTTTGATGAACGTTTg                                                                                                        >  1:1223809/1‑150 (MQ=255)
                                              ggCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTg                                                                                                        >  1:1679632/1‑150 (MQ=255)
                                               gCATATTTGAGCAGGCTCACAAAATGCGATTAAATGACTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGa                                                                                                       <  1:1038066/150‑1 (MQ=255)
                                                 aGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAg                                                                                                     <  2:1158301/150‑1 (MQ=255)
                                                  gATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGTGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGt                                                                                                    >  1:1380586/1‑150 (MQ=255)
                                                     ttGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTc                                                                                                 <  2:1679632/150‑1 (MQ=255)
                                                     ttGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTc                                                                                                 <  1:680012/150‑1 (MQ=255)
                                                               tCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGATGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGa                                                                                       >  2:285590/1‑150 (MQ=255)
                                                                cacaAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGAt                                                                                      >  2:240028/1‑150 (MQ=255)
                                                                  caAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCGTGAGGAAATAACCTGAGTGTTTGTTGAACGTTTGAAGTTTCTACTGAAAGATaa                                                                                    <  1:913386/150‑1 (MQ=255)
                                                                   aaaaTGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATaaa                                                                                   <  1:70482/150‑1 (MQ=255)
                                                                    aaaTGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATACCATCAGTGTTTGATGAACTTTTGAAGTTTCTACTGAAAGATaaaa                                                                                  >  1:923886/1‑150 (MQ=255)
                                                                     aaTGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAAc                                                                                 <  1:1636924/150‑1 (MQ=255)
                                                                        gCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCAt                                                                              >  1:381787/1‑150 (MQ=255)
                                                                           aTTAAATGCCTGAATTACAAAAGCGTGGAATTTGTGACGGGATGCACATAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTa                                                                           <  1:1583089/150‑1 (MQ=255)
                                                                              aaaTGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAAt                                                                        >  2:1259854/1‑150 (MQ=255)
                                                                              aaaTGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAAt                                                                        >  2:384923/1‑150 (MQ=255)
                                                                                 tGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCa                                                                     <  2:1604414/150‑1 (MQ=255)
                                                                                   ccTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATa                                                                   >  1:1220829/1‑150 (MQ=255)
                                                                                   ccTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATa                                                                   >  2:68563/1‑150 (MQ=255)
                                                                                    cTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATaa                                                                  >  2:493330/1‑150 (MQ=255)
                                                                                        aTTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAAtt                                                              <  2:1069176/150‑1 (MQ=255)
                                                                                         ttACAAAAGCGAGGAATTTGTGACGGGTTGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTa                                                             >  2:1345527/1‑150 (MQ=255)
                                                                                         ttACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTa                                                             <  1:1714095/150‑1 (MQ=255)
                                                                                           aCAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGACCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAAACATATAAATACATAAAATTACa                                                           >  2:680346/1‑150 (MQ=255)
                                                                                            cAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACAt                                                          <  2:1046423/150‑1 (MQ=255)
                                                                                             aaaaGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACAtt                                                         >  1:201077/1‑150 (MQ=255)
                                                                                              aaaGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTc                                                        >  1:869514/1‑150 (MQ=255)
                                                                                               aaGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAATAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTca                                                       <  1:168492/150‑1 (MQ=255)
                                                                                                 gCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTcaca                                                     <  1:1456459/150‑1 (MQ=255)
                                                                                                 gCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTcaca                                                     >  2:1220082/1‑150 (MQ=255)
                                                                                                    aGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAAc                                                  >  2:862726/1‑150 (MQ=255)
                                                                                                     ggAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAAcc                                                 <  2:1220829/150‑1 (MQ=255)
                                                                                                     ggAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAAcc                                                 >  2:1176597/1‑150 (MQ=255)
                                                                                                         tttGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTc                                             <  2:1719536/150‑1 (MQ=255)
                                                                                                            gtgACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCAtt                                          >  2:666206/1‑150 (MQ=255)
                                                                                                             tgACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATtc                                         <  2:1342308/150‑1 (MQ=255)
                                                                                                             tgACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATtc                                         <  1:870311/150‑1 (MQ=255)
                                                                                                              gACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCa                                                                     <  2:612713/121‑1 (MQ=255)
                                                                                                              gACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCa                                                                     >  1:612713/1‑121 (MQ=255)
                                                                                                               aCGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATtctc                                       >  2:1567042/1‑150 (MQ=255)
                                                                                                                cGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATATCCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCAtgcgct                                      <  2:67823/150‑6 (MQ=255)
                                                                                                                  ggATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCttt                                    >  2:622265/1‑150 (MQ=255)
                                                                                                                   gATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTg                                   <  2:201077/150‑1 (MQ=255)
                                                                                                                    aTGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTgc                                  <  1:43541/150‑1 (MQ=255)
                                                                                                                      gCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATACCCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTgcgc                                >  1:83066/1‑150 (MQ=255)
                                                                                                                       cacaGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTg                                   <  1:73428/146‑1 (MQ=255)
                                                                                                                       cacaGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTg                                   >  2:73428/1‑146 (MQ=255)
                                                                                                                       cacaGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCt                               >  1:1216984/1‑150 (MQ=255)
                                                                                                                       cacaGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCt                               <  1:428359/150‑1 (MQ=255)
                                                                                                                       cacaGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCt                               >  1:217149/1‑150 (MQ=255)
                                                                                                                       cacaGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCt                               >  1:1267346/1‑150 (MQ=255)
                                                                                                                          agaacagatcAATAAAAAATGGCAGAAAATAATCTGCCCTTTAGGAAATAACCTGAGTGTTTTATGAACGTTTGAAGTTTCTAGTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCtcat                            <  1:680346/150‑1 (MQ=255)
                                                                                                                           gaacagaacaATAAAAAAGGGCAGAAAATATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCtcatca                          >  2:948550/1‑150 (MQ=255)
                                                                                                                              cagaacaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCaaa                        >  1:595349/1‑150 (MQ=255)
                                                                                                                                 aacaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGt                     <  1:971518/150‑1 (MQ=255)
                                                                                                                                  acaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGtt                    >  1:1640826/1‑150 (MQ=255)
                                                                                                                                  acaATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGtt                    >  2:1611950/1‑150 (MQ=255)
                                                                                                                                   caATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAAACTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTa                   <  2:835433/150‑1 (MQ=255)
                                                                                                                                       aaaaaaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacagg               <  1:1567042/150‑1 (MQ=255)
                                                                                                                                        ctcaaGGGCAGATTATTTTCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacagga              >  1:1363068/4‑150 (MQ=255)
                                                                                                                                          caaGGGCAGATTATTTTCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacaggaca            >  1:309802/2‑150 (MQ=255)
                                                                                                                                           aaGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacaggacag           >  2:172300/1‑150 (MQ=255)
                                                                                                                                           aaGGGCAGAAAATAATCTGCACTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTacaggacag           <  1:581244/150‑1 (MQ=255)
                                                                                                                                              ggCAGAAAATAATCTGCCCTTGAGGAAATAACATGACTGAATGATGAACGTTAGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATATAATTACATTCACAAAACGGTCATTCTCTTTGCGCTCATCAAAAGTTACATAACAGGaa        >  1:302963/1‑150 (MQ=255)
                                                                                                                                                  gAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAAttt    >  1:474788/1‑150 (MQ=255)
                                                                                                                                                   aaaaTAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTc   >  1:1265823/1‑150 (MQ=255)
                                                                                                                                                    aaaTAATCTGCCCTTGAGGTAATAACCTGAGTGTTTGATGAACGTTTGTAGTTTCTACTGAAAGATAAAAGCATTTATATTCATAAAATTTCATTCACAAACCGGTCTTTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCt  <  2:302963/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TATGCATGTTTTACGTAAAGTAAATTACCGCCAGAGAGGCAAAATGGGCAGATTTGAGCAGGCTCACAAAATGCGATTAAATGCCTGAATTACAAAAGCGAGGAATTTGTGACGGGATGCACAGAACAGAACAATAAAAAAGGGCAGAAAATAATCTGCCCTTGAGGAAATAACCTGAGTGTTTGATGAACGTTTGAAGTTTCTACTGAAAGATAAAACCATTTAAATTCATAAAATTACATTCACAAACCGGTCATTCTCTTTGCGCTCATCAAATGTTACAGGACAGGAAATTTCT  >  NC_000913/2259401‑2259698

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: