Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,760,361 T→C 6.9% N209D (AAC→GAC)  PP_2413 ← GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,760,3610TC 6.9% 210.6 / 4.0 72N209D (AAC→GAC) PP_2413GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (30/37);  new base C (3/2);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

CCATCTGGTCGATCAGCGCATTGCCCTCCACCACCGTGCCTACCCGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTGGAAACCAGTGGACGGCGTTCGTGCAGGGCTTCC  >  NC_002947/2760214‑2760509
                                                                                                                                                   |                                                                                                                                                    
ccATCTGGTCGATCAGCGCATTGCCCTCCACCACCGTGCCTACCCGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGc                                                                                                                                                    >  1:854546/1‑150 (MQ=255)
   tCTGGTCGATCAGCGCATTGCCCTCCACCACCGTGCCTACCCGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGtt                                                                                                                                                 >  1:1193105/1‑150 (MQ=255)
   tCTGGTCGATCAGCGCATTGCCCTCCACCACCGAGCCTACCCGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGtt                                                                                                                                                 >  2:641549/1‑150 (MQ=255)
       gTCGATCAGCGCATTGCCCTCCACCACCGTGCCTACCCGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACg                                                                                                                                             >  1:1409208/1‑150 (MQ=255)
          gATCAGCGCATTGCCCTCCACCACCGTGCCTACCCGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAg                                                                                                                                          <  2:1193105/150‑1 (MQ=255)
           aTCAGCGCATTGCCCTCCACCACCGTGCCTACCCGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAgg                                                                                                                                         <  2:677099/150‑1 (MQ=255)
               gcgcATTGCCCTCCACCACCGTGCCTACCCGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACa                                                                                                                                     >  1:906772/1‑150 (MQ=255)
                  cATTGCCCTCCACCACCGTGCCTACCCGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGc                                                                                                                                  <  1:120025/150‑1 (MQ=255)
                        ccTCCACCACCGTGCCTACCCGCTCGCTATGAGGGTGGTACTGTGGTAGGGCTTATCTTTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGccac                                                                                                                            <  2:37062/150‑1 (MQ=255)
                        ccTCCACCACCGTGCCTACCCGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGccac                                                                                                                            >  1:224262/1‑150 (MQ=255)
                        ccTCCACCACCGTGCCTACCCGCCCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGccac                                                                                                                            >  1:110229/1‑150 (MQ=255)
                            caccacCGTGCCTACCCGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACa                                                                                                                                     <  1:1422304/137‑1 (MQ=255)
                            caccacCGTGCCTACCCGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACa                                                                                                                                     >  2:1422304/1‑137 (MQ=255)
                               caccGTGCCTACCCGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATa                                                                                                                     >  2:562680/1‑150 (MQ=255)
                               caccGTGCCTACCCGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACCCAGGTACAGGCTGCCACCCACATa                                                                                                                     >  2:318218/1‑150 (MQ=255)
                                           ccGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGc                                                                                                                               >  1:238188/1‑128 (MQ=255)
                                           ccGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGc                                                                                                                               <  2:238188/128‑1 (MQ=255)
                                           ccGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCg                                                                                                         >  2:90859/1‑150 (MQ=255)
                                           ccGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCg                                                                                                         >  2:1056980/1‑150 (MQ=255)
                                           ccGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTCGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCg                                                                                                         >  2:631597/1‑150 (MQ=255)
                                                   tATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTcgc                                                                                                 >  2:784202/1‑150 (MQ=255)
                                                    aTGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTcgcg                                                                                                <  2:1409208/150‑1 (MQ=255)
                                                      gAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGAt                                                                                              >  2:338908/1‑150 (MQ=255)
                                                               aCAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAACATgg                                                                                     <  2:595710/150‑1 (MQ=255)
                                                                  gcagGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGtt                                                                                  >  2:317044/1‑150 (MQ=255)
                                                                  gcagGCGTCGGTTGCGGTCGACCACATAAAGGACCGAACCGTCTGGGGATAAATGCTCGCCCAGCAGGCCGTTGCGGACCTCATAGATATGCACCTGCAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGtt                                                                                  >  1:1033538/1‑150 (MQ=255)
                                                                        gTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACa                                                                            <  1:78070/150‑1 (MQ=255)
                                                                                 ggTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAg                                                                                                               <  1:1365168/106‑1 (MQ=255)
                                                                                 ggTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAg                                                                                                               >  2:1365168/1‑106 (MQ=255)
                                                                                    cGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGAccacca                                                                <  1:306078/150‑1 (MQ=255)
                                                                                     gACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTCGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGac                                                                     <  1:244764/144‑1 (MQ=255)
                                                                                     gACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTCGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGac                                                                     >  2:244764/1‑144 (MQ=255)
                                                                                      aCCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAgg                                                              <  1:90859/150‑1 (MQ=255)
                                                                                      aCCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAgg                                                              <  2:1033538/150‑1 (MQ=255)
                                                                                       ccACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGt                                                             >  2:849661/1‑150 (MQ=255)
                                                                                           aTACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTCGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTc                                                         <  1:631597/150‑1 (MQ=255)
                                                                                            tACAGGTACGAACCGTCCTTGAAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCg                                                        <  1:784202/150‑1 (MQ=255)
                                                                                               aGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTCCAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGgcc                                                     <  2:826006/150‑1 (MQ=255)
                                                                                                    cGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGa                                                <  2:110229/150‑1 (MQ=255)
                                                                                                     gAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGGTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGAc                                               >  2:910406/1‑150 (MQ=255)
                                                                                                      aaCCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACa                                              <  1:427861/150‑1 (MQ=255)
                                                                                                       aCCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAg                                             >  1:466159/1‑150 (MQ=255)
                                                                                                         cGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACGGCGCGACCACCAGGTTGTCGGCCGCCGACAGGt                                           >  1:617028/1‑150 (MQ=255)
                                                                                                         cGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGt                                           <  2:854546/150‑1 (MQ=255)
                                                                                                          gTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGAACACCAGGTTGTCGGCCGCCGACAGGTa                                          <  1:338908/150‑1 (MQ=255)
                                                                                                           tCCTTGTAGAAATTCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAg                                         <  2:670687/150‑1 (MQ=255)
                                                                                                           tCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAg                                         >  2:210811/1‑150 (MQ=255)
                                                                                                             cTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGc                                                                                                                               <  2:1034300/62‑1 (MQ=255)
                                                                                                             cTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGCTACAGGCTGc                                                                                                                               >  1:1034300/1‑62 (MQ=255)
                                                                                                              ttGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGgccgc                                                   <  2:481157/137‑1 (MQ=255)
                                                                                                              ttGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGgccgc                                                   >  1:481157/1‑137 (MQ=255)
                                                                                                                gTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGt                                    <  1:641549/150‑1 (MQ=255)
                                                                                                                 tAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTg                                   >  1:446959/1‑150 (MQ=255)
                                                                                                                 tAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATAAACCAAATAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTg                                   <  1:318218/150‑1 (MQ=255)
                                                                                                                      aTGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTGGAAAc                              <  1:1056980/150‑1 (MQ=255)
                                                                                                                      aTGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTGGAAAc                              <  2:466159/150‑1 (MQ=255)
                                                                                                                        gCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTGGAAACCa                            <  1:1198387/150‑1 (MQ=255)
                                                                                                                        gCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTGGAAACCa                            <  1:849661/150‑1 (MQ=255)
                                                                                                                         cTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATgg                                                                                     <  2:1272753/92‑1 (MQ=255)
                                                                                                                         cTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATgg                                                                                     >  1:1272753/1‑92 (MQ=255)
                                                                                                                          tCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTGGAAACCAGt                          <  2:446959/150‑1 (MQ=255)
                                                                                                                          tCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTGGAAACCAGt                          <  1:910406/150‑1 (MQ=255)
                                                                                                                             cccAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTGGAAACCAGTGGa                       <  2:357427/150‑1 (MQ=255)
                                                                                                                                     gCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTGGAAACCAGTGGACGGCGTTc               >  2:681198/1‑150 (MQ=255)
                                                                                                                                     gCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTGGAAACCAGTGGACGGCGTTc               <  2:906772/150‑1 (MQ=255)
                                                                                                                                           gAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTg                                   <  2:60414/124‑1 (MQ=255)
                                                                                                                                           gAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTg                                   >  1:60414/1‑124 (MQ=255)
                                                                                                                                           gAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTGGAAACCAGTGGACGGCGTTCGTGCAg         >  1:746780/1‑150 (MQ=255)
                                                                                                                                           gAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTGGAAACCAGTGGACGGCGTTCGTGCAg         <  1:210811/150‑1 (MQ=255)
                                                                                                                                            aGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTGGAAACCAGTGGACGGCGTTCGTGCAgg        <  1:317044/150‑1 (MQ=255)
                                                                                                                                              gATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTGGAAACCAGTGGACGGCGTTCGTGCa          >  1:1369333/1‑146 (MQ=255)
                                                                                                                                              gATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTGGAAACCAGTGGACGGCGTTCGTGCa          <  2:1369333/146‑1 (MQ=255)
                                                                                                                                                  ttGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTGGAAACCAGTGGACGGCGTTCGTGCAGGGCTTcc  >  1:762763/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCATCTGGTCGATCAGCGCATTGCCCTCCACCACCGTGCCTACCCGCTCGCTATGAGGGTGGTACAGCAGGCGTCGGTTGCGGTCGACCACATACAGGTACGAACCGTCCTTGTAGAAATGCTCGCCCAGCAGGCTGTTGAGGATGTTGCGTTCACGCAGGTACAGGCTGCCACCCACATACCCCAGGTAGCGGCCGTCGCGATCGAAAATGGGTTGCGACAGCGCGACCACCAGGTTGTCGGCCGCCGACAGGTAGGGCGTGGAAACCAGTGGACGGCGTTCGTGCAGGGCTTCC  >  NC_002947/2760214‑2760509

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: