Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,475,002 T→G 12.7% intergenic (‑775/+26) algD ← / ← PP_1289 GDP‑mannose 6‑dehydrogenase/DUF328 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,475,0020TG12.7% 190.8 / 22.0 80intergenic (‑775/+26)algD/PP_1289GDP‑mannose 6‑dehydrogenase/DUF328 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (34/35);  new base G (4/6);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

GTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCGCTGCTGGAT  >  NC_002947/1474863‑1475150
                                                                                                                                           |                                                                                                                                                    
gTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGCCCTCTTCGCGGCTGAAGCCGCGaa                                                                                                                                            <  2:401359/150‑1 (MQ=255)
gTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGaa                                                                                                                                            >  1:1236594/1‑150 (MQ=255)
gTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTGAAGCCGCGaa                                                                                                                                            <  2:279843/150‑1 (MQ=255)
 tATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAag                                                                                                                                           >  2:2361703/1‑150 (MQ=255)
    gAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGCAGAgg                                                                                                                                        >  1:1027238/1‑150 (MQ=255)
     aaCGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAaga                                                                                                                                          >  1:1737404/1‑147 (MQ=255)
     aaCGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAaga                                                                                                                                          <  2:1737404/147‑1 (MQ=255)
     aaCGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAggg                                                                                                                                       <  1:576092/150‑1 (MQ=255)
     aaCGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAggg                                                                                                                                       >  2:1540238/1‑150 (MQ=255)
     aaCGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAggg                                                                                                                                       >  2:691427/1‑150 (MQ=255)
         gTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGGGAAGAGGGGGTg                                                                                                                                   >  1:155992/1‑150 (MQ=255)
                ggCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCCTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCt                                                                                                                            <  2:1613412/150‑1 (MQ=255)
                     ttGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTc                                                                                                                       <  2:449154/150‑1 (MQ=255)
                           aTAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCggg                                                                                                                 <  2:823551/150‑1 (MQ=255)
                               aGGGAGGTGCGAAGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGa                                                                                                             >  1:1454939/1‑150 (MQ=255)
                                 ggAGGTGCCACGGATTTTTGGCGTCGGATTCGTGTCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGCCGGCTCACCTCTTCTCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATc                                                                                                           <  1:1111641/150‑1 (MQ=255)
                                  gAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATAACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGGGTGACAGGCTCACTCCGCGGGGTGAtgg                                                                                                          >  1:1015011/1‑148 (MQ=255)
                                      tGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCa                                                                                                      <  2:392505/150‑1 (MQ=255)
                                       gCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAg                                                                                                     <  2:1236594/150‑1 (MQ=255)
                                        cGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAgg                                                                                                    >  1:862024/1‑150 (MQ=255)
                                         gACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGa                                                                                                   >  2:1799116/1‑150 (MQ=255)
                                         gACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGa                                                                                                   >  1:1266198/1‑150 (MQ=255)
                                          aCGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGaa                                                                                                  <  1:1799116/150‑1 (MQ=255)
                                          aCGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGaa                                                                                                  >  1:834897/1‑150 (MQ=255)
                                           cGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGCGTGATCGCGCAGGAac                                                                                                 <  2:294034/150‑1 (MQ=255)
                                             ggattttGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACCCCGCGGAGTGagcct                                                                                                         >  1:1879957/4‑136 (MQ=255)
                                             ggattttGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACCCCGCGGAGTGagcct                                                                                                         <  2:1879957/137‑5 (MQ=255)
                                             gATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAacac                                                                                               >  1:1758131/1‑150 (MQ=255)
                                                ttttGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGCCCTCTTCGCGGCTGAAGCCGCGaa                                                                                                                                            <  2:1516270/102‑1 (MQ=39)
                                                ttttGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGCCCTCTTCGCGGCTGAAGCCGCGaa                                                                                                                                            >  1:1516270/1‑102 (MQ=39)
                                                   atctgTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATg                                                                                         <  1:880117/146‑1 (MQ=255)
                                                      cGTCGGATTCGTGGCGTGATGGTGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTc                                                                                      <  2:2055931/150‑1 (MQ=255)
                                                      cGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTc                                                                                      <  2:1266198/150‑1 (MQ=255)
                                                        tCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAggg                                                                                                                                       <  2:2209958/99‑1 (MQ=255)
                                                        tCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAggg                                                                                                                                       >  1:2209958/1‑99 (MQ=255)
                                                        tCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATAACAGCGGCTGAACTATTCGCGGATCAAGCCGCGAAGACGGCGTGTCGGGCTCACGCCGCGGGGTGATCGCGCAGGAACACCAAATggagat                                                                                    >  2:880117/1‑146 (MQ=255)
                                                             ttCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCCCTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTg                                                                               <  1:2361703/150‑1 (MQ=255)
                                                               cGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGa                                                                             <  1:1469036/150‑1 (MQ=255)
                                                                        aTGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGc                                                                    <  1:1540238/150‑1 (MQ=255)
                                                                         tGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCgggg                                                                                                                <  1:2238260/105‑1 (MQ=255)
                                                                         tGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCgggg                                                                                                                >  2:2238260/1‑105 (MQ=255)
                                                                         tGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCtt                                                                                >  2:1536211/1‑137 (MQ=255)
                                                                         tGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCtt                                                                                <  1:1536211/137‑1 (MQ=255)
                                                                         tGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCa                                                                   <  1:401616/150‑1 (MQ=255)
                                                                            cGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTg                                                                <  2:834897/150‑1 (MQ=255)
                                                                              gTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGCGCAATCCCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTgta                                                              <  1:2358527/150‑1 (MQ=255)
                                                                                ttGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTg                                                                          >  1:1522104/1‑136 (MQ=255)
                                                                                ttGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTg                                                                          <  2:1522104/136‑1 (MQ=255)
                                                                                     gCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTg                                                                                    >  2:1987175/1‑121 (MQ=255)
                                                                                     gCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTg                                                                                    <  1:1987175/121‑1 (MQ=255)
                                                                                           ccTGTTGCAGGGCAGAGGCGGATCAAGGCGGCTAACCTCTTCGCGGCTTAAGACGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGa                                                 <  2:508080/150‑1 (MQ=255)
                                                                                            cTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCcgc                                                                                                                    <  1:665373/82‑1 (MQ=255)
                                                                                            cTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCcgc                                                                                                                    >  2:665373/1‑82 (MQ=255)
                                                                                             tGTTGCAGGGCAGAGGCGGATAACGGCGACTGCCCTATTCGCGGATTCAGCCGCGGAGAGGGCGTGACAGGCTCCCTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGc                                               >  2:884335/1‑150 (MQ=255)
                                                                                                   aGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGa                                                                             >  2:777040/1‑114 (MQ=255)
                                                                                                   aGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGa                                                                             <  1:777040/114‑1 (MQ=255)
                                                                                                   aGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGaa                                         >  2:1938739/1‑150 (MQ=255)
                                                                                                     ggCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACt                                       <  2:1454939/150‑1 (MQ=255)
                                                                                                      gCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGTCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTAccc                                                     <  2:573787/135‑1 (MQ=255)
                                                                                                      gCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTAccc                                                     >  1:573787/1‑135 (MQ=255)
                                                                                                             cGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGAGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGAACAGTACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGc                               >  1:1954489/1‑150 (MQ=255)
                                                                                                               gATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGc                                     >  1:512011/1‑142 (MQ=255)
                                                                                                               gATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGc                                     <  2:512011/142‑1 (MQ=255)
                                                                                                                     ggcggcTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCgg                       >  2:668014/1‑150 (MQ=255)
                                                                                                                       cggcTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCgggg                     <  1:932704/150‑1 (MQ=255)
                                                                                                                          cTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCgg                                                                     <  2:1380246/99‑1 (MQ=255)
                                                                                                                          cTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCgg                                                                     >  1:1380246/1‑99 (MQ=255)
                                                                                                                          cTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCggg                      >  1:1616048/1‑146 (MQ=255)
                                                                                                                          cTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCggg                      <  2:1616048/146‑1 (MQ=255)
                                                                                                                           tGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGa                                                                             >  1:1863201/1‑90 (MQ=255)
                                                                                                                           tGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGa                                                                             <  2:1863201/90‑1 (MQ=255)
                                                                                                                            gcccTCTTCGCGGCTGAAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCtt                                                                                >  1:887894/3‑86 (MQ=38)
                                                                                                                            gcccTCTTCGCGGCTGAAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCtt                                                                                <  2:887894/84‑1 (MQ=38)
                                                                                                                            gcccTCTTCGCGGCTGAAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATc                                            <  2:1793576/120‑1 (MQ=255)
                                                                                                                            gcccTCTTCGCGGCTGAAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATc                                            >  1:1793576/3‑122 (MQ=255)
                                                                                                                            gACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTg                >  1:1386345/1‑150 (MQ=255)
                                                                                                                                tctTCGCGGCTGAAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCgg                       >  2:476102/1‑139 (MQ=255)
                                                                                                                                tctTCGCGGCTGAAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCgg                       <  1:476102/139‑1 (MQ=255)
                                                                                                                                     gcgGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCgctgc       >  1:1477137/1‑150 (MQ=255)
                                                                                                                                     gcgGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCgctgc       >  2:2159225/1‑150 (MQ=255)
                                                                                                                                          ttCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCGCTGCTGGAt  >  1:564581/1‑150 (MQ=255)
                                                                                                                                          gtCAGCGGCGAAGAGGGCGTGACAGGCTCACTCCCCGGGGTGATCGAGCAGGAACACCAAATGTTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACACCTGAGCATCGAACTGCTTGACCTGTTCCGGGGTCCTGATGCACTACCGGAt  <  2:155992/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
GTATGAACGGTCTTATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTTGAGCTGTTCCGGGGTGCTGATGCGCTGCTGGAT  >  NC_002947/1474863‑1475150

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: