Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,536,587 C→G 8.8% A61G (GCC→GGC)  PP_3125 → Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,536,5870CG 8.8% 209.9 / 10.7 68A61G (GCC→GGC) PP_3125Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (35/27);  new base G (3/3);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-01

TGGGGATGACCCTGGAAATGCTGGCCGAAAAGACTGGCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGT  >  NC_002947/3536440‑3536736
                                                                                                                                                   |                                                                                                                                                     
tGGGGATGACCCTGGAAATGCTGGCCGAAAAGACTGGCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCg                                                                                                                                                     >  2:658935/1‑150 (MQ=255)
        aCCCTGGAAATGCTGGCCGAAAAGACTGGCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAGACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAg                                                                                                                                             >  1:1318195/1‑150 (MQ=255)
         cccTGGAAATGCTGGCCGAAAAGACTGGCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAgg                                                                                                                                            >  2:756616/1‑150 (MQ=255)
          cgtGGAAATGCTGGCCGAAAAGACTGGCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGc                                                                                                                                           >  1:984966/3‑150 (MQ=255)
          ccgggAAATGCTGGCCGAAAAGACTGGCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGc                                                                                                                                           >  2:1066202/4‑150 (MQ=255)
          ccTGGAAATGCTGGCCGAAAAGACTGGCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGc                                                                                                                                           >  2:134008/1‑150 (MQ=255)
           cTGGAAATGCTGGCCGAAAAGACTGGCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAagag                                                                                                                                                >  2:334019/1‑144 (MQ=255)
           cTGGAAATGCTGGCCGAAAAGACTGGCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAagag                                                                                                                                                <  1:334019/144‑1 (MQ=255)
           cTGGAAATGCTGGCCGAAAAGACTGGCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGcc                                                                                                                                          >  2:137315/1‑150 (MQ=255)
              gAAATGCTGGCCGAAAAGACTGGCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCggg                                                                                                                                       <  1:134008/150‑1 (MQ=255)
                aaTGCTGGCCGAAAAGACTGGCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCa                                                                                                                                     <  1:137315/150‑1 (MQ=255)
                  tGCTGGCCGAAAAGACTGGCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCaga                                                                                                                                   <  1:756616/150‑1 (MQ=255)
                      ggCCGAAAAGACTGGCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCg                                                                                                                               <  1:1333747/150‑1 (MQ=255)
                                aCTGGCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTg                                                                                                                     <  2:939437/150‑1 (MQ=255)
                                aCTGGCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTg                                                                                                                     >  1:1321712/1‑150 (MQ=255)
                                    gCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTgc                                                                                                                 >  1:875899/1‑150 (MQ=255)
                                    gCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTgc                                                                                                                 >  2:1272894/1‑150 (MQ=255)
                                    gCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTgc                                                                                                                 >  2:1089806/1‑150 (MQ=255)
                                      aTGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTgcgc                                                                                                               >  2:858641/1‑150 (MQ=255)
                                           cAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTAc                                                                                                          >  1:321398/1‑150 (MQ=255)
                                           cAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTAc                                                                                                          >  1:762426/1‑150 (MQ=255)
                                             aaaGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCg                                                                                                        <  1:1272894/150‑1 (MQ=255)
                                                 cTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAgcg                                                                                                    >  2:1150657/1‑150 (MQ=255)
                                                          cAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGt                                                                                                            <  2:584044/133‑1 (MQ=255)
                                                          cAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGt                                                                                                            >  1:584044/1‑133 (MQ=255)
                                                          cAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCt                                                                                           <  2:921199/150‑1 (MQ=255)
                                                           aaGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGc                                                                                                                              <  2:1018748/114‑1 (MQ=255)
                                                           aaGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGc                                                                                                                              >  1:1018748/1‑114 (MQ=255)
                                                            aGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTgg                                                                                         <  1:1080780/150‑1 (MQ=255)
                                                                cGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTgg                                                                                     <  2:984966/150‑1 (MQ=255)
                                                                       ggCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATggg                                                                              <  1:81052/150‑1 (MQ=255)
                                                                         ccTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCa                                                                            <  2:152877/150‑1 (MQ=255)
                                                                          cTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCaa                                                                           >  2:1145208/1‑150 (MQ=255)
                                                                            cAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAgg                                                                         <  2:972247/150‑1 (MQ=255)
                                                                             aaCACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAggg                                                                        >  2:887597/1‑150 (MQ=255)
                                                                                acGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAgg                                                                                               >  1:960195/1‑124 (MQ=255)
                                                                                acGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAgg                                                                                               <  2:960195/124‑1 (MQ=255)
                                                                                 cGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCgg                                                                    <  1:887597/150‑1 (MQ=255)
                                                                                      tCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTAc                                                               <  2:875899/150‑1 (MQ=255)
                                                                                          tCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCgg                                                           <  2:788271/150‑1 (MQ=255)
                                                                                            gctcggccgcgGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGcc                                                         >  2:708123/8‑150 (MQ=255)
                                                                                                   cgcTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCa                                                  <  1:1089806/150‑1 (MQ=255)
                                                                                                      tgagaCTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTc                                               <  2:1318195/146‑1 (MQ=255)
                                                                                                      taaaaCTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTc                                               >  1:179182/3‑150 (MQ=255)
                                                                                                      tGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGAGCAAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTc                                               >  2:156698/1‑150 (MQ=255)
                                                                                                                             gTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATg                                                                                <  1:737639/94‑1 (MQ=255)
                                                                                                                             gTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATg                                                                                >  2:737639/1‑94 (MQ=255)
                                                                                                                                aaTGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGcc                                                                                                                       >  1:589743/1‑52 (MQ=255)
                                                                                                                                aaTGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGcc                                                                                                                       <  2:589743/52‑1 (MQ=255)
                                                                                                                                  tGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGtt                   >  2:1388583/1‑150 (MQ=255)
                                                                                                                                  tGTCGAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGttt                  >  2:120541/1‑150 (MQ=255)
                                                                                                                                      gAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTg               >  2:759776/1‑150 (MQ=255)
                                                                                                                                      gAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTg               >  2:111940/1‑150 (MQ=255)
                                                                                                                                         ctgCTCTTCGGCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGa              <  2:81757/146‑1 (MQ=255)
                                                                                                                                         ctgCTCTTCGGCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGa              >  1:81757/3‑148 (MQ=255)
                                                                                                                                          tgCTCTTCGGCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCctgct                         <  1:159703/135‑1 (MQ=255)
                                                                                                                                          tgCTCTTCGGCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCctgct                         >  2:159703/2‑136 (MQ=255)
                                                                                                                                          tgCTCTTCGGCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTc           >  1:169165/2‑150 (MQ=255)
                                                                                                                                           gCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCa          <  1:1145208/150‑1 (MQ=255)
                                                                                                                                            ctctTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAg         >  2:1036290/1‑150 (MQ=255)
                                                                                                                                              ctTCGGCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAgcc       <  2:169165/150‑1 (MQ=255)
                                                                                                                                              ctTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAgcc       <  1:1150657/150‑1 (MQ=255)
                                                                                                                                               ttCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCAc                                                    <  2:533478/104‑1 (MQ=255)
                                                                                                                                               ttCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCAc                                                    >  1:533478/1‑104 (MQ=255)
                                                                                                                                               ttCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGcc                      >  2:1361139/1‑134 (MQ=255)
                                                                                                                                               ttCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAgccg      >  2:1306397/1‑150 (MQ=255)
                                                                                                                                               ttCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCGGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGcc                      <  1:1361139/134‑1 (MQ=255)
                                                                                                                                                 cGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAgccgcc    >  1:1308759/1‑150 (MQ=255)
                                                                                                                                                   ccGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGACAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGt  >  1:18568/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGGGGATGACCCTGGAAATGCTGGCCGAAAAGACTGGCATGACCAAAAGCTACCTGTCCAAGGTCGAGCGCGGCCTCAACACGCCCTCGATCGCTGCGGCGCTGAAACTGGCGAGGGCGCTGAACGTCAATGTCGAAGAGCTCTTCGCCGAAGAGCAGGCCGGGCAGAGCCGCTACAGCCTGGTGCGCCGTACCGAGCGCCAGGCCCTGGTTGGCGATGGGCAAGGGCCGGGGTACGCGGCCCTCACCAGTCAGATTGGCCAGCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGT  >  NC_002947/3536440‑3536736

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: