Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,033,811 G→C 7.5% R143P (CGG→CCG)  PP_0891 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,033,8110GC 7.5% 291.6 / 13.0 94R143P (CGG→CCG) PP_0891hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (47/40);  new base C (3/4);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

CCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCAC  >  NC_002947/1033664‑1033955
                                                                                                                                                    |                                                                                                                                                  
ccACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCgcg‑gc                                                                                                                                                 >  2:1876541/1‑150 (MQ=255)
 cacaGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCgcg‑gcg                                                                                                                                                >  2:206861/1‑150 (MQ=255)
 cacaGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCACCTGGA‑AGCgcc‑gcg                                                                                                                                                >  1:1203605/1‑150 (MQ=255)
  acaGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑gcgc                                                                                                                                               >  1:250393/1‑150 (MQ=255)
  acaGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑gcgc                                                                                                                                               >  2:1471297/1‑150 (MQ=255)
        aaGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAg                                                                                                                                         <  1:1876541/150‑1 (MQ=255)
         aGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCAACTTGCATGGCCAGTCGCACAGCTACAATTTCCTGTTGCTGGCCACCCGCCTGGA‑AGCGCC‑GCGCATCCAgg                                                                                                                                        <  1:2124062/150‑1 (MQ=255)
          ggCACAGATCGACCTGCAGCAGATCAGTGCCCTTGCCCTCGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACTCACCTGGA‑AGCGCC‑GCGCATCCAggc                                                                                                                                       <  2:1203605/150‑2 (MQ=255)
             acaGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAggt‑ggt                                                                                                                                   >  1:1176721/1‑150 (MQ=255)
             acaGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAggt                                                                                                                                       <  2:892460/147‑1 (MQ=255)
             acaGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAggt                                                                                                                                       >  1:892460/1‑147 (MQ=255)
              caGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGTCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑AGCGCC‑GCGCATCCAggcgggt                                                                                                                                   <  2:1943575/150‑6 (MQ=255)
               aGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGa                                                                                                                                 >  2:2388392/1‑150 (MQ=255)
                gATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGAc                                                                                                                                <  2:1417259/150‑1 (MQ=255)
                gATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGAc                                                                                                                                <  2:1446558/150‑1 (MQ=255)
                    gACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGG‑CAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTg                                                                                                                           <  2:185181/150‑1 (MQ=255)
                        tGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGc                                                                                                                        >  1:49169/1‑150 (MQ=255)
                        tGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGc                                                                                                                        >  2:67925/1‑150 (MQ=255)
                           aGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCctc                                                                                                                     <  2:2141330/150‑1 (MQ=255)
                           aGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCctc                                                                                                                     >  2:285328/1‑150 (MQ=255)
                            gCCTATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCctct                                                                                                                    >  2:1603562/1‑150 (MQ=255)
                              cGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTgg                                                                                                                  <  1:1275807/150‑1 (MQ=255)
                                     tGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATATCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAggt                                                                                                                                       <  2:1810479/123‑1 (MQ=255)
                                     tGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATATCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAggt                                                                                                                                       >  1:1810479/1‑123 (MQ=255)
                                            gCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGaa                                                                                                    <  1:1603562/150‑1 (MQ=255)
                                                 cGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGAttt                                                                                               <  1:2372426/150‑1 (MQ=255)
                                                   gTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTg                                                                                             >  2:2213438/1‑150 (MQ=255)
                                                   gTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTg                                                                                             <  2:2398109/150‑1 (MQ=255)
                                                    tGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTgg                                                                                            >  2:1617862/1‑150 (MQ=255)
                                                     gCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGc                                                                                           <  1:285328/150‑1 (MQ=255)
                                                       acaGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCtt                                                                                         >  2:2004150/1‑150 (MQ=255)
                                                       acaGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑AGCGCC‑GCGCATCCAggcgg                                                                                                                                     >  2:1777603/1‑104 (MQ=255)
                                                       acaGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑AGCGCC‑GCGCATCCAggc‑gg                                                                                                                                    <  1:1777603/107‑1 (MQ=255)
                                                         aGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTtgc                                                                                       <  1:1227307/150‑1 (MQ=255)
                                                           aTCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCtt                                                                                                                            >  2:1509782/1‑111 (MQ=255)
                                                           aTCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCtt                                                                                                                            <  1:1509782/111‑1 (MQ=255)
                                                            tagAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTtgctgc                                                                                    <  1:2072513/148‑1 (MQ=255)
                                                                 ttGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCgggg                                                                               >  2:519331/1‑150 (MQ=255)
                                                                   gcgcCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCt                                                                             <  2:1120171/150‑1 (MQ=255)
                                                                        tGCCGCTGCCGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGGGAACCG‑‑GGGAGTGGGTGTTCCCTTCCAGGT‑GTTGACCATTGAGCCTCTGGTGCTGCCCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGaaa                                                                        <  1:2388392/150‑1 (MQ=255)
                                                                          ccGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGTGTCTAAGGGTCGGCGCGCCGAATATTTTGGCTTGCTGCCGGGGCTGGAAAGc                                                                      <  1:1617862/150‑1 (MQ=255)
                                                                           cGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGcc                                                                     >  1:8520/1‑150 (MQ=255)
                                                                            gCTGACGGTCACCCTGCATGGCCAGTCGCAGAGCGACAATGTCCTGTTGCTGGCCACCCGCCTGGC‑TGGGCG‑GCACTTCTAGGT‑GGTGAACCTTGCGCTTCTGGTGCTACGAGCCGAAGATTTTGGGTTGATTGCGCGGCTGGCATGCCt                                                                    >  2:1141218/1‑150 (MQ=255)
                                                                            gCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGaaa                                                                        >  2:296276/1‑146 (MQ=255)
                                                                            gCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGaaa                                                                        <  1:296276/146‑1 (MQ=255)
                                                                            gCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCt                                                                    <  2:713373/150‑1 (MQ=255)
                                                                              tGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTgc                                                                  <  1:1162777/150‑1 (MQ=255)
                                                                                  ggTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCaa                                                              <  1:67925/150‑1 (MQ=255)
                                                                                  ggTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCaa                                                              <  1:1562919/150‑1 (MQ=255)
                                                                                    tCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAGGCCTGCGCAAGt                                                            >  2:417401/1‑150 (MQ=255)
                                                                                      aCCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTc                                                          >  2:345606/1‑150 (MQ=255)
                                                                                      aCCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTc                                                          >  2:2056908/1‑150 (MQ=255)
                                                                                       cccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGGTGCGGCCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCg                                                         >  2:1244769/1‑150 (MQ=255)
                                                                                        ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGc                                                        >  2:777814/1‑150 (MQ=255)
                                                                                        ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGc                                                        >  2:1854301/1‑150 (MQ=255)
                                                                                        ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGc                                                        >  2:1826561/1‑150 (MQ=255)
                                                                                        ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGc                                                        <  1:2177248/150‑1 (MQ=255)
                                                                                         cTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGcc                                                       >  1:922297/1‑150 (MQ=255)
                                                                                            cATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAgg                                                    >  2:1974833/1‑150 (MQ=255)
                                                                                             aTGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGt                                                   <  1:1826561/150‑1 (MQ=255)
                                                                                                 ccAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGTGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGaa                                               >  1:160211/1‑150 (MQ=255)
                                                                                                  cAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAg                                              <  1:1832983/150‑1 (MQ=255)
                                                                                                   aGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGt                                             >  1:358560/1‑150 (MQ=255)
                                                                                                       gCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCAt                                         <  2:1176721/150‑1 (MQ=255)
                                                                                                         acaGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCa                                       >  1:1148348/1‑150 (MQ=255)
                                                                                                           aGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAAc                                     <  2:677993/150‑1 (MQ=255)
                                                                                                            gCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAAcc                                    <  1:206861/150‑1 (MQ=255)
                                                                                                             cTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCa                                                      <  1:1363061/131‑1 (MQ=255)
                                                                                                             cTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCa                                                      >  2:1363061/1‑131 (MQ=255)
                                                                                                             cTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAAccc                                   <  2:160211/150‑1 (MQ=255)
                                                                                                              tACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCg                                  >  2:2118378/1‑150 (MQ=255)
                                                                                                              tACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCg                                  >  1:1527196/1‑150 (MQ=255)
                                                                                                               aCAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGt                                 <  1:777814/150‑1 (MQ=255)
                                                                                                               aCAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGt                                 <  2:49169/150‑1 (MQ=255)
                                                                                                                cAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTc                                >  1:690291/1‑150 (MQ=255)
                                                                                                                   tGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGACTTTAAATCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGt                             >  1:2090924/1‑150 (MQ=255)
                                                                                                                   gacattGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCAt                                         <  1:803015/133‑1 (MQ=255)
                                                                                                                   gacattGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCAt                                         >  2:803015/6‑138 (MQ=255)
                                                                                                                     tCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAAc                                     <  2:467940/140‑1 (MQ=255)
                                                                                                                     tCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAAc                                     >  1:467940/1‑140 (MQ=255)
                                                                                                                       cTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCCTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCg                         >  2:226752/1‑150 (MQ=255)
                                                                                                                                 ccACCCGCCTGGA‑TGCGCG‑GCGCCTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCgg               <  1:226752/150‑1 (MQ=255)
                                                                                                                                  cACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTTCGGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCggg              >  1:159480/1‑149 (MQ=255)
                                                                                                                                  cACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGt              >  1:2384137/1‑150 (MQ=255)
                                                                                                                                  cACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCAGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGcc                                                                     >  1:733818/1‑95 (MQ=255)
                                                                                                                                  cACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCAGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGcc                                                                     <  2:733818/95‑1 (MQ=255)
                                                                                                                                  cACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTCGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTTCATCAACCCGTCGGTACCGGTTAGCGCGGt              <  2:1409629/150‑1 (MQ=255)
                                                                                                                                   aCCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTg             >  1:253327/1‑150 (MQ=255)
                                                                                                                                     ccGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTg                                                                                             <  1:2096261/68‑1 (MQ=255)
                                                                                                                                     ccGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTg                                                                                             >  2:2096261/1‑68 (MQ=255)
                                                                                                                                     ccGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCt           >  2:1481331/1‑150 (MQ=255)
                                                                                                                                      caccTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTg          >  1:1709166/3‑150 (MQ=255)
                                                                                                                                      cGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTg          <  1:2004150/150‑1 (MQ=255)
                                                                                                                                         cTGGA‑AGCGCC‑GCGCATCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATc       >  2:1219275/1‑150 (MQ=255)
                                                                                                                                          tGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGATAAAAGTAATGATAGTTGTTCTCTGTGTCTTGCTGCCGGTGCTGGAAAGCCTGCGCAAGTTCGTCAGTTTGAAGTCCATCAACCCGTCGGTGCCGGTTAGCGCGGTGCTGATCt      >  1:1225886/1‑150 (MQ=255)
                                                                                                                                             a‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCa   >  2:68324/1‑150 (MQ=255)
                                                                                                                                                gcgcC‑GCGCATCCAGGCGGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCAc  >  1:1257203/1‑150 (MQ=255)
                                                                                                                                                  gcg‑gcgCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTgc                                                                  <  2:854892/83‑1 (MQ=255)
                                                                                                                                                  gcg‑gcgCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTgc                                                                  >  1:854892/1‑83 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGA‑TGCGCG‑GCGCTTCCAGGT‑GGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCAC  >  NC_002947/1033664‑1033955

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