Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,042,403 A→T 9.0% intergenic (+110/+86) ygfZ → / ← yqfA iron‑sulfur cluster repair protein, plumbagin resistance/hemolysin III family HylIII inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,042,4030AT 9.0% 407.0 / 25.6 122intergenic (+110/+86)ygfZ/yqfAiron‑sulfur cluster repair protein, plumbagin resistance/hemolysin III family HylIII inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (62/49);  new base T (4/7);  total (66/56)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.88e-01

Only 100 of 125 total aligned reads displayed.

TGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAA  >  NC_000913/3042256‑3042551
                                                                                                                                                   |                                                                                                                                                    
tGCATATCGAGCCGCTGCCGTATTCGCTCGAGGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATa                                                                                                                                                      <  3:924225/148‑1 (MQ=255)
tGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATaa                                                                                                                                                     <  1:434218/149‑1 (MQ=255)
 gCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAg                                                                                                                                                    >  1:421234/1‑149 (MQ=255)
   atatCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGc                                                                                                                                                  <  2:261657/149‑1 (MQ=255)
    tatCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCg                                                                                                                                                 >  2:985518/1‑149 (MQ=255)
       cGAGCCGCTGCCGTATTCGCTCGGAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGtt                                                                                                                                               >  3:789123/1‑148 (MQ=255)
       cGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGtt                                                                                                                                               <  3:29985/148‑1 (MQ=255)
         aGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCa                                                                                                                                             >  3:895886/1‑148 (MQ=255)
         aGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGACTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCa                                                                                                                                             >  3:740927/1‑148 (MQ=255)
           ccGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACg                                                                                                                                           <  3:63466/148‑1 (MQ=255)
           ccGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACg                                                                                                                                           <  3:107531/148‑1 (MQ=255)
           ccGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACg                                                                                                                                           <  4:138126/148‑1 (MQ=255)
           ccGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACg                                                                                                                                           >  3:475128/1‑148 (MQ=255)
           ccGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGACATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACg                                                                                                                                           <  4:740927/148‑1 (MQ=255)
              cTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCg                                                                                                                                        >  4:807173/1‑148 (MQ=255)
                gCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGccttat                                                                                                                                     <  2:983290/149‑5 (MQ=39)
                gCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGcctta                                                                                                                                      <  3:11854/148‑4 (MQ=39)
                gCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGcctta                                                                                                                                      <  3:557205/148‑4 (MQ=39)
                 ccGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATc                                                                                                                                    >  2:778214/1‑149 (MQ=255)
                   gTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCg                                                                                                                                  >  1:421956/1‑149 (MQ=255)
                    tATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCgg                                                                                                                                 <  1:706614/149‑1 (MQ=255)
                         gCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATcc                                                                                                                                   >  2:529423/1‑142 (MQ=255)
                         gCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATcc                                                                                                                                   <  1:529423/142‑1 (MQ=255)
                         gCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGt                                                                                                                             >  4:176303/1‑148 (MQ=255)
                          cTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCgat                                                                                                                                >  3:695359/1‑142 (MQ=255)
                            cGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTtgt                                                                                                                          <  4:475954/148‑1 (MQ=255)
                            cGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTtgt                                                                                                                          >  4:793923/1‑148 (MQ=255)
                             gAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGc                                                                                                                        >  1:866400/1‑149 (MQ=255)
                              aaGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGc                                                                                                                        >  4:1001253/1‑148 (MQ=255)
                                  gTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCg                                                                                                                    <  3:391027/148‑1 (MQ=255)
                                   tAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGAt                                                                                                                  <  1:367772/149‑1 (MQ=255)
                                     aaTCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGc                                                                                                                >  2:460576/1‑149 (MQ=255)
                                      aTCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGc                                                                                                                >  4:182574/1‑148 (MQ=255)
                                       tctcTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGcc                                                                                                               <  3:763133/148‑1 (MQ=255)
                                       tctcTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCt                                                                                                              >  2:220125/1‑149 (MQ=255)
                                         tctTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGa                                                                                                            <  1:617796/149‑1 (MQ=255)
                                          ctTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGa                                                                                                            >  4:208844/1‑148 (MQ=255)
                                          ctTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGa                                                                                                            >  3:533787/1‑148 (MQ=255)
                                            ttATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGccttat                                                                                                                                     >  4:705651/1‑117 (MQ=38)
                                            ttATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGccttat                                                                                                                                     <  3:705651/121‑5 (MQ=38)
                                               tCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGa                                                                                                       <  3:465685/148‑1 (MQ=255)
                                                      aGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGccttatc                                                                                                                                    >  2:528141/1‑107 (MQ=38)
                                                      aGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGgcg                                                                                               >  2:687709/1‑149 (MQ=255)
                                                       ggCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGcctta                                                                                                                                      <  3:807856/109‑4 (MQ=38)
                                                       ggCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGcctta                                                                                                                                      <  3:748218/109‑4 (MQ=38)
                                                       ggCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGcctta                                                                                                                                      >  4:748218/1‑106 (MQ=38)
                                                       ggCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGcctta                                                                                                                                      >  4:807856/1‑106 (MQ=38)
                                                       ggCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGgcgc                                                                                              >  2:263359/1‑149 (MQ=255)
                                                        gCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGcgcg                                                                                             >  1:370988/1‑149 (MQ=255)
                                                          aTTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTc                                                                                           >  1:972766/1‑149 (MQ=255)
                                                                 gCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCa                                                                                     <  3:913312/148‑1 (MQ=255)
                                                                   cTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAgg                                                                                   >  3:410040/1‑148 (MQ=255)
                                                                     gATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCt                                                                                >  1:817845/1‑149 (MQ=255)
                                                                          gACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACaaa                                                                           <  1:460576/149‑1 (MQ=255)
                                                                          gACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACaa                                                                            >  4:1036669/1‑148 (MQ=255)
                                                                           aCGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAg                                                                          <  1:595999/149‑1 (MQ=255)
                                                                           aCGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAg                                                                          <  1:220125/149‑1 (MQ=255)
                                                                            cGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAg                                                                          <  4:815126/148‑1 (MQ=255)
                                                                              cTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGa                                                                                                                   <  3:273961/105‑1 (MQ=255)
                                                                                ggCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATa                                                                     <  1:677936/149‑1 (MQ=255)
                                                                                  cgcgTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCATACAAAGGCATAc                                                                    <  3:74984/148‑1 (MQ=255)
                                                                                   gcgTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATAccc                                                                  <  1:687709/149‑1 (MQ=255)
                                                                                    cgTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATAccc                                                                  <  3:316525/148‑1 (MQ=255)
                                                                                       cTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCAtt                                                               >  3:434923/1‑148 (MQ=255)
                                                                                        ttATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTAc                                                             <  1:250067/149‑1 (MQ=255)
                                                                                         tATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTAc                                                             >  4:28796/1‑148 (MQ=255)
                                                                                            cATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGcc                                                          >  3:321157/1‑148 (MQ=255)
                                                                                            cATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGcc                                                          >  3:393799/1‑148 (MQ=255)
                                                                                              tGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTg                                                        <  3:208844/148‑1 (MQ=255)
                                                                                                  tGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCaa                                                   >  1:563050/1‑149 (MQ=255)
                                                                                                   gggATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCaa                                                   >  4:1036470/1‑148 (MQ=255)
                                                                                                   gggATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTCCCCAat                                                  <  2:207699/149‑1 (MQ=255)
                                                                                                    ggATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAata                                                 <  1:263359/149‑1 (MQ=255)
                                                                                                              ccTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGaa                                                                                                                                      <  2:984489/54‑1 (MQ=255)
                                                                                                              ccTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGaa                                                                                                                                      >  1:984489/1‑54 (MQ=255)
                                                                                                               cTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGAt                                      >  1:792816/1‑149 (MQ=255)
                                                                                                                tttATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGAt                                      <  4:533787/148‑1 (MQ=255)
                                                                                                                tttATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATc                                     <  2:792816/149‑1 (MQ=255)
                                                                                                                tttATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATc                                     >  1:332790/1‑149 (MQ=255)
                                                                                                                 ttATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATc                                     >  4:586031/1‑148 (MQ=255)
                                                                                                                    tCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAg                                        >  1:767968/1‑142 (MQ=255)
                                                                                                                    tCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCa                                 >  2:179582/1‑149 (MQ=255)
                                                                                                                    tCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGACAGCTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAg                                        <  2:767968/142‑1 (MQ=255)
                                                                                                                     cacaACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCa                                 >  4:361950/1‑148 (MQ=255)
                                                                                                                       caacaGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGaa                              >  2:257023/1‑149 (MQ=255)
                                                                                                                        aacaGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCTGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGaaa                             <  2:421234/149‑1 (MQ=255)
                                                                                                                        aacaGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACaa                                            >  4:298278/1‑134 (MQ=255)
                                                                                                                        aacaGGGACGAGGCCCGGAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGaaa                             >  1:737599/1‑149 (MQ=255)
                                                                                                                            gggACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGc                                   >  4:456120/1‑139 (MQ=255)
                                                                                                                            gggACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTgg                         <  1:968892/149‑1 (MQ=255)
                                                                                                                             ggACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTgg                         <  4:644641/148‑1 (MQ=255)
                                                                                                                             ggACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTgg                         >  3:298502/1‑148 (MQ=255)
                                                                                                                                   ggCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCacacac                  >  2:927870/1‑149 (MQ=255)
                                                                                                                                            aTCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGa         <  2:332790/149‑1 (MQ=255)
                                                                                                                                             tCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAg        >  2:413557/1‑149 (MQ=255)
                                                                                                                                                gATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCa      >  4:207319/1‑148 (MQ=255)
                                                                                                                                                 aTAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCAc     >  4:75318/1‑148 (MQ=255)
                                                                                                                                                 aTAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCAc     >  3:1011519/1‑148 (MQ=255)
                                                                                                                                                  tAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGcac                      >  1:185258/1‑130 (MQ=255)
                                                                                                                                                   aaGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGaa  >  1:34615/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGCATATCGAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAA  >  NC_000913/3042256‑3042551

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