Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,971,926 G→A 10.0% H256Y (CAT→TAT)  tar ← methyl‑accepting chemotaxis protein II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,971,9260GA10.0% 302.9 / 20.3 100H256Y (CAT→TAT) tarmethyl‑accepting chemotaxis protein II
Reads supporting (aligned to +/- strand):  ref base G (42/48);  new base A (4/6);  total (46/54)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

Only 100 of 110 total aligned reads displayed.

CGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGC  >  NC_000913/1971779‑1972060
                                                                                                                                                   |                                                                                                                                      
cggcAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATAtgt                                                                                                                                       <  4:573286/149‑1 (MQ=255)
cggcAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATAtgt                                                                                                                                       <  4:1176819/149‑1 (MQ=255)
      tttCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACg                                                                                                                                 >  1:1263801/1‑149 (MQ=255)
       ttcttcCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGc                                                                                                                                >  1:1371983/1‑149 (MQ=255)
         cttcCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTCACATATGCAACGctc                                                                                                                              <  2:1040612/149‑1 (MQ=255)
         cttcCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGctc                                                                                                                              <  2:1263801/149‑1 (MQ=255)
         cttcCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGctc                                                                                                                              <  4:557340/149‑1 (MQ=255)
          ttcCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTCACATATGCAACGctct                                                                                                                             <  2:1282314/149‑1 (MQ=255)
          ttcCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTCACATATGCAACGctct                                                                                                                             >  2:966713/1‑149 (MQ=255)
          ttcCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGctct                                                                                                                             >  1:338247/1‑149 (MQ=255)
             cAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGaaa                                                                                                                                   >  3:1100514/1‑140 (MQ=255)
             cAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGaaa                                                                                                                                   <  4:1100514/140‑1 (MQ=255)
             cAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTgcg                                                                                                                          >  2:434200/1‑149 (MQ=255)
              aGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTgcgc                                                                                                                         >  4:1253752/1‑149 (MQ=255)
               gcgcGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGcc                                                                                                                        >  3:624222/1‑149 (MQ=255)
                   ggATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGt                                                                                                                    >  1:1172764/1‑149 (MQ=255)
                        ccTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCGTATGTGAAACGCTCTGCGCCAGGTCGccc                                                                                                               <  4:1164562/149‑1 (MQ=255)
                          tgctgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTgc                                                                                                                           >  3:952440/1‑135 (MQ=255)
                          tgctgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTgc                                                                                                                           >  3:1251061/1‑135 (MQ=255)
                          tgctgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTgc                                                                                                                           <  4:952440/135‑1 (MQ=255)
                          tgctgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCAt                                                                                                             >  1:8523/1‑149 (MQ=255)
                          tgctgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCAt                                                                                                             >  4:480117/1‑149 (MQ=255)
                            ctgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCAttt                                                                                                           <  2:974770/149‑1 (MQ=255)
                            ctgTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCAttt                                                                                                           <  2:1383859/149‑1 (MQ=255)
                              gTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCa                                                                                                         <  3:480117/149‑1 (MQ=255)
                               ttCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCAc                                                                                                        >  2:1151467/1‑149 (MQ=255)
                                tCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACt                                                                                                       <  1:1070722/149‑1 (MQ=255)
                                    tACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCg                                                                                                                  >  1:1050517/1‑134 (MQ=255)
                                    tACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCg                                                                                                                  <  2:1050517/134‑1 (MQ=255)
                                      cGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGcct                                                                                                 <  4:639708/149‑2 (MQ=255)
                                            gaAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAAt                                                                                           >  1:7535/1‑149 (MQ=255)
                                              aaGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATgg                                                                                         <  3:1253752/149‑1 (MQ=255)
                                               aGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGt                                                                                        >  3:536896/1‑149 (MQ=255)
                                                 aTCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATAtgt                                                                                                                                       <  4:1256080/100‑1 (MQ=255)
                                                 aTCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATAtgt                                                                                                                                       >  3:1256080/1‑100 (MQ=255)
                                                 aTCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCa                                                                                      <  2:1172764/149‑1 (MQ=255)
                                                  tCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAg                                                                                     <  2:1371983/149‑1 (MQ=255)
                                                  tCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAg                                                                                     >  4:467612/1‑149 (MQ=255)
                                                     gtgtTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGt                                                                                  <  3:1461537/149‑1 (MQ=255)
                                                      tgtTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTa                                                                                 <  1:1151467/149‑1 (MQ=255)
                                                          gCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTTCACATATGCAACGctcttt                                                                                                                           <  4:1203629/103‑3 (MQ=38)
                                                           cccGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGc                                                                                                                 >  4:460381/1‑112 (MQ=255)
                                                            ccgccgCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTgcgc                                                                                                   <  4:143727/125‑1 (MQ=255)
                                                            ccgccgCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTgcgc                                                                                                   >  3:143727/1‑125 (MQ=255)
                                                            ccgccgCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGcc                                                                           >  3:215251/1‑149 (MQ=255)
                                                              gccgcAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAg                                                                         >  3:1026786/1‑149 (MQ=255)
                                                                             gTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATAtgt                                                                                                                                       >  1:756978/1‑72 (MQ=255)
                                                                             gTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATAtgt                                                                                                                                       <  2:756978/72‑1 (MQ=255)
                                                                             gTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGAt                                                          >  2:1077190/1‑149 (MQ=255)
                                                                                   gCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTaccac                                                                 <  4:1428891/136‑1 (MQ=255)
                                                                                   gCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTaccac                                                                 >  3:1428891/1‑136 (MQ=255)
                                                                                   gCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTacca                                                                  <  4:100544/135‑1 (MQ=255)
                                                                                   gCATAGATGGCATCAGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTacca                                                                  >  3:100544/1‑135 (MQ=255)
                                                                                     aTAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGaa                                                  <  4:1026786/149‑1 (MQ=255)
                                                                                      tAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAAt                                                                                           >  2:684961/1‑107 (MQ=255)
                                                                                      tAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAAt                                                                                           <  1:684961/107‑1 (MQ=255)
                                                                                      tAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAAt                                                 <  2:338247/149‑1 (MQ=255)
                                                                                         aTGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAAtgtg                                              <  3:467612/149‑1 (MQ=255)
                                                                                         aTGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAAtgtg                                              <  4:215251/149‑1 (MQ=255)
                                                                                         aTGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAAtgtg                                              <  2:8523/149‑1 (MQ=255)
                                                                                          tGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCa                                                                                                         >  2:632543/1‑89 (MQ=255)
                                                                                          tGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCa                                                                                                         <  1:632543/89‑1 (MQ=255)
                                                                                          tGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGa                                             <  4:624222/149‑1 (MQ=255)
                                                                                             cATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCa                                          <  4:536896/149‑1 (MQ=255)
                                                                                              aTCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCaa                                         <  1:794192/149‑1 (MQ=255)
                                                                                                      cTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAAt                                                                                           >  2:352517/1‑91 (MQ=255)
                                                                                                      cTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAAt                                                                                           <  1:352517/91‑1 (MQ=255)
                                                                                                      cTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAAtttt                                  >  1:740336/1‑148 (MQ=255)
                                                                                                      cTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAAtttt                                  <  2:740336/148‑1 (MQ=255)
                                                                                                          gcgGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGTGCCAGGTCGCCCATTTCACTGCGCCCGTCAAt                                                                                           <  4:931715/87‑1 (MQ=255)
                                                                                                          gcgGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAAt                                                                                           >  3:931715/1‑87 (MQ=255)
                                                                                                               cATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCg                          <  4:1055273/147‑1 (MQ=255)
                                                                                                               cATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCg                          >  3:1055273/1‑147 (MQ=255)
                                                                                                               cATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGa                        <  2:7535/149‑1 (MQ=255)
                                                                                                                  gAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTaccacc                                                                <  2:954449/106‑1 (MQ=255)
                                                                                                                  gAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTaccacc                                                                >  1:954449/1‑106 (MQ=255)
                                                                                                                  gAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGc                                           <  4:1184975/127‑1 (MQ=255)
                                                                                                                  gAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGc                                           >  3:1184975/1‑127 (MQ=255)
                                                                                                                                 tcaAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGa                                             >  3:790988/1‑110 (MQ=255)
                                                                                                                                 tcaAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGa                                             <  4:790988/110‑1 (MQ=255)
                                                                                                                                 tcaAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAg                       <  4:1051082/132‑1 (MQ=255)
                                                                                                                                 tcaAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAg                       >  3:1051082/1‑132 (MQ=255)
                                                                                                                                   aaaGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAAtttt                                  >  3:638747/1‑119 (MQ=255)
                                                                                                                                   aaaGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAAtttt                                  <  4:638747/119‑1 (MQ=255)
                                                                                                                                   aaaGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGt                      <  4:1419186/131‑1 (MQ=255)
                                                                                                                                   aaaGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGt                      >  3:1419186/1‑131 (MQ=255)
                                                                                                                                      gagCGTTTCACATATGCAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGa                        >  1:59161/1‑126 (MQ=255)
                                                                                                                                      gagCGTTTCACATATGCAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGa                        <  2:59161/126‑1 (MQ=255)
                                                                                                                                      gagCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAAttttt                                 <  4:1246756/117‑1 (MQ=255)
                                                                                                                                      gagCGTTGCATATGTGAAACGCTCTGCGCCAGGGCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAAttttt                                 >  3:1246756/1‑117 (MQ=255)
                                                                                                                                       agCGTTTCACATATGCAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGc                                                             >  1:150255/1‑88 (MQ=38)
                                                                                                                                           ttGCATATGTGGAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGaa                                                  >  4:133849/1‑95 (MQ=255)
                                                                                                                                           ttGCATATGTGGAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGaa                                                  <  3:133849/95‑1 (MQ=255)
                                                                                                                                           ttGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAg                                                                         <  1:589940/72‑1 (MQ=255)
                                                                                                                                           ttGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAg                                                                         >  2:589940/1‑72 (MQ=255)
                                                                                                                                                tatGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGa                        >  3:499947/1‑116 (MQ=255)
                                                                                                                                                  tgtgAAACGCTCTGCGCCAGGTCTCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGa                        <  1:1254832/114‑1 (MQ=255)
                                                                                                                                                  tgtgAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGa                        <  4:903678/114‑1 (MQ=255)
                                                                                                                                                  tgtgAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGa                        >  3:903678/1‑114 (MQ=255)
                                                                                                                                                  tgtgAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGAATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGa                        >  2:1254832/1‑114 (MQ=255)
                                                                                                                                                   |                                                                                                                                      
CGGCAGTTTCTTCCAGCGCGGATGCCTGCTGTTCAGTACGGGAGGAAAGATCGGTGTTGCCCGCCGCAATTTCACGGGTACCGGCATAGATGGCATCTGAACCTTCGCGGACATGAGTGACGGTGTCAGTCAAAGAGCGTTGCATATGTGAAACGCTCTGCGCCAGGTCGCCCATTTCACTGCGCCCGTCAATGGTCAGGGTATTCGCCAGGTTACCACCGGCGATTTCGCGAATGTGAGCAATAATTTTTGCCAGCGGAGTAAGCAACATACGGCGAATGC  >  NC_000913/1971779‑1972060

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