Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,123,355 A→G 5.5% L258P (CTG→CCG)  PP_5371 ← rubredoxin/rubredoxin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,123,3550AG 5.5% 522.8 / 11.8 164L258P (CTG→CCG) PP_5371rubredoxin/rubredoxin reductase
Reads supporting (aligned to +/- strand):  ref base A (87/68);  new base G (4/5);  total (91/73)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

Only 100 of 175 total aligned reads displayed.

CCAGGTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGCGGAACGCCTGG  >  NC_002947/6123206‑6123504
                                                                                                                                                     |                                                                                                                                                     
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)
       tGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCg                                                                                                                                                <  1:2794588/150‑1 (MQ=255)
        gCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCTGCCCGAGGCGGCAGGcca                                                                                                                                               <  2:1903924/150‑3 (MQ=255)
        gCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCTGCCCGAGGCGGCAGGcca                                                                                                                                               <  1:2522872/150‑3 (MQ=255)
        gCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCTGCCCGAGGCGGCAGGcca                                                                                                                                               <  1:246879/150‑3 (MQ=255)
        gCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGt                                                                                                                                               <  2:2848045/150‑1 (MQ=255)
           ccGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTcc                                                                                                                                            <  1:4623917/150‑1 (MQ=255)
           ccGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTcc                                                                                                                                            <  2:129539/150‑1 (MQ=255)
            cgacgaCCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCa                                                                                                                                           >  1:2317681/1‑150 (MQ=255)
            cgacgaCCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCa                                                                                                                                           >  1:4888729/1‑150 (MQ=255)
             gacgacCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCgctatcg                                                                                                                                          <  1:5067736/150‑7 (MQ=255)
             gacgacCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCAt                                                                                                                                          >  2:5005419/1‑150 (MQ=255)
                gacCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATggg                                                                                                                                       >  2:1527649/1‑150 (MQ=255)
                  ccTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGcc                                                                                                                                     <  1:3400100/150‑1 (MQ=255)
                     gACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAgg                                                                                                                                  <  2:4435989/150‑1 (MQ=255)
                      aCCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGc                                                                                                                                 >  1:2248110/1‑150 (MQ=255)
                        cGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCAc                                                                                                                               >  1:3007948/1‑150 (MQ=255)
                        cGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCAc                                                                                                                               >  1:3083032/1‑150 (MQ=255)
                                 gggTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGa                                                                                                                      <  1:4130537/150‑1 (MQ=255)
                                   gTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGTCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGAtt                                                                                                                    >  2:418610/1‑150 (MQ=255)
                                       gTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTAccc                                                                                                                >  2:1365144/1‑150 (MQ=255)
                                           gTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAAc                                                                                                            >  1:1286046/1‑150 (MQ=255)
                                               ccTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGa                                                                                                        <  1:1527649/150‑1 (MQ=255)
                                                 ttCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATa                                                                                                      >  2:4192554/1‑150 (MQ=255)
                                                 ttCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATa                                                                                                      <  1:4192082/150‑1 (MQ=255)
                                                 ttCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATa                                                                                                      >  2:874796/1‑150 (MQ=255)
                                                 ttCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATa                                                                                                      >  1:1652123/1‑150 (MQ=255)
                                                  tCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAg                                                                                                     >  1:551047/1‑150 (MQ=255)
                                                       gTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCAtt                                                                                                >  1:3255265/1‑150 (MQ=255)
                                                       gTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGGTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCAtt                                                                                                >  2:3580896/1‑150 (MQ=255)
                                                                  tGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGt                                                                                     <  1:4278741/150‑1 (MQ=255)
                                                                   gCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTc                                                                                    <  2:2487347/150‑1 (MQ=255)
                                                                      cTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGtt                                                                                 <  1:791387/150‑1 (MQ=255)
                                                                      cTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGtt                                                                                 <  2:3137683/150‑1 (MQ=255)
                                                                       tGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGtt                                                                                 <  1:353097/149‑1 (MQ=255)
                                                                       tGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGtt                                                                                 >  2:353097/1‑149 (MQ=255)
                                                                       tGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTg                                                                                >  2:188076/1‑150 (MQ=255)
                                                                       tGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTg                                                                                <  2:114726/150‑1 (MQ=255)
                                                                       tGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTg                                                                                >  2:4915059/1‑150 (MQ=255)
                                                                       tGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTg                                                                                >  1:4065464/1‑150 (MQ=255)
                                                                         cAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGc                                                                              >  1:2025841/1‑150 (MQ=255)
                                                                         cAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGc                                                                              <  2:4334295/150‑1 (MQ=255)
                                                                         cAAGGTCGTGCGCAGGCTCCAGTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGc                                                                              <  2:4501462/150‑1 (MQ=255)
                                                                          aaGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTgg                                                                               <  1:2191607/148‑1 (MQ=255)
                                                                          aaGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTgg                                                                               >  2:2191607/1‑148 (MQ=255)
                                                                                 tGCGCAGGCTCCATTGCACCCCCAACCCCGATCGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGc                                                                      <  2:5082997/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)
                                                                                        gCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGAt                                                               >  2:57924/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)
                                                                                          tCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCa                                                             >  2:237871/1‑150 (MQ=255)
                                                                                          tCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGAGCa                                                             >  2:4162521/1‑150 (MQ=255)
                                                                                              ttGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAAcc                                                         >  1:2353075/1‑150 (MQ=255)
                                                                                              ttGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCAGCCTCGGGCAGCAGGCGTTCCAAGGGCCAGGCCCCCAGGCCGATTCCCCGACCCTGATAGCCATTGCGCGCACGGTCGGGGGCCAGCTCGCAGCCGTTCaagcg                                                         >  2:4210700/1‑147 (MQ=255)
                                                                                                 cACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTc                                                      <  1:565422/150‑1 (MQ=255)
                                                                                                 cACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTc                                                      >  2:3599622/1‑150 (MQ=255)
                                                                                                  aCCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCACCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCg                                                     <  2:2680563/150‑1 (MQ=255)
                                                                                                    ccccAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGcc                                                   >  1:3430353/1‑150 (MQ=255)
                                                                                                    ccccAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGcc                                                   >  1:3373184/1‑150 (MQ=255)
                                                                                                      ccAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCa                                                 <  1:514959/150‑1 (MQ=255)
                                                                                                       cAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAg                                                >  1:3646845/1‑150 (MQ=255)
                                                                                                           cccGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATg                                            <  2:3780946/150‑1 (MQ=255)
                                                                                                           cccGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATg                                            <  1:3824625/150‑1 (MQ=255)
                                                                                                             cGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGc                                          <  2:1979229/150‑1 (MQ=255)
                                                                                                              gATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGcc                                         >  1:3663087/1‑150 (MQ=255)
                                                                                                                    gcTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTg                                   >  1:475852/1‑150 (MQ=255)
                                                                                                                     cTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTgc                                  <  2:3413061/150‑1 (MQ=255)
                                                                                                                        gctgcAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCa                               <  2:4888729/150‑1 (MQ=255)
                                                                                                                         ctgcAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCAc                              <  2:1636170/150‑1 (MQ=255)
                                                                                                                          tgcAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACg                             >  1:4442632/1‑150 (MQ=255)
                                                                                                                            cAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAGACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGcc                           >  2:287725/1‑150 (MQ=255)
                                                                                                                            cAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGcc                           <  2:743671/150‑1 (MQ=255)
                                                                                                                            cAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGcc                           >  2:1412250/1‑150 (MQ=255)
                                                                                                                                   tGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAgg                    >  2:3327791/1‑150 (MQ=255)
                                                                                                                                     gCCTGCTGCCCGAGGCGGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCg                  >  1:4717894/1‑150 (MQ=255)
                                                                                                                                     gCCTGCTGCCCGAGGCGGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCg                  >  1:4945363/1‑150 (MQ=255)
                                                                                                                                       cTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGtt                >  1:1386806/1‑150 (MQ=255)
                                                                                                                                       cTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGtt                >  2:2347965/1‑150 (MQ=255)
                                                                                                                                        tGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTg               >  1:504488/1‑150 (MQ=255)
                                                                                                                                          ccgccTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAATCTCGCAGCCGTTCAATCCGTCGCCCAGAATGGCCACTCTGCGCGCGCCATCAAGGCGTTGGc             >  2:3936852/1‑150 (MQ=255)
                                                                                                                                          ccgccTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGc             >  1:893990/1‑150 (MQ=255)
                                                                                                                                          ccgccTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGc             <  2:2792869/150‑1 (MQ=255)
                                                                                                                                             ccTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGCGGa          >  2:3785627/1‑150 (MQ=255)
                                                                                                                                             ccTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGCGGa          >  1:3144421/1‑150 (MQ=255)
                                                                                                                                              cTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGCGGaa         <  2:1387015/150‑1 (MQ=255)
                                                                                                                                              cTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGCGGaa         >  1:596527/1‑150 (MQ=255)
                                                                                                                                               tCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGCGGAAc        <  1:1565864/150‑1 (MQ=255)
                                                                                                                                                 gggCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGCGGAACGc      >  2:1230758/1‑150 (MQ=255)
                                                                                                                                                 gggCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGCGGAACGc      >  2:3179417/1‑150 (MQ=255)
                                                                                                                                                 aggcggcAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGCGGAACGc      >  2:2719349/2‑150 (MQ=255)
                                                                                                                                                   gcagcaGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGCGGAACGCCt    <  2:1593989/150‑1 (MQ=255)
                                                                                                                                                    cagcagGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGCGGAACGCCTg   <  1:3302581/150‑1 (MQ=255)
                                                                                                                                                     agcagGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCCCTCTGCGCACGCCATCAAGGCGTTGGCGGAACGCCTgg  <  2:4821982/150‑1 (MQ=255)
                                                                                                                                                     agcagGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGCGGAACGCCTgg  <  2:3449447/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAGGTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGCGGAACGCCTGG  >  NC_002947/6123206‑6123504

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: