Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 591,085 T→C 5.4% V110A (GTG→GCG)  PP_0501 → NAD‑dependent epimerase/dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947591,0850TC 5.4% 323.5 / 3.9 105V110A (GTG→GCG) PP_0501NAD‑dependent epimerase/dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base T (49/50);  new base C (4/2);  total (53/52)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

TGCAGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCT  >  NC_002947/590936‑591234
                                                                                                                                                     |                                                                                                                                                       
tGCAGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGt                                                                                                                                                         <  1:180651/150‑1 (MQ=255)
tGCAGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGt                                                                                                                                                         <  2:233903/150‑1 (MQ=255)
 gCAGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTg                                                                                                                                                        <  2:421413/150‑1 (MQ=255)
  cAGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGc                                                                                                                                                       <  2:1142464/150‑1 (MQ=255)
   aGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCa                                                                                                                                                      <  1:1596571/150‑1 (MQ=255)
   aGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTAGCGGCGGTAGCGTAGGTTCAGGCATCGGTCGAAGACACGGTGCGTACCCACCAGAAGAACTTCATCGGTACCCTCAGCGATTGAGAAGACATGCGCGTGCa                                                                                                                                                      >  2:14429/1‑150 (MQ=255)
     ggCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTACCTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATg                                                                                                                                                    <  2:885154/150‑1 (MQ=255)
        tGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCg                                                                                                                                                 >  1:378597/1‑150 (MQ=255)
         ggTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGt                                                                                                                                                <  1:14429/150‑1 (MQ=255)
           tCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTgc                                                                                                                                              <  1:1139639/150‑1 (MQ=255)
                  ggCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGt                                                                                                                                     <  1:2118719/150‑1 (MQ=255)
                   gCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTg                                                                                                                                    >  1:1381971/1‑150 (MQ=255)
                     tGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGct                                                                                                                                  >  2:1743829/1‑150 (MQ=255)
                      gCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑Cgg                                                                                                                                       >  1:1407149/1‑144 (MQ=255)
                      gCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑Cgg                                                                                                                                       <  2:1407149/144‑1 (MQ=255)
                        ggctggctGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCtt                                                                                                                               >  1:1467170/1‑150 (MQ=255)
                               tgCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTcc                                                                                                                        >  2:1932296/1‑150 (MQ=255)
                               tgCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTcc                                                                                                                        >  2:520232/1‑150 (MQ=255)
                                 ccGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCcag                                                                                                                      >  1:79931/1‑150 (MQ=255)
                                      cTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCcagcgcag                                                                                                                 <  1:1743829/150‑1 (MQ=255)
                                         tggtgCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCgg                                                                                                              >  2:218535/1‑150 (MQ=255)
                                                tCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACgg                                                                                                       >  1:2190258/1‑150 (MQ=255)
                                                  tGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGc                                                                                                                                                       >  2:1799355/1‑102 (MQ=255)
                                                  tGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGc                                                                                                                                                       <  1:1799355/102‑1 (MQ=255)
                                                          gTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTc                                                                                                                         <  1:1372052/122‑1 (MQ=255)
                                                          gTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTc                                                                                                                         >  2:1372052/1‑122 (MQ=255)
                                                             gCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGaa                                                                                          <  1:1928739/150‑1 (MQ=255)
                                                               gTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGc                                                                                                                            >  1:1399469/1‑114 (MQ=255)
                                                               gTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGc                                                                                                                            <  2:1399469/114‑1 (MQ=255)
                                                               gTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAgg                                                                                        <  1:617918/150‑1 (MQ=255)
                                                               gTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAgg                                                                                        <  1:1477916/150‑1 (MQ=255)
                                                                tCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAggg                                                                                       >  1:2066690/1‑150 (MQ=255)
                                                                        ggCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGAt                                                                               >  2:429676/1‑150 (MQ=255)
                                                                        ggCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGAt                                                                               <  2:752424/150‑1 (MQ=255)
                                                                         gCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCCCTGCTACCGCTATCGATGCCATGGCGCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGAtt                                                                              >  1:274444/1‑150 (MQ=255)
                                                                           aTCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTc                                                                            <  2:1808803/150‑1 (MQ=255)
                                                                           aTCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTc                                                                            <  2:114720/150‑1 (MQ=255)
                                                                           aTCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTc                                                                            >  2:2144667/1‑150 (MQ=255)
                                                                               gTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGaa                                                                        >  2:18101/1‑150 (MQ=255)
                                                                                  gAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGcaac                                                                                                   <  1:596594/120‑1 (MQ=255)
                                                                                  gAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGcaac                                                                                                   >  2:596594/1‑120 (MQ=255)
                                                                                  gAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGac                                                                     >  1:1439841/1‑150 (MQ=255)
                                                                                        ccGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCcagcgcagc                                                                                                                <  2:1411611/101‑1 (MQ=255)
                                                                                        ccGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCcagcgcagc                                                                                                                >  1:1411611/1‑101 (MQ=255)
                                                                                        ccGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGc                                                                                                      <  2:528958/111‑1 (MQ=255)
                                                                                        ccGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGc                                                                                                      >  1:528958/1‑111 (MQ=255)
                                                                                         cGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCaa                                                              >  1:470164/1‑150 (MQ=255)
                                                                                          ggTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCaaa                                                             <  1:2213417/150‑1 (MQ=255)
                                                                                             gCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGgc                                                          >  1:361315/1‑150 (MQ=255)
                                                                                                tACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGcc                                                       <  2:378597/150‑1 (MQ=255)
                                                                                                tACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGcc                                                       <  1:2144667/150‑1 (MQ=255)
                                                                                                    caccaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTg                                                   >  1:131173/1‑150 (MQ=255)
                                                                                                     accaGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGa                                                  >  1:211197/1‑150 (MQ=255)
                                                                                                       caGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGGCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGAcc                                                >  1:1890380/1‑150 (MQ=255)
                                                                                                         gagTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGAcccc                                              >  2:1842868/1‑150 (MQ=255)
                                                                                                             aaCTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGCAGACACGCCAAAGGCGCCGTTGACCCCTTAt                                          <  2:1439841/150‑1 (MQ=255)
                                                                                                             aaCTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTAt                                          <  1:1891641/150‑1 (MQ=255)
                                                                                                                ttCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGaa                                                                                          >  2:1301329/1‑99 (MQ=255)
                                                                                                                ttCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGaa                                                                                          <  1:1301329/99‑1 (MQ=255)
                                                                                                                ttCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGcc                                       >  1:1081024/1‑150 (MQ=255)
                                                                                                                 tCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCg                                      <  2:131173/150‑1 (MQ=255)
                                                                                                                  cATCGGTACGCTCAACGTTTGGGAGGCCATGCGGGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGt                                     <  2:13698/150‑1 (MQ=255)
                                                                                                                  cATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCTGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGt                                     <  2:79931/150‑1 (MQ=255)
                                                                                                                  cATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGt                                     <  2:1304142/150‑1 (MQ=255)
                                                                                                                   aTCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCATCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTg                                    <  2:1431187/150‑1 (MQ=255)
                                                                                                                   aTCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAGAGGCGCCGTTGACCCCTTATGCCGTg                                    <  2:1467170/150‑1 (MQ=255)
                                                                                                                   aTCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTg                                    >  2:443559/1‑150 (MQ=255)
                                                                                                                   aTCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTg                                    <  1:1932296/150‑1 (MQ=255)
                                                                                                                       gTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGt                                     <  2:2221495/145‑1 (MQ=255)
                                                                                                                             tCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTgg                          >  1:1836964/1‑150 (MQ=255)
                                                                                                                             tCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTgg                          >  1:1502027/1‑150 (MQ=255)
                                                                                                                               aaCGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTg                           >  2:1222752/1‑147 (MQ=255)
                                                                                                                               aaCGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTg                           <  1:1222752/147‑1 (MQ=255)
                                                                                                                                aCGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCa                       <  2:1381971/150‑1 (MQ=255)
                                                                                                                                     gcgcacgCCATGCACGCGCATGGCTTCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGc                 >  2:900771/7‑150 (MQ=255)
                                                                                                                                       gcacgCCATGCACGCGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCa                >  1:1526191/5‑150 (MQ=255)
                                                                                                                                        gAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGc                             <  2:157297/136‑1 (MQ=255)
                                                                                                                                        gAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGc                             >  1:157297/1‑136 (MQ=255)
                                                                                                                                        cacgCCATGCACGCGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATg                                                                                               >  2:1633342/4‑70 (MQ=38)
                                                                                                                                        cacgCCATGCACGCGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATg                                                                                               <  1:1633342/67‑1 (MQ=38)
                                                                                                                                         acgCCATGCACGCGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGGAACAATGGTGAAGGGGAATCGATTTCCGAAGAAACGCCAAAGGCGCCGTTGACCCCTTATGCCGTg                                    >  1:892963/3‑128 (MQ=255)
                                                                                                                                         acgCCATGCACGCGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTg                                    <  2:892963/126‑1 (MQ=255)
                                                                                                                                         aaGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGt              <  1:1842868/150‑1 (MQ=255)
                                                                                                                                           gCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTAc            <  1:18101/150‑1 (MQ=255)
                                                                                                                                            ccATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTAcc           >  1:1987203/1‑150 (MQ=255)
                                                                                                                                              aTGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATg                                                                                               <  1:2220371/64‑1 (MQ=255)
                                                                                                                                              aTGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATg                                                                                               <  1:2220367/64‑1 (MQ=255)
                                                                                                                                              aTGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATg                                                                                               >  2:2220371/1‑64 (MQ=255)
                                                                                                                                              aTGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATg                                                                                               >  2:2220367/1‑64 (MQ=255)
                                                                                                                                              aTGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGAccc                                               <  2:811265/112‑1 (MQ=255)
                                                                                                                                              aTGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTg         <  2:1836964/150‑1 (MQ=255)
                                                                                                                                              aTGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTg         >  1:2196166/1‑150 (MQ=255)
                                                                                                                                              aTGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTg         >  2:255896/1‑150 (MQ=255)
                                                                                                                                                gcgcGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGa       <  1:218535/150‑1 (MQ=255)
                                                                                                                                                 cgcgTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGAc      <  2:2196166/150‑1 (MQ=255)
                                                                                                                                                  gcgTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTActaccgcg     >  2:2053629/1‑143 (MQ=255)
                                                                                                                                                   cgTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACtt    <  1:443559/150‑1 (MQ=255)
                                                                                                                                                    gtgtATGGCGTGCTC‑‑CGGGTGCTCTTCGCCTCCAGGGCAGAGGTGTATGTGAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACTTGGACTTc   >  2:329524/1‑150 (MQ=255)
                                                                                                                                                     tGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGTTTTTTTCTTCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCt  <  1:520232/150‑1 (MQ=255)
                                                                                                                                                     tGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCt  <  2:470164/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
TGCAGGGCTGGTCACTCAGGCTGCGGCTGGCTGCCGTGCTGTGGTGCATCTGGCGGCGGTAGCGTCGGTTCAGGCATCGGTCGAAGACCCGGTGCGTACCCACCAGAGTAACTTCATCGGTACCCTCAACGTTTGCGAAGCCATGCGCGTGCATGGCGTGCGC‑‑CGGGTGCTCTTCGCCTCCAGCGCAGCGGTGTACGGCAACAATGGTGAAGGGGAATCGATTTCCGAAGACACGCCAAAGGCGCCGTTGACCCCTTATGCCGTGGACAAGCTGGCCAGTGAGCAGTACCTGGACTTCT  >  NC_002947/590936‑591234

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: