Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,001,199 A→T 8.0% Q58L (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,1990AT 8.0% 259.3 / 12.3 87Q58L (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 = 4.12e-01

GGAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGCCCGTGCGG  >  NC_002947/1001050‑1001345
                                                                                                                                                     |                                                                                                                                                    
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                                                                                                                 >  1:297664/1‑150 (MQ=255)
                                     gAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCt                                                                                                                 <  2:1043998/150‑1 (MQ=255)
                                        gAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGca                                                                                                              >  1:892825/1‑150 (MQ=255)
                                        gAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGca                                                                                                              >  1:233822/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:181497/150‑1 (MQ=255)
                                                                   aTTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCacaac                                                                                  <  1:1054300/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                                           <  1:298703/146‑1 (MQ=255)
                                                                                                             ctgctgGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGAATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGa                                           >  2:298703/1‑146 (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                                                                                              <  2:1735325/68‑1 (MQ=255)
                                                                                                                                          tGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCa                                                                                              >  1:1735325/1‑68 (MQ=255)
                                                                                                                                          tGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGc          <  1:1506443/150‑1 (MQ=255)
                                                                                                                                          tGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGc          >  2:709407/1‑150 (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)
                                                                                                                                                     |                                                                                                                                                    
GGAAGTGGTGGATGTGGCCGGCTTGCAGGCCATGCTCGAGGAAACTCCCGGCAAGGCCGCCCAGGCGATTCTGGCGGCAGCGGGGCAGGGTGTGGTCGAGGCGCAGTTGCTGCTGGGGCAGATCCTGCTGGATGGGCGTGGCATCCAGCAGGATGCCACGGTTGCCCGGCGTTGGTTCGGTATCGCTGCACAGGGTGGCAGTGCCATG‑GCGCACAACATGCTCGGGCGTTGCCTGGAG‑CATGGCTGGGGTGGCGAGGTGAATTTGGTGCAGGCGGCCGTTCACTACGCCCGTGCGG  >  NC_002947/1001050‑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: ‑