Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,904,595 T→G 8.2% I234S (ATT→AGT)  manZ → mannose‑specific enzyme IID component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,904,5950TG 8.2% 427.6 / 30.4 134I234S (ATT→AGT) manZmannose‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base T (64/59);  new base G (2/9);  total (66/68)
Fisher's exact test for biased strand distribution p-value = 5.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

Only 100 of 138 total aligned reads displayed.

CTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTACGCTTGCGGCCTGCTGGGACTG  >  NC_000913/1904447‑1904743
                                                                                                                                                    |                                                                                                                                                    
 tGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTAtt                                                                                                                                                     >  4:500033/1‑149 (MQ=255)
 tGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTAtt                                                                                                                                                     >  4:382512/1‑149 (MQ=255)
  gCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTc                                                                                                                                                    >  1:889549/1‑149 (MQ=255)
    aaaaaCTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTg                                                                                                                                                  <  3:853069/149‑1 (MQ=255)
     aaaaCTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTgg                                                                                                                                                 >  3:262270/1‑149 (MQ=255)
      aaaCTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGa                                                                                                                                                <  1:253311/149‑1 (MQ=255)
       aaCTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAc                                                                                                                                               <  4:364282/149‑1 (MQ=255)
        aCTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAca                                                                                                                                              <  2:514448/149‑2 (MQ=255)
         cTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcag                                                                                                                                             <  2:319453/149‑3 (MQ=255)
          tGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAg                                                                                                                                            >  3:686617/1‑149 (MQ=255)
          tGACGGAAGGGGCGTCTATACTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTGACTACTGTCCAGAATAGTCTGGAcagt                                                                                                                                            <  4:631319/149‑4 (MQ=255)
              ggAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTaa                                                                                                                                        <  1:237213/149‑1 (MQ=255)
              ggAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcagta                                                                                                                                           <  2:473569/146‑5 (MQ=255)
                   gggCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCa                                                                                                                                   >  3:486389/1‑149 (MQ=255)
                         cTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTg                                                                                                                             <  2:891520/149‑1 (MQ=255)
                          tATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTgg                                                                                                                            <  4:23618/149‑1 (MQ=255)
                             ccTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTAc                                                                                                                         <  2:887815/149‑1 (MQ=255)
                              cTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTAcc                                                                                                                        >  4:448376/1‑149 (MQ=255)
                                cGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCAc                                                                                                                      <  3:815339/149‑1 (MQ=255)
                                  gCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActg                                                                                                                    >  4:903170/1‑149 (MQ=255)
                                    cTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgct                                                                                                                  >  3:632560/1‑149 (MQ=255)
                                     tGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCa                                                                                                                                   >  4:426368/1‑131 (MQ=255)
                                     tGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCa                                                                                                                                   <  3:426368/131‑1 (MQ=255)
                                           tCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAcc                                                                                                           >  2:332839/1‑149 (MQ=255)
                                             aTGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCtt                                                                                                         >  3:384216/1‑149 (MQ=255)
                                                ggggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGc                                                                                                      <  4:632560/149‑1 (MQ=255)
                                                ggggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGc                                                                                                      <  3:142701/149‑1 (MQ=255)
                                                 gggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCt                                                                                                     >  3:630536/1‑149 (MQ=255)
                                                      ttGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTAt                                                                                                >  2:258934/1‑149 (MQ=255)
                                                        ggTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTAtgt                                                                                              <  1:52664/149‑1 (MQ=255)
                                                         gTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTAtgtg                                                                                             >  2:366063/1‑149 (MQ=255)
                                                          ttAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCCGACTATTCTGGACCAGTTAATGCCAGTCCTGGTCCCACTGCTGCTGACCTTTGCTTGTATGTgg                                                                                            >  3:215971/1‑149 (MQ=255)
                                                          ttAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGGCCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTgg                                                                                            <  2:930089/149‑1 (MQ=255)
                                                          ttAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTgg                                                                                            >  1:790468/1‑149 (MQ=255)
                                                           tAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGc                                                                                           <  1:443084/149‑1 (MQ=255)
                                                             aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTa                                                                                         >  2:431630/1‑149 (MQ=255)
                                                               aaGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCa                                                                                                                       >  2:775640/1‑117 (MQ=255)
                                                               aaGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACt                                                                                       >  3:868690/1‑149 (MQ=255)
                                                                         ctGTCAACATCACGCTGGTTGTCTCTCGCACTACTGACCAGACGGGCAAAGAACACATTACTACTGTCCAGACTATTCTCGACCAGTTAATGCCAGGCCGGGTACCACTGCTGCTGACCTTTGCTTGTATGCGGCTACTGCGCaaaaaa                                                                             <  4:571319/148‑1 (MQ=255)
                                                                         aTGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCaaaaaa                                                                             <  2:790468/149‑1 (MQ=255)
                                                                          tGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGGACACGTTACTACTGTCCAGACTATTCTGGAc                                                                                                                                               <  1:141063/82‑1 (MQ=255)
                                                                          tGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGGACACGTTACTACTGTCCAGACTATTCTGGAc                                                                                                                                               >  2:141063/1‑82 (MQ=255)
                                                                           gTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGt                                                                           >  3:14664/1‑149 (MQ=255)
                                                                           gTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGt                                                                           <  4:384216/149‑1 (MQ=255)
                                                                            tCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGATGGGAAAAGAACACGTTGCTACTATCCAGACTATTCTTGTCCAGTTAATGCCAGGCCTGGTACCACTGATGCTGACCTGTGCTTGTATCTGGCTAATGCGCGAAAAAgag                                                                          >  3:571319/1‑147 (MQ=255)
                                                                            tCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCaaaaaa                                                                             >  1:529824/1‑146 (MQ=255)
                                                                            tCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGtt                                                                          <  1:183072/149‑1 (MQ=255)
                                                                            tCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGtt                                                                          >  1:358264/1‑149 (MQ=255)
                                                                              aaCATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTaa                                                                        >  1:303115/1‑149 (MQ=255)
                                                                                   cccGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGc                                                                   <  4:701455/149‑1 (MQ=255)
                                                                                    acGCTGGTTGTCTCTCGCCTTACTGACCAGACGGGCCAAGAACCCGGAACTAATGTCAAGACTATTGTGGACAAGTTAATGCCAAGCCTGGTACCACTGCTGCTGACCTTTGCATGTATGTGGCTACTGCGCAAAAAAATTAACCCGCt                                                                  >  4:395023/2‑149 (MQ=255)
                                                                                        tGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTgg                                                              <  3:377748/149‑1 (MQ=255)
                                                                                          gTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGat                                                            <  2:788124/149‑1 (MQ=255)
                                                                                          gTTGTCTCTCGCATGACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGat                                                            >  3:63953/1‑149 (MQ=255)
                                                                                           ttGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAg                                                                                                                                  >  1:136248/1‑78 (MQ=255)
                                                                                           ttGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGatc                                                           <  4:63953/149‑1 (MQ=255)
                                                                                           ttGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGatc                                                           >  4:742773/1‑149 (MQ=255)
                                                                                            tGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGatca                                                          >  4:621130/1‑149 (MQ=255)
                                                                                                tctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgtg                                                               >  1:418799/1‑140 (MQ=255)
                                                                                                tctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGat                                                            <  2:868575/143‑1 (MQ=255)
                                                                                                tctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGat                                                            >  1:868575/1‑143 (MQ=255)
                                                                                                tctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGt                                                      >  1:576515/1‑149 (MQ=255)
                                                                                                  tcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTg                                                    >  1:822125/1‑149 (MQ=255)
                                                                                                  tcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTg                                                    <  1:814983/149‑1 (MQ=255)
                                                                                                  tcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTg                                                    >  3:361106/1‑149 (MQ=255)
                                                                                                   cGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTgg                                                   <  4:262270/149‑1 (MQ=255)
                                                                                                        tACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGcttct                                              <  4:215971/149‑1 (MQ=255)
                                                                                                          cTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCttcttc                                            <  2:822125/149‑1 (MQ=255)
                                                                                                             aCCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTc                                         >  4:152963/1‑149 (MQ=255)
                                                                                                               cAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCAt                                       >  3:10644/1‑149 (MQ=255)
                                                                                                               cAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCAt                                       >  1:679251/1‑149 (MQ=255)
                                                                                                                aGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATc                                      >  3:593220/1‑149 (MQ=255)
                                                                                                                aGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCCTCTTCGTCATc                                      <  1:693947/149‑1 (MQ=255)
                                                                                                                 gACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgt                                                                <  3:866518/122‑1 (MQ=255)
                                                                                                                 gACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgt                                                                >  4:866518/1‑122 (MQ=255)
                                                                                                                    gggCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTa                                  <  3:692278/149‑1 (MQ=255)
                                                                                                                     ggCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTAt                                 <  3:170235/149‑1 (MQ=255)
                                                                                                                       cAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTa                                                                                                 >  4:265829/1‑83 (MQ=255)
                                                                                                                       cAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTa                                                                                                 <  3:265829/83‑1 (MQ=255)
                                                                                                                           gAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCa                                                                                                                       >  2:546340/1‑57 (MQ=255)
                                                                                                                           gAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTgg                           <  4:552598/149‑1 (MQ=255)
                                                                                                                                 gTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTAc                       <  3:844880/147‑1 (MQ=255)
                                                                                                                                 gTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTACGc                     <  4:473952/149‑1 (MQ=255)
                                                                                                                                 gTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGATATCGCTGGTTAc                       >  4:844880/1‑147 (MQ=255)
                                                                                                                                   tactacTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTACGCtt                   >  1:55161/1‑149 (MQ=255)
                                                                                                                                      tactGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTACGCTTGCg                >  1:692241/1‑149 (MQ=255)
                                                                                                                                      tactGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTACGCTTGCg                >  4:692996/1‑149 (MQ=255)
                                                                                                                                        ctGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGACATCGGTATCGCTGGTTACGCTTGCggt              >  4:594191/1‑148 (MQ=255)
                                                                                                                                           tCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGtt                                                                          <  4:879748/86‑1 (MQ=255)
                                                                                                                                            ccAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTACTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTACGCTTGCGGCctgc          <  3:594191/149‑1 (MQ=255)
                                                                                                                                             cAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTTCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTACGCTTGCGGCctgct         >  4:321022/1‑149 (MQ=255)
                                                                                                                                             cAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTACGCTTGCGGCctgct         >  1:409079/1‑149 (MQ=255)
                                                                                                                                             cAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTACGCTTGCGGCctgct         <  2:508090/149‑1 (MQ=255)
                                                                                                                                              aGACTATTCTGGACCAGTTAATGCCAGGCGTGGTACCACTGCTGCTTACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGCTCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTACGCTTGCGGCctgctg        >  2:572324/1‑149 (MQ=255)
                                                                                                                                               gACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTACGCTTGCGGCCTGCTgg       >  2:471541/1‑149 (MQ=255)
                                                                                                                                               gACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTACGCTTGCGGCCTGCTgg       <  3:382512/149‑1 (MQ=255)
                                                                                                                                                 cTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTACGCTTGCGGCCTGCTGGGa     <  3:870212/149‑1 (MQ=255)
                                                                                                                                                  tATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTACGCTTGCGGCCTGCTGGGAc    >  1:508281/1‑149 (MQ=255)
                                                                                                                                                   aTTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTATATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTACGCTTGCGGCCTGCTGGGACt   >  2:490558/1‑149 (MQ=255)
                                                                                                                                                    ttCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTACGCTTGCGGCCTGCTGGGACTg  >  1:34726/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGCAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCGCTGGTTACGCTTGCGGCCTGCTGGGACTG  >  NC_000913/1904447‑1904743

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