Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,123,345 C→T 5.5% A261A (GCG→GCA PP_5371 ← rubredoxin/rubredoxin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,123,3450CT 5.5% 514.4 / 12.7 164A261A (GCG→GCAPP_5371rubredoxin/rubredoxin reductase
Reads supporting (aligned to +/- strand):  ref base C (83/71);  new base T (4/5);  total (88/76)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 170 total aligned reads displayed.

GCACTGAGCACCAGGTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGC  >  NC_002947/6123196‑6123493
                                                                                                                                                     |                                                                                                                                                    
gcactgAGCACCAGGTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTgcc                                                                                                                                                      <  2:2874157/150‑1 (MQ=255)
gcactgAGCACCAGGTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTgcc                                                                                                                                                      <  1:1707499/150‑1 (MQ=255)
  actgAGCACCAGGTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTgccgc                                                                                                                                                    >  1:728295/1‑150 (MQ=255)
    tgAGCACCAGGTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCt                                                                                                                                                  >  2:3302541/1‑150 (MQ=255)
      aGCACCAGGTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGAAAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCg                                                                                                                                                >  1:4241479/1‑150 (MQ=255)
       gCACCAGGTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCgg                                                                                                                                               <  2:503264/150‑1 (MQ=255)
        cACCAGGTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCggg                                                                                                                                              <  2:3872616/150‑1 (MQ=255)
          ccAGGTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGgca                                                                                                                                            <  2:1586492/150‑1 (MQ=255)
              gTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATCGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGcagcag                                                                                                                                        <  2:3777117/150‑1 (MQ=255)
               tCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAgg                                                                                                                                       <  2:3034430/150‑1 (MQ=255)
                 tGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCTGCCCGAGGCGGCAGGcc                                                                                                                                     <  2:2593784/150‑2 (MQ=255)
                  gCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCTGCCCGAGGCGGCAGGcca                                                                                                                                    <  1:246879/150‑3 (MQ=255)
                  gCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCTGCCCGAGGCGGCAGGcca                                                                                                                                    <  1:2522872/150‑3 (MQ=255)
                  gCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCTGCCCGAGGCGGCAGGcca                                                                                                                                    <  2:796997/150‑3 (MQ=255)
                  gCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGt                                                                                                                                    <  2:2847867/150‑1 (MQ=255)
                     ccGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTcc                                                                                                                                 <  1:4623917/150‑1 (MQ=255)
                     ccGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTcc                                                                                                                                 >  1:3449447/1‑150 (MQ=255)
                      cgacgaCCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCa                                                                                                                                >  1:4888729/1‑150 (MQ=255)
                      cgacgaCCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCa                                                                                                                                >  1:2317681/1‑150 (MQ=255)
                       gacgacCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCgctatcg                                                                                                                               <  1:5067788/150‑7 (MQ=255)
                       gacgacCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCAt                                                                                                                               >  2:5005419/1‑150 (MQ=255)
                            ccTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGcc                                                                                                                          <  1:3400100/150‑1 (MQ=255)
                            ccTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGACCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGcc                                                                                                                          >  2:1907907/1‑150 (MQ=255)
                            ccTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGACCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGcc                                                                                                                          >  2:3623011/1‑150 (MQ=255)
                            ccTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCAATGCACCCCCAACCCCGATAGAGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGAACCAAGGGcc                                                                                                                          >  2:565422/1‑150 (MQ=255)
                              tGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAg                                                                                                                        <  2:3427144/150‑1 (MQ=255)
                                aCCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGc                                                                                                                      >  1:2248122/1‑150 (MQ=255)
                                  cGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCAc                                                                                                                    >  1:3007948/1‑150 (MQ=255)
                                  cGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCAc                                                                                                                    >  1:3083032/1‑150 (MQ=255)
                                   gTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCAcc                                                                                                                   <  1:907328/150‑1 (MQ=255)
                                             gTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGTCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGAtt                                                                                                         >  2:418610/1‑150 (MQ=255)
                                             gTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGAtt                                                                                                         <  1:440431/150‑1 (MQ=255)
                                              tCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTa                                                                                                        >  2:1266847/1‑150 (MQ=255)
                                                 gTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTAccc                                                                                                     >  2:1365144/1‑150 (MQ=255)
                                                     gTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAAc                                                                                                 <  2:2871292/150‑1 (MQ=255)
                                                     gTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAAc                                                                                                 >  1:1286046/1‑150 (MQ=255)
                                                         ccTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGa                                                                                             <  1:1527649/150‑1 (MQ=255)
                                                           ttCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATa                                                                                           >  2:874796/1‑150 (MQ=255)
                                                           ttCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATa                                                                                           >  2:4192554/1‑150 (MQ=255)
                                                           ttCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATa                                                                                           <  2:1610739/150‑1 (MQ=255)
                                                           ttCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATa                                                                                           >  2:2786862/1‑150 (MQ=255)
                                                           ttCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATa                                                                                           >  1:1652123/1‑150 (MQ=255)
                                                           ttCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATa                                                                                           <  1:4192442/150‑1 (MQ=255)
                                                             cGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGc                                                                                         <  1:3755428/150‑1 (MQ=255)
                                                                 gTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCAtt                                                                                     >  1:3255265/1‑150 (MQ=255)
                                                                 gTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGGTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCAtt                                                                                     >  2:3580896/1‑150 (MQ=255)
                                                                         cGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCaa                                                                             <  2:4278015/150‑1 (MQ=255)
                                                                             gCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTc                                                                         <  1:3865041/150‑1 (MQ=255)
                                                                             gCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTc                                                                         <  2:2487347/150‑1 (MQ=255)
                                                                               aCTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGt                                                                       <  1:2979667/150‑1 (MQ=255)
                                                                                cTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGtt                                                                      >  1:484058/1‑150 (MQ=255)
                                                                                cTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGtt                                                                      <  2:3139672/150‑1 (MQ=255)
                                                                                cTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGtt                                                                      <  2:3137683/150‑1 (MQ=255)
                                                                                cTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGtt                                                                      <  1:791387/150‑1 (MQ=255)
                                                                                 tGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGtt                                                                      >  2:353097/1‑149 (MQ=255)
                                                                                 tGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTg                                                                     >  2:4915059/1‑150 (MQ=255)
                                                                                 tGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTg                                                                     >  2:188076/1‑150 (MQ=255)
                                                                                 tGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTg                                                                     <  2:114726/150‑1 (MQ=255)
                                                                                 tGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTg                                                                     >  1:4065464/1‑150 (MQ=255)
                                                                                   cAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGc                                                                   >  1:2025841/1‑150 (MQ=255)
                                                                                   cAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGc                                                                   <  2:4334295/150‑1 (MQ=255)
                                                                                    aaGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTgg                                                                    <  1:2191607/148‑1 (MQ=255)
                                                                                    aaGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTgg                                                                    >  2:2191607/1‑148 (MQ=255)
                                                                                       gTCGTGCGCAGGCTCCATTGCACCCCCACCCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAAc                                                               <  1:4967700/150‑1 (MQ=255)
                                                                                                aGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATCGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCg                                                      >  2:4042486/1‑150 (MQ=255)
                                                                                                aGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCg                                                      <  1:1005680/150‑1 (MQ=255)
                                                                                                aGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCg                                                      >  1:3884147/1‑150 (MQ=255)
                                                                                                aGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCg                                                      >  2:1462935/1‑150 (MQ=255)
                                                                                                   cTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATc                                                   >  1:1398/1‑150 (MQ=255)
                                                                                                    tCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCa                                                  >  2:724244/1‑150 (MQ=255)
                                                                                                    tCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCa                                                  >  2:3634969/1‑150 (MQ=255)
                                                                                                    tCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGAGCa                                                  >  2:4162521/1‑150 (MQ=255)
                                                                                                        ttGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAAcc                                              >  1:2353075/1‑150 (MQ=255)
                                                                                                            aCCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCACCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCg                                          <  2:2680563/150‑1 (MQ=255)
                                                                                                              ccccAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGcc                                        >  1:3373184/1‑150 (MQ=255)
                                                                                                                ccAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCa                                      <  1:514959/150‑1 (MQ=255)
                                                                                                                 cAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAg                                     >  1:3646845/1‑150 (MQ=255)
                                                                                                                    ccccGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAAt                                  >  1:3595498/1‑150 (MQ=255)
                                                                                                                     cccGATAGCGCTTGCTGCTAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGACCAGAATg                                 <  2:4486123/150‑1 (MQ=255)
                                                                                                                     cccGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATg                                 <  2:3780946/150‑1 (MQ=255)
                                                                                                                       cGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGc                               <  2:1979229/150‑1 (MQ=255)
                                                                                                                       cGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGc                               >  2:2482866/1‑150 (MQ=255)
                                                                                                                       cGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGc                               >  1:567383/1‑150 (MQ=255)
                                                                                                                             cgcTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCAGTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCActct                         <  2:5060508/150‑1 (MQ=255)
                                                                                                                                  gctgcAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCa                    >  2:4791452/1‑150 (MQ=255)
                                                                                                                                  gctgcAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCa                    <  2:4888729/150‑1 (MQ=255)
                                                                                                                                   ctgcAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCAc                   <  2:1636170/150‑1 (MQ=255)
                                                                                                                                    tgcAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACg                  >  1:556418/1‑150 (MQ=255)
                                                                                                                                      cAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGcc                >  2:1412250/1‑150 (MQ=255)
                                                                                                                                      cAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGcc                >  1:1319220/1‑150 (MQ=255)
                                                                                                                                       aaGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCa               >  1:4156081/1‑150 (MQ=255)
                                                                                                                                          ggaacgCCTGCTGCCCGAGGCGGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCa            >  1:2785127/6‑150 (MQ=255)
                                                                                                                                             tGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAgg         >  2:3327791/1‑150 (MQ=255)
                                                                                                                                               gCCTGCTGCCCGAGGCGGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCg       >  1:4717980/1‑150 (MQ=255)
                                                                                                                                                 cTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGtt     >  2:2347965/1‑150 (MQ=255)
                                                                                                                                                 cTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGtt     <  1:4682886/150‑1 (MQ=255)
                                                                                                                                                 cTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGtt     >  1:1387015/1‑150 (MQ=255)
                                                                                                                                                 cTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGtt     >  2:740530/1‑150 (MQ=255)
                                                                                                                                                  tGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTg    >  1:504512/1‑150 (MQ=255)
                                                                                                                                                    ccgccTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAATCTCGCAGCCGTTCAATCCGTCGCCCAGAATGGCCACTCTGCGCGCGCCATCAAGGCGTTGGc  >  2:3936852/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCACTGAGCACCAGGTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGC  >  NC_002947/6123196‑6123493

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