Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,001,196 A→T 8.0% Q57L (CAG→CTG)  PP_0863 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,001,1960AT 8.0% 254.9 / 12.3 87Q57L (CAG→CTG) PP_0863hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (39/41);  new base T (5/2);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-01

AGGAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGCCCGTGCGG  >  NC_002947/1001049‑1001345
                                                                                                                                                   |                                                                                                                                                       
aGGAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCcagc                                                                                                                                                       >  2:1286577/1‑150 (MQ=255)
 ggAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCcagca                                                                                                                                                      <  1:1920414/150‑1 (MQ=255)
  gAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCcagcag                                                                                                                                                     >  1:1133994/1‑150 (MQ=255)
  gAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCcagcag                                                                                                                                                     >  1:630557/1‑150 (MQ=255)
    aGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGa                                                                                                                                                   <  2:1157663/150‑1 (MQ=255)
        gtggATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGCTTCTGGCGGCCGAGGCGGAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGcc                                                                                                                                               >  1:230561/1‑150 (MQ=255)
        gtggATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGcc                                                                                                                                               >  1:1505299/1‑150 (MQ=255)
          ggATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCTGCTGGATGCCa                                                                                                                                              >  2:1629753/1‑149 (MQ=255)
          ggATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCTGCTGGATGCCa                                                                                                                                              <  1:1629753/149‑1 (MQ=255)
          ggATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCAc                                                                                                                                             >  2:1687242/1‑150 (MQ=255)
          ggATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCAc                                                                                                                                             >  2:1302137/1‑150 (MQ=255)
                 gCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGcc                                                                                                                                      >  2:590971/1‑150 (MQ=255)
                 gCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGcc                                                                                                                                      >  1:669760/1‑150 (MQ=255)
                   cGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCg                                                                                                                                    >  2:1727317/1‑150 (MQ=255)
                       ttGCAGGCCATGCTCGAGGAAACGCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGt                                                                                                                                <  2:55359/150‑1 (MQ=255)
                        tGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGtt                                                                                                                               <  2:1133994/150‑1 (MQ=255)
                         gCAGGCCATGCTCGAGGAAACTCCCTGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAACCGCTGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTg                                                                                                                              <  2:1315084/150‑1 (MQ=255)
                         gCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGTGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTg                                                                                                                              <  1:419461/150‑1 (MQ=255)
                         gCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATACAGAAGAATCGCATGGTTGCCCGGCGTTg                                                                                                                              <  2:1113466/150‑1 (MQ=255)
                          cAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTgg                                                                                                                             <  1:1878627/150‑1 (MQ=255)
                           aGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAAGTTTTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTggg                                                                                                                            >  2:1485286/1‑149 (MQ=255)
                            ggCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGtt                                                                                                                           <  1:599694/150‑1 (MQ=255)
                               cATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCGGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCTGCTGGATGCCACg                                                                                                                                            >  2:1173289/1‑130 (MQ=255)
                               cATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCGGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCTGCTGGATGCCACg                                                                                                                                            <  1:1173289/130‑1 (MQ=255)
                                 tGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTa                                                                                                                      <  1:57771/150‑1 (MQ=255)
                                      gAGGAAACTCCCGGCAAGGCCGGCCAGGCGATTCTGGCGGCAGCGGGGCAGGGCGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCCTCCAGCAGGATGGCACGGTTGCCCGGCGTTGGTTCGGTATCGCt                                                                                                                 >  1:542580/1‑150 (MQ=255)
                                      gAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCt                                                                                                                 <  2:1043998/150‑1 (MQ=255)
                                      gAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCt                                                                                                                 >  1:297664/1‑150 (MQ=255)
                                         gAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGca                                                                                                              >  1:233822/1‑150 (MQ=255)
                                         gAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGca                                                                                                              >  1:892825/1‑150 (MQ=255)
                                         gAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGca                                                                                                              >  2:1471150/1‑150 (MQ=255)
                                           aaCTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGcaca                                                                                                            >  1:1303635/1‑150 (MQ=255)
                                              tCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAggg                                                                                                         <  2:630557/150‑1 (MQ=255)
                                                ccGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTg                                                                                                       >  1:1570084/1‑150 (MQ=255)
                                                ccGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTg                                                                                                       >  1:853025/1‑150 (MQ=255)
                                                   gCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCa                                                                                                    <  2:1570084/150‑1 (MQ=255)
                                                        gccgccCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGcc                                                                                               >  2:725892/1‑150 (MQ=255)
                                                        gccgccCAGGCGATTCTGGCGGCAGCGGCGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGCTCGGTATCGCTGCACAGGGTGGCAGTGcc                                                                                               >  2:1881677/1‑150 (MQ=255)
                                                          cgccCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCAt                                                                                             <  1:590971/150‑1 (MQ=255)
                                                           gccCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATg                                                                                            <  2:806000/150‑1 (MQ=255)
                                                             ccAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑gc                                                                                         <  2:487333/150‑1 (MQ=255)
                                                                 gCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGcac                                                                                     >  2:804807/1‑150 (MQ=255)
                                                                    aTTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCacaac                                                                                  <  1:1054300/150‑1 (MQ=255)
                                                                    aTTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCacaac                                                                                  <  1:181497/150‑1 (MQ=255)
                                                                     ttCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCacaaca                                                                                 >  2:41206/1‑150 (MQ=255)
                                                                      tCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCATTTGCTGCTCGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACAt                                                                                >  2:531017/1‑150 (MQ=255)
                                                                       cTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATg                                                                               <  2:1289366/150‑1 (MQ=255)
                                                                           cggcAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCg                                                                           <  1:909993/150‑1 (MQ=255)
                                                                             gcAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCggg                                                                         <  1:599269/150‑1 (MQ=255)
                                                                              cAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCggtc                                                                        <  2:1505299/150‑3 (MQ=255)
                                                                              cAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGc                                                                        <  2:1612225/150‑1 (MQ=255)
                                                                              cAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGc                                                                        <  2:1608254/150‑1 (MQ=255)
                                                                                 cGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGtt                                                                     <  1:1727317/150‑1 (MQ=255)
                                                                                  ggggCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTg                                                                    >  2:836858/1‑150 (MQ=255)
                                                                                              ggTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATgg                                                       <  1:1285437/150‑1 (MQ=255)
                                                                                               gTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGc                                                      >  1:1475568/1‑150 (MQ=255)
                                                                                               gTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGc                                                      <  2:1475568/150‑1 (MQ=255)
                                                                                                tCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCt                                                     >  2:764233/1‑150 (MQ=255)
                                                                                                  gAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTgg                                                   <  2:892825/150‑1 (MQ=255)
                                                                                                       gcAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTgg                                              >  1:653505/1‑150 (MQ=255)
                                                                                                         aGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCg                                            >  2:1470083/1‑150 (MQ=255)
                                                                                                           ttGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAg                                          >  1:1834540/1‑150 (MQ=255)
                                                                                                           tggctgctGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAg                                          <  2:393570/148‑1 (MQ=255)
                                                                                                             gctgctGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGt                                        <  2:297664/150‑1 (MQ=255)
                                                                                                              ctgctgGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGAATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGa                                           >  2:298703/1‑146 (MQ=255)
                                                                                                              ctgctgGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGAATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGa                                           <  1:298703/146‑1 (MQ=255)
                                                                                                               tgctgGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGa                                      >  1:642240/1‑150 (MQ=255)
                                                                                                               tgctgGGGAAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGGTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGGGGCGAGGTgg                                      >  1:1599707/1‑149 (MQ=255)
                                                                                                                gctgGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGaa                                     <  1:725892/150‑1 (MQ=255)
                                                                                                                   ggggCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCACGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAAttt                                  <  2:698233/150‑1 (MQ=255)
                                                                                                                    gggCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTg                                                                                                 >  2:917693/1‑88 (MQ=255)
                                                                                                                    gggCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTg                                                                                                 <  1:917693/88‑1 (MQ=255)
                                                                                                                    gggCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTg                                 <  2:233822/150‑1 (MQ=255)
                                                                                                                     ggCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTgg                                <  2:230561/150‑1 (MQ=255)
                                                                                                                      gCAGATCCTGCTGGATGGGCGTGGCATCAAGCAGAAGGCCATGGTTGCCCGGCGGTGGTTAGGTATCGCGGAAAACGCCAG‑TGTTTCACGTACGGCCTAGATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGGGAATTTGGt                               >  1:1876587/1‑150 (MQ=255)
                                                                                                                             ctgctgGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAggcg                        >  2:1249217/1‑150 (MQ=255)
                                                                                                                               gctgGATGGGCGTGGCCTCCAGCAGGATGCCAGCGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGACTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAggcggc                      <  2:309062/150‑1 (MQ=255)
                                                                                                                                   gATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGtt                  >  2:339587/1‑150 (MQ=255)
                                                                                                                                    agggtcGTGGCATCCAGCAGGATGCCACGGTTGTCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTc                 <  2:865927/145‑1 (MQ=255)
                                                                                                                                     tGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCa                >  2:180586/1‑150 (MQ=255)
                                                                                                                                       ccacgcccATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACt              <  2:54036/143‑1 (MQ=255)
                                                                                                                                        ccGTGGCATCCTGCTGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTa             >  1:1836122/2‑150 (MQ=255)
                                                                                                                                           tGGCATCCTGCTGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGc          >  1:396457/1‑150 (MQ=255)
                                                                                                                                           tGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCa                                                                                              >  1:1735325/1‑68 (MQ=255)
                                                                                                                                           tGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCa                                                                                              <  2:1735325/68‑1 (MQ=255)
                                                                                                                                           tGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGc          >  2:709407/1‑150 (MQ=255)
                                                                                                                                           tGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGc          <  1:1506443/150‑1 (MQ=255)
                                                                                                                                             gCATCCTGCTGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACATCATGCTCGGGCGTTGCAAGGAGACACGTCTGAACT‑CCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGccc        >  1:1868474/1‑150 (MQ=255)
                                                                                                                                               aTCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGCCCGt      >  1:1197504/1‑150 (MQ=255)
                                                                                                                                               aTCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGCCCGt      >  2:500051/1‑150 (MQ=255)
                                                                                                                                                 ccAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGCCCGTGc    <  1:1687242/150‑1 (MQ=255)
                                                                                                                                                   agcagGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGGACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGCCCGTGCgg  >  1:1797022/1‑150 (MQ=255)
                                                                                                                                                   agcagGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGAGAATTTGGTGCAGGCGGCCGTTCACTACGCCCGTGCgg  >  1:1855225/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                       
AGGAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGCCCGTGCGG  >  NC_002947/1001049‑1001345

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑