Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 584,237 A→G 5.0% D186G (GAT→GGT)  appY → global transcriptional activator; DLP12 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913584,2370AG 5.0% 315.2 / 4.4 100D186G (GAT→GGT) appYglobal transcriptional activator; DLP12 prophage
Reads supporting (aligned to +/- strand):  ref base A (40/54);  new base G (2/3);  total (43/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 103 total aligned reads displayed.

CTGTAAAATTACTGGAATAATCAGTTTTAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTT  >  NC_000913/584094‑584381
                                                                                                                                               |                                                                                                                                                
cTGTAAAATTACTGGAATAATCAGTTTTAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACt                                                                                                                                              <  7:30779/148‑1 (MQ=255)
  gTAAAATTACTGGAATAATCAGTTTTAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAg                                                                                                                                            >  4:47813/1‑148 (MQ=255)
   tAAAATTACTGGAATAATCAGTTTTAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTa                                                                                                                                             <  1:212499/146‑1 (MQ=255)
   tAAAATTACTGGAATAATCAGTTTTAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTa                                                                                                                                             <  5:32556/146‑1 (MQ=255)
   tAAAATTACTGGAATAATCAGTTTTAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAgg                                                                                                                                           <  3:196646/148‑1 (MQ=255)
   tAAAATTACTGGAATAATCAGTTTTAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAg                                                                                                                                            >  3:26280/1‑147 (MQ=255)
   tAAAATTACTGGAATAATCAGTTTTAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAg                                                                                                                                            <  7:275282/147‑1 (MQ=255)
   tAAAATTACTGGAATAATCAGTTTTAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAgg                                                                                                                                           <  1:325433/147‑1 (MQ=255)
     aaaTTACTGGAATAATCATTTTTAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGa                                                                                                                                          <  4:248898/147‑1 (MQ=255)
         tACTGGAATAATCAGTTTTAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAgg                                                                                                                                     >  5:348132/1‑148 (MQ=255)
             ggAATAATCAGTTTTAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATg                                                                                                                                 >  1:375921/1‑148 (MQ=255)
              gAATAATCAGTTTTAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATg                                                                                                                                 <  7:126609/147‑1 (MQ=255)
               aataatCAGTTTTAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATg                                                                                                                                 >  5:112214/1‑146 (MQ=255)
               aataatCAGTTTTAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATg                                                                                                                                 <  6:112214/146‑1 (MQ=255)
                         tttAATATTGAGCGGCAATGGCATTTTAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATAc                                                                                                                                               <  2:311586/122‑1 (MQ=255)
                         tttAATATTGAGCGGCAATGGCATTTTAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATAc                                                                                                                                               >  1:311586/1‑122 (MQ=255)
                         tttAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGAt                                                                                                                                                 >  5:265879/1‑120 (MQ=255)
                         tttAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGAt                                                                                                                                                 <  6:265879/120‑1 (MQ=255)
                            aaTATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATa                                                                                                                   <  7:197348/147‑1 (MQ=255)
                            aaTATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATa                                                                                                                   <  6:2682/147‑1 (MQ=255)
                                 tGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTc                                                                                                              >  1:97536/1‑147 (MQ=255)
                                   aGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCaa                                                                                                            <  6:348132/147‑1 (MQ=255)
                                     cgtcAATGTCATTTAAAAGATATTGCGTAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACTTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAAct                                                                                                         <  8:167799/145‑1 (MQ=255)
                                     ccgtaaTTGCATTTAAAAGATATTGCGGAATTTATTTATACGAGTGAAATTTTAATAAAAAAAAGATTAAGGGATGAAGCAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAAct                                                                                                         <  6:120525/144‑1 (MQ=255)
                                       gCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAActct                                                                                                       <  2:97536/148‑1 (MQ=255)
                                         aaTGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTa                                                                                                     >  8:164564/1‑148 (MQ=255)
                                            ggCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAgg                                                                                                                                           >  7:370887/1‑107 (MQ=255)
                                            ggCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAgg                                                                                                                                           <  8:370887/107‑1 (MQ=255)
                                            ggCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTc                                                                                                  >  5:357107/1‑148 (MQ=255)
                                              cATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTa                                                                                                <  6:290337/148‑1 (MQ=255)
                                               aTTTAAAAGATATTGCGGAATTGATTTATACGACGGAACTGTTAATAAAAAAACGATTAAGGGATGAAGGACCGTCATTTACTGAAATATTGAGAGGTACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTAt                                                                                               <  6:143717/148‑1 (MQ=255)
                                                tttAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTAt                                                                                               >  2:225128/1‑147 (MQ=255)
                                                 ttAAAAGATATTGCGGATTTGAGTTATACGAGTGAAAGTTTATTAAAAAAGAGAGTATGGGATGAAGGGACGTCGTTTACTGAATTATTGAGAGAGACGAGGATGCGGGATGCATAATGACTCATGACTTCGAACTCTTTTtctctc                                                                                              >  4:35525/1‑147 (MQ=255)
                                                 ttAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATc                                                                                              >  8:52580/1‑147 (MQ=255)
                                                  tAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGAt                                                                                                                                                 <  7:51370/95‑1 (MQ=255)
                                                  tAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGAt                                                                                                                                                 >  8:51370/1‑95 (MQ=255)
                                                  tAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCaa                                                                                            >  3:229912/1‑148 (MQ=255)
                                                   aaaaGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCaa                                                                                            <  5:39554/147‑1 (MQ=255)
                                                   aaaaGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATACTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCaa                                                                                            <  5:370337/147‑1 (MQ=255)
                                                    aaaGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaaaaa                                                                                                                          >  4:346114/1‑116 (MQ=255)
                                                    aaaGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaaaaa                                                                                                                          <  3:346114/116‑1 (MQ=255)
                                                    aaaGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTa                                                                                                >  7:136093/1‑142 (MQ=255)
                                                    aaaGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTa                                                                                                <  8:136093/142‑1 (MQ=255)
                                                     aaGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGt                                                                                         >  6:298579/1‑148 (MQ=255)
                                                      aGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTc                                                                                        >  1:299969/1‑148 (MQ=255)
                                                              cGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGcc                                                                                  >  2:86815/1‑146 (MQ=255)
                                                              cGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCa                                                                                <  3:47813/148‑1 (MQ=255)
                                                               ggAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGACCGTCATTTAGTGAAATATTGATAGCTACTAGGATGAGGTATGCAAAAACACTCATAACTTCAAACTCTTATTCTCTCAATGTCGTAGcc                                                                                  >  2:333222/1‑145 (MQ=255)
                                                                gAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGccc                                                                                 >  1:7933/1‑145 (MQ=255)
                                                                gAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGccc                                                                                 <  2:7933/145‑1 (MQ=255)
                                                                 aaTTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCt                                                                                                 <  8:4722/128‑1 (MQ=255)
                                                                 aaTTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCt                                                                                                 >  7:4722/1‑128 (MQ=255)
                                                                  aTTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAATACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGaaa                                                                            >  7:284102/1‑148 (MQ=255)
                                                                  aTTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGaaa                                                                            >  3:18075/1‑148 (MQ=255)
                                                                      aTTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtgt                                                                         <  8:299352/147‑1 (MQ=255)
                                                                      aTTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtgt                                                                         <  4:18075/147‑1 (MQ=255)
                                                                       tttATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTgg                                                                       <  7:164564/148‑1 (MQ=255)
                                                                       tttATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTgg                                                                       >  2:281621/1‑148 (MQ=255)
                                                                       tttATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTgg                                                                       <  1:281621/148‑1 (MQ=255)
                                                                            aCGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGGTACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCtat                                                                   >  5:60988/1‑147 (MQ=255)
                                                                            aCGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGGTACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCtat                                                                   >  2:135105/1‑147 (MQ=255)
                                                                             cGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCtata                                                                  <  5:392596/147‑1 (MQ=255)
                                                                             cGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCtata                                                                  <  5:1074/147‑1 (MQ=255)
                                                                                  gAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTa                                                                                                >  3:31851/1‑112 (MQ=255)
                                                                                  gAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTa                                                                                                <  4:31851/112‑1 (MQ=255)
                                                                                   aatgTTTAAAAAAAAAAAGATTAAGGACTCAAGCACCTTAATTTACTTAAATTTTGAGATATACTCGGATGAGTTATGCAAAAAAACTCATAACTTAAAACTCTTATTCTATCAATGTCGTCGCCCATAAAtgt                                                                         <  6:12147/131‑1 (MQ=255)
                                                                                   aaaGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAATACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTa                                                            <  8:284102/147‑1 (MQ=255)
                                                                                   aaaGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTa                                                            <  7:52580/147‑1 (MQ=255)
                                                                                    aaGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTc                                                                                        <  3:204152/118‑1 (MQ=255)
                                                                                    aaGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACt                                                          >  7:28937/1‑148 (MQ=255)
                                                                                     aGTTTAATAAAAAAAAGATTAAACGATGAAGGAACGACATTTACTTAAATATTCAGAGATACTAGGATGAGGTATGCAAAAAAACTCAAAACAACAAACTCATATTCTATAAAtcaag                                                                                       >  4:204152/1‑114 (MQ=255)
                                                                                         taataaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcat                                                      <  3:225606/147‑1 (MQ=255)
                                                                                              aaaaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGGTACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttca                                                <  6:60988/148‑1 (MQ=255)
                                                                                              aaaaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGGTACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttca                                                <  1:135105/148‑1 (MQ=255)
                                                                                                         aaGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCAtt                                      >  3:384441/1‑147 (MQ=255)
                                                                                                          aGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCAttt                                     <  2:375921/147‑1 (MQ=255)
                                                                                                              aTGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGa                                <  8:28937/148‑1 (MQ=255)
                                                                                                                gAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttc                                                        >  8:15989/1‑122 (MQ=255)
                                                                                                                gAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttc                                                        <  7:15989/122‑1 (MQ=255)
                                                                                                                gAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGat                               >  6:93696/1‑147 (MQ=255)
                                                                                                                   ggAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACATCAAACACTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATTTGCATATAAAGattat                            >  8:223656/1‑147 (MQ=255)
                                                                                                                    gAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGAttatta                          <  5:298579/148‑1 (MQ=255)
                                                                                                                     aaCGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATtattat                         <  4:26280/148‑1 (MQ=255)
                                                                                                                      aCGTCATTTACTGATATATTGATAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATTTCGTAGCCCATAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTATAGATTATTATg                        <  7:223656/148‑1 (MQ=255)
                                                                                                                      aCGTCATTTAATGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATg                        <  1:225128/148‑1 (MQ=255)
                                                                                                                         tCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAATATTATTATGgtg                     >  6:327909/1‑148 (MQ=255)
                                                                                                                          cATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATaa                                                                 >  6:58666/1‑103 (MQ=255)
                                                                                                                          cATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATaa                                                                 <  5:58666/103‑1 (MQ=255)
                                                                                                                           aTTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTc                   <  2:299969/148‑1 (MQ=255)
                                                                                                                           aTTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTc                   <  5:93696/148‑1 (MQ=255)
                                                                                                                             ttACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcata                                              <  5:135427/119‑1 (MQ=255)
                                                                                                                             ttACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcata                                              >  6:135427/1‑119 (MQ=255)
                                                                                                                                  gAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACg                <  7:441888/144‑1 (MQ=255)
                                                                                                                                  gAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACg                >  8:441888/1‑144 (MQ=255)
                                                                                                                                  gAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCAt            >  5:127656/1‑148 (MQ=255)
                                                                                                                                   aaaTATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTAc                                                           <  8:440827/100‑1 (MQ=255)
                                                                                                                                   aaaTATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTAc                                                           >  7:440827/1‑100 (MQ=255)
                                                                                                                                            gagaTACTAGGATGAGGTATGCAAAAAAACTCATATCTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATAGGTGCATTTAAAGATTATTATGGGGTGCCGACATCTCAATAttt  >  2:68659/1‑148 (MQ=255)
                                                                                                                                            gagaTACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTAttt  <  6:127656/148‑1 (MQ=255)
                                                                                                                                            gagaTACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTAttt  <  1:86815/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CTGTAAAATTACTGGAATAATCAGTTTTAATATTGAGCGGCAATGGCATTTAAAAGATATTGCGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAAAGATTATTATGGTGTCACGCCATCTCATTATTT  >  NC_000913/584094‑584381

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