Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 579,127 T→C 100% intergenic (‑234/+57) borD ← / ← ybcV DLP12 prophage; putative lipoprotein/DLP12 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913579,1270TC98.7% 10270.5 / 31.5 3609intergenic (‑234/+57)borD/ybcVDLP12 prophage; putative lipoprotein/DLP12 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (24/28);  new base C (1743/1792);  total (1774/1835)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.56e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 3875 total aligned reads displayed.

ATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTG  >  NC_000913/578980‑579273
                                                                                                                                                   |                                                                                                                                                  
    ttAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTTCGGGAGTGTCCGGGAATAATTAAAACGATGCACACATTGTTTAGCGCTTACAcgtattg                                                                                                                                              >  4:211953/1‑149 (MQ=255)
    ttAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACAcgtattg                                                                                                                                              >  6:27596/1‑149 (MQ=255)
       atatTTCTTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGACCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGCTGCACACAGGGTTTAGCGCGTACAcgtattgcat                                                                                                                                           <  5:334147/149‑1 (MQ=255)
       atatTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACAcgtattgcat                                                                                                                                           >  4:28077/1‑149 (MQ=255)
        tatttattAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACgtattgcatt                                                                                                                                          <  6:16427/149‑1 (MQ=255)
        tatttattAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACgtattgcatt                                                                                                                                          <  3:216647/149‑1 (MQ=255)
          tttattAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGGGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGAGCGTACACGTATTGcattat                                                                                                                                        >  1:27896/1‑149 (MQ=255)
          tttattAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGCCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGcattat                                                                                                                                        <  2:10942/149‑1 (MQ=255)
           ttattaATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCAATATg                                                                                                                                       >  1:100015/1‑149 (MQ=255)
               taATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCa                                                                                                                                    <  4:25077/148‑1 (MQ=255)
                aaTGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCTGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCaa                                                                                                                                   <  8:136474/148‑1 (MQ=255)
                aaTGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCACTTACGGGGAACTAGTCTGCGTGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCa                                                                                                                                    >  2:302584/1‑147 (MQ=255)
                 atgtatgtCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAacg                                                                                                                                 >  1:162354/1‑148 (MQ=255)
                 atgtatgtCAGGTGCGAGGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGGGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAac                                                                                                                                  <  2:164701/148‑1 (MQ=255)
                       gtcAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCgg                                                                                                                           >  2:140294/3‑149 (MQ=255)
                        tcAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGt                                                                                                                          >  5:52364/2‑149 (MQ=255)
                         cAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTg                                                                                                                         >  1:35174/1‑149 (MQ=255)
                               cGATGAATCGTCATTGTCTTCCCGGATTAACTATGGCCACCGCCCTGCCGGGGAAGTTCTCTGGGGGAGTGTCCGGGAATAATTATAACGATGCACACCGGGTTTAGCGCGTAGACGTATTGCATTTTGCCAACGCCCCGGTGCTGaga                                                                                                                   <  8:301827/149‑3 (MQ=255)
                                 aTGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACg                                                                                                                 <  4:61291/149‑1 (MQ=255)
                                  tGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACg                                                                                                                 >  4:335816/1‑148 (MQ=255)
                                    aaTCGTCATTGTATTCCCGTATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCAGGGGCTGACACGga                                                                                                               <  5:142137/148‑1 (MQ=255)
                                     aTCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACgg                                                                                                                >  6:131611/1‑146 (MQ=255)
                                      tCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGacac                                                                                                                  >  7:55217/1‑143 (MQ=255)
                                        gTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACAcgtatt                                                                                                                                               <  8:194869/112‑1 (MQ=255)
                                        gTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGaaa                                                                                                          >  4:365793/1‑149 (MQ=255)
                                               tATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACg                                                                                                   >  6:70458/1‑149 (MQ=255)
                                               tATTCCCGGATTAACTATGTCCACAGACCTGACGGGGAACTTCACTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTAAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGAGCTGACACGAAAGAAACCGGACg                                                                                                   >  8:356185/1‑149 (MQ=255)
                                                 ttCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGtt                                                                                                 >  7:122354/1‑149 (MQ=255)
                                                   cccGGTTTAACTATGTCCACAGCACTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACTGGGTTTAGCGCGTACTCGTATTGCATTATGCCAACGCGCCGGTGCTGACACGGAAGAAACCGGACGTTa                                                                                                <  4:144770/148‑1 (MQ=255)
                                                   cccGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGATTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTa                                                                                                <  7:182785/148‑1 (MQ=255)
                                                   cccGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTa                                                                                                <  5:78293/148‑1 (MQ=255)
                                                   cccGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTAt                                                                                               >  5:248407/1‑149 (MQ=255)
                                                    ccGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTAt                                                                                               >  6:30315/1‑148 (MQ=255)
                                                    ccGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTAt                                                                                               >  3:83651/1‑148 (MQ=255)
                                                       gactAACTATGTCCACAGCACTGACGGTTAACTTCTCGGCGGGAGGGTCCGTTAATAATTAAAACGATGCACCCAGGGGTTAGCGCTTACACGTATTTCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGAtt                                                                                           <  8:274627/146‑1 (MQ=255)
                                                         ttAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTa                                                                                         <  3:243329/149‑1 (MQ=255)
                                                            aCTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGa                                                                                                                     <  6:345677/118‑1 (MQ=255)
                                                            aCTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCg                                                                                      >  2:85276/1‑149 (MQ=255)
                                                            aCTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCg                                                                                      >  1:217570/1‑149 (MQ=255)
                                                             cTATGTCCACAGCCCTGACGGGGAAATTATCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGCGTTTAGCGCGTACACGTATTGCTTTATGCCCAAGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATAATTTAGc                                                                                       >  7:336866/1‑147 (MQ=255)
                                                              tATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGAc                                                                                                    >  6:222036/1‑133 (MQ=255)
                                                               aTGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGCTGCACACAGGGTGTAGCGCGTACACGTATTGCATTATGCCAACTCCCCGGTGCTGACACGGAAGAACCCGGACGTTATGATTTAGCggg                                                                                    >  7:82648/1‑146 (MQ=255)
                                                                tGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTg                                                                                    <  5:27596/147‑1 (MQ=255)
                                                                   ccACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGa                                                                                  >  6:19225/1‑146 (MQ=255)
                                                                   ccACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGa                                                                                  >  3:130146/1‑146 (MQ=255)
                                                                    cacaGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGa                                                                                  >  5:293778/1‑145 (MQ=255)
                                                                    cacaGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGa                                                                                  >  3:19206/1‑145 (MQ=255)
                                                                     acaGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGAt                                                                             >  4:19479/1‑149 (MQ=255)
                                                                     acaGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGAt                                                                             >  4:310061/1‑149 (MQ=255)
                                                                     acaGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGAt                                                                             >  1:325194/1‑149 (MQ=255)
                                                                       aGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGAttt                                                                           <  4:88767/149‑1 (MQ=255)
                                                                          ccTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTtgtg                                                                        <  4:11230/149‑1 (MQ=255)
                                                                            tGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTtgtg                                                                        >  7:242638/1‑147 (MQ=255)
                                                                              aCGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTa                                                                      <  5:222886/147‑1 (MQ=255)
                                                                              aCGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAgt                                                                    >  8:383143/1‑149 (MQ=255)
                                                                               cGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTATAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAgtg                                                                   <  2:351649/149‑1 (MQ=255)
                                                                                 gggAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGtt                                                                 <  8:256431/149‑1 (MQ=255)
                                                                                  ggAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTc                                                                <  3:365950/149‑1 (MQ=255)
                                                                                   gAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTACTGTTCt                                                               >  2:245228/1‑149 (MQ=255)
                                                                                   gAACTTCTCTGCGGGAGGGTCCTGGAAAAATTAACACGAGGCACACAGGGTTTAGCGCGTAAACGTATTGCATTATGCCATCGCCCCGGTGCTGACACGGAAGAAACGGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCt                                                               <  6:94339/149‑1 (MQ=255)
                                                                                      cTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGaa                                                            <  8:163243/149‑1 (MQ=255)
                                                                                      cTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGaa                                                            >  3:54096/1‑149 (MQ=255)
                                                                                       ttCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAAt                                                           >  7:211228/1‑149 (MQ=255)
                                                                                       ttCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAAt                                                           <  1:8459/149‑1 (MQ=255)
                                                                                        tctcTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATg                                                          >  1:378248/1‑149 (MQ=255)
                                                                                         ctctGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATg                                                          <  7:306910/148‑1 (MQ=255)
                                                                                          tctGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTa                                                                      >  7:151168/1‑135 (MQ=255)
                                                                                            tGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGctct                                                      >  3:167556/1‑149 (MQ=255)
                                                                                            tGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGct                                                        <  8:114000/147‑1 (MQ=255)
                                                                                              cGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGTGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGctct                                                      <  5:194936/147‑1 (MQ=255)
                                                                                              cGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGctct                                                      >  5:360436/1‑147 (MQ=255)
                                                                                              cGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGctc                                                       <  5:173316/146‑1 (MQ=255)
                                                                                                       ccGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCACCGGAGCTGACACGGAAGAAACAGGACGTTATGAt                                                                                            >  2:274527/1‑100 (MQ=255)
                                                                                                         gggAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTaa                                         >  2:377793/1‑149 (MQ=255)
                                                                                                          ggAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAAt                                        <  6:249622/149‑1 (MQ=255)
                                                                                                           gAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAAt                                        >  7:288461/1‑148 (MQ=255)
                                                                                                           gAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAAt                                        >  8:46741/1‑148 (MQ=255)
                                                                                                             ataatTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCTTTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGaa                                     <  6:10524/149‑1 (MQ=255)
                                                                                                                 attaaaaCGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAAtt                                   <  8:224013/145‑1 (MQ=255)
                                                                                                                 attaaaaCGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTa                                  >  8:7147/4‑149 (MQ=255)
                                                                                                                  ttaaaaCGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCAGTTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGAATTGTGTAGTGTTCTGAATTCTCTCAGATAATAGTAATGAATTAt                                 >  1:133696/3‑149 (MQ=255)
                                                                                                                    aaaaCGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGATGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTATAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATc                                >  4:171710/1‑148 (MQ=255)
                                                                                                                     aaaCGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCaa                              <  3:331802/149‑1 (MQ=255)
                                                                                                                       aCGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCTGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAg                            <  1:244472/149‑1 (MQ=255)
                                                                                                                         gATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGt                          <  1:303344/149‑1 (MQ=255)
                                                                                                                           tGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGta                         >  3:299205/1‑148 (MQ=255)
                                                                                                                             cacacaGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGtata                       <  7:268362/148‑1 (MQ=255)
                                                                                                                              acacacGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATa                                             >  6:383997/1‑125 (MQ=255)
                                                                                                                                acaGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAg                      >  2:54715/1‑146 (MQ=255)
                                                                                                                                acaGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAg                      <  1:54715/146‑1 (MQ=255)
                                                                                                                                acaGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTaa                   >  2:323804/1‑149 (MQ=255)
                                                                                                                                 caGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAat                  <  7:387730/149‑1 (MQ=255)
                                                                                                                                  aGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGACCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAata                 >  2:218902/1‑149 (MQ=255)
                                                                                                                                  aGGGTTTAGCGCGTACACGTATTGCATTATGACAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATAAAAGGTATAGTAata                 >  5:113616/1‑149 (MQ=25)
                                                                                                                                    ggTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATg                                                          >  8:327055/1‑106 (MQ=255)
                                                                                                                                    ggTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATc               <  7:359431/149‑1 (MQ=255)
                                                                                                                                      tttAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCtt             >  3:107303/1‑149 (MQ=255)
                                                                                                                                         aGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCtttt           <  2:192665/148‑1 (MQ=255)
                                                                                                                                             cgTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTaaa                                               <  1:188798/108‑1 (MQ=255)
                                                                                                                                              gTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGtt      >  3:275260/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTG  >  NC_000913/578980‑579273

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: