Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,123,351 G→A 5.4% P259P (CCC→CCT PP_5371 ← rubredoxin/rubredoxin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,123,3510GA 5.4% 528.1 / 12.3 167P259P (CCC→CCTPP_5371rubredoxin/rubredoxin reductase
Reads supporting (aligned to +/- strand):  ref base G (86/72);  new base A (4/5);  total (90/77)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 174 total aligned reads displayed.

AGCACCAGGTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGCGGAACGC  >  NC_002947/6123202‑6123500
                                                                                                                                                     |                                                                                                                                                     
aGCACCAGGTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGAAAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCg                                                                                                                                                       >  1:4241479/1‑150 (MQ=255)
 gCACCAGGTCTGCACCGACGACCTGACCGTCGGCTAGTGTCAGTCGGTAGCCTTCGCCGTCAGGATCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCTCCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCgg                                                                                                                                                      <  2:503130/150‑1 (MQ=255)
  cACCAGGTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCggg                                                                                                                                                     <  2:3872616/150‑1 (MQ=255)
    ccAGGTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGgca                                                                                                                                                   <  2:1586495/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)
          cTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGc                                                                                                                                             <  1:2320008/150‑1 (MQ=255)
           tGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCTGCCCGAGGCGGCAGGcc                                                                                                                                            <  2:2593784/150‑2 (MQ=255)
            gCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCTGCCCGAGGCGGCAGGcca                                                                                                                                           <  2:796997/150‑3 (MQ=255)
            gCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCTGCCCGAGGCGGCAGGcca                                                                                                                                           <  1:246879/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)
                 gacgacCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCgctatcg                                                                                                                                      <  1:5067788/150‑7 (MQ=255)
                 gacgacCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCAt                                                                                                                                      >  2:5005419/1‑150 (MQ=255)
                    gacCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTACCATggg                                                                                                                                   >  1:1593989/1‑150 (MQ=255)
                      ccTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGACCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGcc                                                                                                                                 >  2:3623011/1‑150 (MQ=255)
                      ccTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGACCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGcc                                                                                                                                 >  2:1907907/1‑150 (MQ=255)
                      ccTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCAATGCACCCCCAACCCCGATAGAGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGAACCAAGGGcc                                                                                                                                 >  2:565422/1‑150 (MQ=255)
                          aCCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGc                                                                                                                             >  1:2248122/1‑150 (MQ=255)
                            cGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCAc                                                                                                                           >  1:3007948/1‑150 (MQ=255)
                            cGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCAc                                                                                                                           >  1:3083032/1‑150 (MQ=255)
                                       gTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGTCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGAtt                                                                                                                >  2:418610/1‑150 (MQ=255)
                                       gTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGAtt                                                                                                                <  1:440431/150‑1 (MQ=255)
                                        tCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTa                                                                                                               >  2:1266847/1‑150 (MQ=255)
                                          aGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTAcc                                                                                                             <  1:3153771/150‑1 (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:4191159/150‑1 (MQ=255)
                                                     ttCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATa                                                                                                  >  2:874796/1‑150 (MQ=255)
                                                     ttCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATa                                                                                                  <  1:4192442/150‑1 (MQ=255)
                                                      tCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAg                                                                                                 >  1:551047/1‑150 (MQ=255)
                                                       cGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGc                                                                                                <  1:3755428/150‑1 (MQ=255)
                                                           gTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGGTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCAtt                                                                                            >  2:3580896/1‑150 (MQ=255)
                                                                    gATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAg                                                                                   >  2:3959222/1‑150 (MQ=255)
                                                                       gCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTc                                                                                <  2:2487347/150‑1 (MQ=255)
                                                                       gCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTc                                                                                <  1:3865041/150‑1 (MQ=255)
                                                                         aCTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGt                                                                              <  1:2979667/150‑1 (MQ=255)
                                                                         aCTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGt                                                                              >  2:1387750/1‑150 (MQ=255)
                                                                          cTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGtt                                                                             >  1:484058/1‑150 (MQ=255)
                                                                          cTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGtt                                                                             >  2:3185710/1‑150 (MQ=255)
                                                                          cTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGtt                                                                             <  2:3139672/150‑1 (MQ=255)
                                                                          cTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGtt                                                                             <  2:3137683/150‑1 (MQ=255)
                                                                           tGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGtt                                                                             >  2:353097/1‑149 (MQ=255)
                                                                           tGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTg                                                                            >  2:4915059/1‑150 (MQ=255)
                                                                           tGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTg                                                                            >  1:3360742/1‑150 (MQ=255)
                                                                           tGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTg                                                                            >  2:188076/1‑150 (MQ=255)
                                                                           tGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTg                                                                            >  1:4065464/1‑150 (MQ=255)
                                                                           tGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTg                                                                            <  2:114726/150‑1 (MQ=255)
                                                                             cAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGc                                                                          <  2:4334295/150‑1 (MQ=255)
                                                                             cAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGc                                                                          >  1:2025841/1‑150 (MQ=255)
                                                                             cAAGGTCGTGCGCAGGCTCCAGTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGc                                                                          <  2:4501462/150‑1 (MQ=255)
                                                                              aaGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTgg                                                                           <  1:2191607/148‑1 (MQ=255)
                                                                              aaGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTgg                                                                           >  2:2191607/1‑148 (MQ=255)
                                                                              aaGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCa                                                                         >  1:75597/1‑150 (MQ=255)
                                                                                       cgcAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAg                                                                <  1:1814847/150‑1 (MQ=255)
                                                                                          aGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCg                                                             >  1:3884147/1‑150 (MQ=255)
                                                                                          aGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCg                                                             <  1:1005680/150‑1 (MQ=255)
                                                                                          aGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCg                                                             >  2:1462935/1‑150 (MQ=255)
                                                                                             cTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATc                                                          >  1:1398/1‑150 (MQ=255)
                                                                                              tCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCa                                                         >  2:3634969/1‑150 (MQ=255)
                                                                                              tCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCa                                                         >  2:724244/1‑150 (MQ=255)
                                                                                               ccATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCaa                                                        >  2:3824625/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)
                                                                                                      aCCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCACCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCg                                                 <  2:2680563/150‑1 (MQ=255)
                                                                                                        ccccAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCACACTCGCAGCCGATCAAACCGTCGcc                                               >  2:4075493/1‑150 (MQ=255)
                                                                                                        ccccAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGcc                                               >  1:3373184/1‑150 (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)
                                                                                                               cccGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAACGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATg                                        <  2:1161821/150‑1 (MQ=255)
                                                                                                                 cGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGc                                      >  2:2482866/1‑150 (MQ=255)
                                                                                                                 cGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGc                                      >  1:567383/1‑150 (MQ=255)
                                                                                                                 cGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGc                                      <  2:1979229/150‑1 (MQ=255)
                                                                                                                       cgcTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCAGTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCActct                                <  2:5060508/150‑1 (MQ=255)
                                                                                                                         cTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTgc                              <  2:3413061/150‑1 (MQ=255)
                                                                                                                             ctgcAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCAc                          <  2:1636170/150‑1 (MQ=255)
                                                                                                                               gcAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGc                        <  2:484058/150‑1 (MQ=255)
                                                                                                                                cAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGcc                       >  1:1319220/1‑150 (MQ=255)
                                                                                                                                 aaGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCa                      >  1:4156081/1‑150 (MQ=255)
                                                                                                                                    ggaacgCCTGCTGCCCGAGGCGGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCa                   >  1:2785127/6‑150 (MQ=255)
                                                                                                                                         gCCTGCTGCCCGAGGCGGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCg              >  1:4717894/1‑150 (MQ=255)
                                                                                                                                         gCCTGCTGCCCGAGGCGGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCg              >  1:4945363/1‑150 (MQ=255)
                                                                                                                                         gCCTGCTGCCCGAGGCGGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCg              >  1:4717980/1‑150 (MQ=255)
                                                                                                                                           cTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGtt            >  2:2347965/1‑150 (MQ=255)
                                                                                                                                           cTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGtt            >  1:1387015/1‑150 (MQ=255)
                                                                                                                                           cTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGtt            >  2:740530/1‑150 (MQ=255)
                                                                                                                                            tGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTg           >  1:504512/1‑150 (MQ=255)
                                                                                                                                            tGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTg           >  1:504488/1‑150 (MQ=255)
                                                                                                                                              ccgccTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGc         <  2:2792869/150‑1 (MQ=255)
                                                                                                                                              ccgccTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGc         >  1:893990/1‑150 (MQ=255)
                                                                                                                                                 ccTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCCTCAAGGCGTTGGCGGa      >  2:3528238/1‑150 (MQ=255)
                                                                                                                                                  cTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGCGGaa     <  2:1387015/150‑1 (MQ=255)
                                                                                                                                                   tCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGCGGAAc    >  1:3417505/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCACCAGGTCTGCACCGACGACCTGACCGTCGGCAAGGGTCAGTCGGTAGCCTTCGCCGTCCGGCTCGATGCACTGCAAGGTCGTGCGCAGGCTCCATTGCACCCCCAACCCCGATAGCGCTTGCTGCAAGTGCTGGCCTGCCGCCTCGGGCAGCAGGCGTTCCATGGGCCAGGCACCCAGGCCGATTACCCGAACCTGATAGCCATTGCTCGCAAGGTCGTTGGCAAACTCGCAGCCGATCAAACCGTCGCCCAGAATGGCCACTCTGCGCACGCCATCAAGGCGTTGGCGGAACGC  >  NC_002947/6123202‑6123500

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: