Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,355,480:1 +T 11.5% coding (5/264 nt) PP_1180 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,355,4801.T11.5% 207.4 / 24.5 61coding (5/264 nt)PP_1180hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (29/25);  new base T (3/4);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

CCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAGGCCTGGCCTGCAACACCCAT  >  NC_002947/1355334‑1355629
                                                                                                                                                     |                                                                                                                                                      
ccTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAatgcatgc                                                                                                                                                     <  1:81774/150‑1 (MQ=255)
ccTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTGTGATATTGTCCCGCCTCGATATGGTCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGCTACGCAatgcatgc                                                                                                                                                     >  1:461808/1‑150 (MQ=255)
 cTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATgcg                                                                                                                                                    >  1:874981/1‑150 (MQ=255)
  tGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATgcgc                                                                                                                                                   <  2:294313/150‑1 (MQ=255)
   gCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAatgca                                                                                                                                                         <  2:1270721/144‑1 (MQ=255)
   gCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAatgca                                                                                                                                                         >  1:1270721/1‑144 (MQ=255)
    cGGCACTTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCATGCATTGCGCAt                                                                                                                                                 <  2:476110/150‑1 (MQ=255)
    cGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCAt                                                                                                                                                 <  2:1246146/150‑1 (MQ=255)
        aCCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCGCGCGGGCGGATGCGCATGCATTGCGCATcc                                                                                                                                               >  2:52619/1‑150 (MQ=255)
              gCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCGCGCGGGCGGATGCGCATGCATTGCGCATccccgcga                                                                                                                                         <  1:52619/150‑7 (MQ=38)
                cTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCgg                                                                                                                                     >  1:849285/1‑150 (MQ=255)
                 tGCGGGTTCTCGCGCGGCGTGTCGGATTGCATGGCAGGGAGTCTCCATAGTCTCTGTATTGCTCAGGTCCCTTATGGGCGCGCACCGTTCCGTCCACTTGCTCGATAACCTCGCCG‑‑GCATGCGCCCGGCATGCGCACTAGCCCGCGCggg                                                                                                                                    >  2:410280/1‑150 (MQ=255)
                          tcGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTAGCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGAc                                                                                                                           >  1:1157957/1‑150 (MQ=255)
                               aGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCgg                                                                                                                      <  1:82232/150‑1 (MQ=255)
                                 tagtgtCGGGTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCTCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGa                                                                                                                    <  1:410280/148‑1 (MQ=255)
                                 tCGTGTCTGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGa                                                                                                                    <  1:81631/150‑1 (MQ=255)
                                   gtgtCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACa                                                                                                                  >  2:16580/1‑150 (MQ=255)
                                     gtCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGa                                                                                                                <  2:461808/150‑1 (MQ=255)
                                            tGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGa                                                                                                         <  1:988870/150‑1 (MQ=255)
                                              cATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTa                                                                                                                                <  1:504882/125‑1 (MQ=255)
                                              cATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTa                                                                                                                                >  2:504882/1‑125 (MQ=255)
                                                      ggtgTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCg                                                                                                                       >  2:741968/1‑126 (MQ=255)
                                                      ggtgTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCg                                                                                                                       <  1:741968/126‑1 (MQ=255)
                                                       gtgtCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCaa                                                                                              >  2:217268/1‑150 (MQ=255)
                                                          tctcCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GCATGCGCCATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATgg                                                                                           >  2:185451/1‑150 (MQ=255)
                                                            tcCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCt                                                                                         >  2:1365077/1‑150 (MQ=255)
                                                            tcCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCACCTGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCt                                                                                         >  1:86790/1‑150 (MQ=255)
                                                              cacaGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAACGGCTTg                                                                                       >  2:1363205/1‑150 (MQ=255)
                                                                 aGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCCCTCGCGCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCg                                                                                    <  1:797629/150‑1 (MQ=255)
                                                                  gTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGc                                                                                   >  2:16160/1‑150 (MQ=255)
                                                                     tCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCTCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGCAATAAGGGGCAATGGCTTGCCGCTcc                                                                                <  1:403141/150‑1 (MQ=255)
                                                                             gTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGa                                                                        >  2:1320271/1‑150 (MQ=255)
                                                                              ttCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGa                                                                                                                    <  2:1050659/105‑1 (MQ=255)
                                                                              ttCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGa                                                                                                                    >  1:1050659/1‑105 (MQ=255)
                                                                                    cgcCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCGCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCaaa                                                                 <  1:127046/150‑1 (MQ=255)
                                                                                     gcCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAg                                                                <  1:281063/150‑1 (MQ=255)
                                                                                       cTTATGGGCGCCCACCGTTCCGTCCATTTGCTAGATAACCTCGCGG‑‑GGATGCGGAATGCATGCGCATCCGCCCGAGCGGGCGTACAGACTTCGGGACAGACTAAGGAATAAGGGGCAATGGATTTCCGCTCCCTGGCAGAAGTTAAAGGt                                                              >  1:930977/1‑150 (MQ=255)
                                                                                        ttATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGa                                                                                                                    >  2:82472/1‑95 (MQ=255)
                                                                                        ttATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGa                                                                                                                    <  1:82472/95‑1 (MQ=255)
                                                                                        ttATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAg                                                                      >  2:282721/1‑141 (MQ=255)
                                                                                        ttATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAg                                                                      <  1:282721/141‑1 (MQ=255)
                                                                                        ttATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGtt                                                            >  1:1219823/1‑150 (MQ=255)
                                                                                        ttATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGtt                                                            >  1:729989/1‑150 (MQ=255)
                                                                                            gggCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTgg                                                        >  2:1029466/1‑150 (MQ=255)
                                                                                               cgcCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGcc                                                     <  1:1365077/150‑1 (MQ=255)
                                                                                                        ttCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATaa                                                                                                     <  1:562966/94‑1 (MQ=255)
                                                                                                        ttCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATaa                                                                                                     >  2:562966/1‑94 (MQ=255)
                                                                                                            gTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACt                                        >  1:1196836/1‑150 (MQ=255)
                                                                                                                aCTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGCCAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTg                                    >  1:348717/1‑150 (MQ=255)
                                                                                                                 cTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGc                                   <  2:874981/150‑1 (MQ=255)
                                                                                                                  ctGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGCCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCt                                  >  2:974591/2‑150 (MQ=255)
                                                                                                                         aTAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGt                           <  1:217268/150‑1 (MQ=255)
                                                                                                                              cTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAg                      <  2:1219823/150‑1 (MQ=255)
                                                                                                                                  cgggcggATGCGCATGCATTGCGCATC‑‑CCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAAGAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGTCAggcct                  >  2:1024481/6‑150 (MQ=255)
                                                                                                                                   gggcggATGCGCATGCATTGCGCATC‑‑CCCGCGAGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAggcctg                 <  2:1343128/146‑1 (MQ=255)
                                                                                                                                   gg‑‑ggATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAggcctg                 <  1:1320271/150‑1 (MQ=255)
                                                                                                                                       tgATGCGCATGCATTGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTg                                                         <  1:1260496/107‑1 (MQ=255)
                                                                                                                                       tgATGCGCATGCATTGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTg                                                         >  2:1260496/2‑108 (MQ=255)
                                                                                                                                       ggATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAACGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAggcctggc               <  1:1363205/150‑1 (MQ=255)
                                                                                                                                        gATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGAGTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAggcctggcc              >  1:259555/1‑150 (MQ=255)
                                                                                                                                             gcAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGaa                                                                       <  1:87905/89‑1 (MQ=255)
                                                                                                                                             gcAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGaa                                                                       >  2:87905/1‑89 (MQ=255)
                                                                                                                                             gcAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCt                                                          >  1:266521/1‑101 (MQ=255)
                                                                                                                                             gcAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCt                                                          <  2:266521/101‑1 (MQ=255)
                                                                                                                                             gcAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAGGCCTGGCCTGCaa         >  1:147546/1‑150 (MQ=255)
                                                                                                                                                atgcattgCGCATC‑‑CCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGTCAGGCCTGGCCTGCAACAcc     <  1:1024481/150‑1 (MQ=255)
                                                                                                                                                atgcatgcGCATCCGCCAGCGCGGGCAGACCGACGTCGGGGCAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTAACTGGCCGACCTCCAACATCACTGACTCTACCCTGAAACTTCCGCTTATGGGTGCCAGGCCTGGCCTGCAacac      >  2:312118/1‑150 (MQ=255)
                                                                                                                                                atgcatgcACATCAGCACGCGCGGGAGTACCGACATCGGTGCAGGCGCAGGAATAAGGGGCAAGGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCAGAGCCCAAAAACTGTTGCTGATGTGTGCCAGGCCGGGCCTGCAacac      >  2:311677/1‑150 (MQ=255)
                                                                                                                                                   catgcGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAGGCCTGGCCTGCAACACCCa   >  2:864410/1‑150 (MQ=255)
                                                                                                                                                    atgcGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAGGCCTGGCCTGCAACACCCAt  >  2:937637/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGG‑‑GGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGTACCGACTTCGGGACAGACTCAGGAATAAGGGGCAATGGCTTGCCGCTCCCTGGCCGAAGTCAAAGGTTCTGGGCCGAGCCCAAAAACTGTTGCTGATGGGTGCCAGGCCTGGCCTGCAACACCCAT  >  NC_002947/1355334‑1355629

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