Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,127,019 T→A 5.6% K46* (AAA→TAA) ‡ PP_5371 ← rubredoxin/rubredoxin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,127,0190TA 5.6% 1334.4 / 43.8 394K46* (AAA→TAA) ‡PP_5371rubredoxin/rubredoxin reductase
Reads supporting (aligned to +/- strand):  ref base T (169/203);  new base A (5/17);  total (174/220)
Fisher's exact test for biased strand distribution p-value = 4.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 407 total aligned reads displayed.

CTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTC  >  NC_002947/6126870‑6127168
                                                                                                                                                     |                                                                                                                                                     
cTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCttt                                                                                                                                                       <  1:1188443/150‑1 (MQ=255)
cTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCttt                                                                                                                                                       >  1:7168861/1‑150 (MQ=255)
  gtgCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGc                                                                                                                                                     >  2:4834456/1‑150 (MQ=255)
   tgCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGcc                                                                                                                                                    >  1:2140707/1‑150 (MQ=255)
    gCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCg                                                                                                                                                   >  2:1391661/1‑150 (MQ=255)
     ccAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGa                                                                                                                                                  >  2:8652104/1‑150 (MQ=255)
      cAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGAc                                                                                                                                                 <  2:4138783/150‑1 (MQ=255)
       aGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGAcc                                                                                                                                                <  2:5110233/150‑1 (MQ=255)
       aGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGAcc                                                                                                                                                <  2:5806565/150‑1 (MQ=255)
       aGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGAcc                                                                                                                                                <  2:707982/150‑1 (MQ=255)
        ggCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCTGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCAAAGCCGActt                                                                                                                                               <  2:3831199/150‑3 (MQ=255)
         gCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCtt                                                                                                                                              >  1:2323743/1‑150 (MQ=255)
          ccGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTg                                                                                                                                             <  1:7543734/150‑1 (MQ=255)
          ccGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTg                                                                                                                                             <  1:7334086/150‑1 (MQ=255)
           cGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAGGTCGGCAAAGCCGACTTTga                                                                                                                                            <  1:8206552/150‑2 (MQ=255)
            gTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCa                                                                                                                                           <  1:1828386/150‑1 (MQ=255)
            gTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCAAAGCCGACTTTgaa                                                                                                                                           <  1:1384502/150‑3 (MQ=255)
             tAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCAc                                                                                                                                          <  2:576657/150‑1 (MQ=255)
                 cGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAGGTCGGCAAAGCCGACTTTgaaatgat                                                                                                                                      <  2:8486836/150‑8 (MQ=255)
                 cGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAg                                                                                                                                      <  2:3283111/150‑1 (MQ=255)
                 cGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAg                                                                                                                                      <  1:6614606/150‑1 (MQ=255)
                       ggCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATa                                                                                                                                <  2:8713529/150‑1 (MQ=255)
                       ggCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATa                                                                                                                                >  2:4605257/1‑150 (MQ=255)
                       ggCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATa                                                                                                                                <  1:2626778/150‑1 (MQ=255)
                         ccGCTACCGATGATGACCACCGGGTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGc                                                                                                                              <  1:7415177/150‑1 (MQ=255)
                           gCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCa                                                                                                                            <  2:1287442/150‑1 (MQ=255)
                           gCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCa                                                                                                                            <  2:4003604/150‑1 (MQ=255)
                                  atgatgACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCg                                                                                                                     >  2:7191761/1‑150 (MQ=255)
                                  atgatgACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCg                                                                                                                     >  2:5635767/1‑150 (MQ=255)
                                    gatgaCCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCggg                                                                                                                   <  1:1070900/150‑1 (MQ=255)
                                      tgaCCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCa                                                                                                                 >  2:4963748/1‑150 (MQ=255)
                                       gaCCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCAc                                                                                                                <  1:881430/150‑1 (MQ=255)
                                         ccaccGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGt                                                                                                              >  1:6918107/1‑150 (MQ=255)
                                         ccaccGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGt                                                                                                              >  1:1937323/1‑150 (MQ=255)
                                           accGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGtcc                                                                                                            >  1:640140/1‑150 (MQ=255)
                                              ggTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGtcctcc                                                                                                         <  1:3980099/150‑1 (MQ=255)
                                                           gAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCgg                                                                                            >  2:3106022/1‑150 (MQ=255)
                                                              cGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCggggg                                                                                         <  2:1128601/150‑1 (MQ=255)
                                                               ggCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGc                                                                                        <  1:334116/150‑1 (MQ=255)
                                                                  aaCCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAAt                                                                                     >  1:1422309/1‑150 (MQ=255)
                                                                     cAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGcc                                                                                  >  2:8368436/1‑150 (MQ=255)
                                                                     cAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGcc                                                                                  <  1:6317040/150‑1 (MQ=255)
                                                                     cAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGcc                                                                                  <  1:2552618/150‑1 (MQ=255)
                                                                     cAGTGGTTGCGCCTGGAGAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGcc                                                                                  <  2:6269749/150‑1 (MQ=255)
                                                                       gTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCcgt                                                                                >  2:2129843/1‑150 (MQ=255)
                                                                       gTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCcgt                                                                                <  1:3498008/150‑1 (MQ=255)
                                                                       gTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCcgt                                                                                >  1:4900202/1‑150 (MQ=255)
                                                                       gTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCcgt                                                                                >  1:3638400/1‑150 (MQ=255)
                                                                              cgcCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCggg                                                                         <  1:3795404/150‑1 (MQ=255)
                                                                               gcCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGc                                                                        >  1:3021945/1‑150 (MQ=255)
                                                                                  tGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCggtgacc                                                                     >  1:8538434/1‑145 (MQ=255)
                                                                                    gATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAAcc                                                                   >  1:326510/1‑150 (MQ=255)
                                                                                       aaCCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCttt                                                                <  1:5502243/150‑1 (MQ=255)
                                                                                         ctggCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTgg                                                              >  2:2565240/3‑150 (MQ=255)
                                                                                         ccGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTgg                                                              <  2:592766/150‑1 (MQ=255)
                                                                                         ccGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTgg                                                              <  1:5642778/150‑1 (MQ=255)
                                                                                         ccGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTgg                                                              <  1:6798238/150‑1 (MQ=255)
                                                                                          cGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGc                                                             >  1:8453792/1‑150 (MQ=255)
                                                                                          cGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGc                                                             <  2:7666784/150‑1 (MQ=255)
                                                                                             cgcgGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTtc                                                          <  1:4554901/150‑1 (MQ=255)
                                                                                                    gCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAg                                                   <  1:176711/150‑1 (MQ=255)
                                                                                                       aCTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATc                                                <  1:4598107/150‑1 (MQ=255)
                                                                                                        ctgctgCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCa                                               <  1:5266161/150‑1 (MQ=255)
                                                                                                        ctgctgCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCa                                               <  1:7847394/150‑1 (MQ=255)
                                                                                                          gctgcGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGa                                             <  1:6207182/150‑1 (MQ=255)
                                                                                                           ctgcGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGAc                                            <  2:3014406/150‑1 (MQ=255)
                                                                                                           ctgcGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGAc                                            <  2:6712282/150‑1 (MQ=255)
                                                                                                                    gggTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAgatga                                   <  2:6824150/150‑1 (MQ=255)
                                                                                                                    gggTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAgatga                                   <  2:6142525/150‑1 (MQ=255)
                                                                                                                     ggTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAgatgat                                  >  2:7549904/1‑150 (MQ=255)
                                                                                                                     ggTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGGCCTCCCAGCGGGTCCCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAgatgat                                  >  2:2269555/1‑150 (MQ=255)
                                                                                                                      gTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGatgatg                                 <  2:5187655/150‑1 (MQ=255)
                                                                                                                        gCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCa                               <  1:8638288/150‑1 (MQ=255)
                                                                                                                          tGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCaga                             <  2:1858258/150‑1 (MQ=255)
                                                                                                                           gAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCagag                            >  1:4312902/1‑150 (MQ=255)
                                                                                                                             gTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTCCCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGcc                          >  2:3671455/1‑150 (MQ=255)
                                                                                                                              ttCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCa                         >  2:2861525/1‑150 (MQ=255)
                                                                                                                              ttCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCa                         <  2:2747717/150‑1 (MQ=255)
                                                                                                                               tCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAg                        <  2:5294570/150‑1 (MQ=255)
                                                                                                                               tCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAg                        <  2:8094611/150‑1 (MQ=255)
                                                                                                                                 gATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGt                      >  2:600243/1‑150 (MQ=255)
                                                                                                                                 gATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGt                      >  1:7048108/1‑150 (MQ=255)
                                                                                                                                 gATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGt                      >  1:5950482/1‑150 (MQ=255)
                                                                                                                                  ctgagtgcAAGGTCGGCAAAGCCGACTTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGtt                     >  1:1593864/8‑150 (MQ=255)
                                                                                                                                  aTCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGtt                     >  1:5063308/1‑150 (MQ=255)
                                                                                                                                     aTTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTc                  <  1:2884577/150‑1 (MQ=255)
                                                                                                                                     aTTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTc                  <  2:8052303/150‑1 (MQ=255)
                                                                                                                                      tttCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTct                 >  2:6464174/1‑149 (MQ=255)
                                                                                                                                      tttCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCa                 <  1:8112238/150‑1 (MQ=255)
                                                                                                                                      tttCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCa                 >  2:1649167/1‑150 (MQ=255)
                                                                                                                                       ttCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGTAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCAt                >  1:5373810/1‑150 (MQ=255)
                                                                                                                                           aaGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTAtg            >  1:5857200/1‑150 (MQ=255)
                                                                                                                                           aaGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGAAGAGCCAGGTTTTCATTAtg            >  2:7254858/1‑150 (MQ=255)
                                                                                                                                              tCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTAtgctg         >  2:4364452/1‑150 (MQ=255)
                                                                                                                                               cGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTAtgctgc        <  1:8049674/150‑1 (MQ=255)
                                                                                                                                                ggCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGcc       <  2:8332707/150‑1 (MQ=255)
                                                                                                                                                  cTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGcctc     <  2:554049/150‑1 (MQ=255)
                                                                                                                                                   tttGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGcctcc    >  1:2461656/1‑150 (MQ=255)
                                                                                                                                                    ttGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGcctcct   <  2:3538446/150‑1 (MQ=255)
                                                                                                                                                    ttGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGcctcct   <  1:8322718/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCCTCCTC  >  NC_002947/6126870‑6127168

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: