Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,181,925 C→T G170G (GGG→GGA U712_10925 ← Menaquinol‑cytochrome c reductase cytochrome b subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,181,9250CT98.6% 222.3 / ‑4.8 73G170G (GGG→GGAU712_10925Menaquinol‑cytochrome c reductase cytochrome b subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (39/33);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

ATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAGAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCAGGGTGCCCCGCAAGGAGTGTTTTGACTTGCGTCCCGATCAAAGGCGTGGCTTCAGCGATCTGCAGCCCGACTTTTGTTGCAAACAGAGCCTTCATGTCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGCCCGAGCATGACAAAAAAGATCAGTACACCGACAATCCAG  >  CP006881/2181794‑2182060
                                                                                                                                   |                                                                                                                                       
ataataAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGATCCGATCa                                                                                                                                  <  8:327493/139‑1 (MQ=255)
 taataaAATGTGCAGCCATCAGCCCGACCAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCa                                                                                                                                  >  2:274586/1‑138 (MQ=255)
 taataaAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGATGATAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCGTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCa                                                                                                                                  >  6:20722/1‑138 (MQ=255)
      aaaTGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGct                                                                                                                            >  8:175446/1‑138 (MQ=255)
      aaaTGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTGTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGc                                                                                                                             <  5:368828/138‑1 (MQ=255)
           tgCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCt                                                                                                                       >  3:381148/1‑139 (MQ=255)
           tgCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCt                                                                                                                       >  5:333233/1‑139 (MQ=255)
           tgCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCt                                                                                                                       >  2:184577/1‑139 (MQ=255)
            gCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCtt                                                                                                                      >  8:115332/1‑139 (MQ=255)
             cAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCtt                                                                                                                      <  3:243863/138‑1 (MQ=255)
                     aGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTcg                                                                                                                    <  3:440753/132‑2 (MQ=255)
                     aGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTcg                                                                                                                    >  4:440753/1‑131 (MQ=255)
                          gAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAgt                                                                                                        >  7:128301/1‑138 (MQ=255)
                          gAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAgt                                                                                                        >  1:416097/1‑138 (MQ=255)
                           aaCAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAgtc                                                                                                       >  2:94355/1‑137 (MQ=255)
                               aTGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGAc                                                                                                   >  1:391875/1‑139 (MQ=255)
                                tGCGGCAGGAAGAAAGAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCAGGGTGCACCGCAAGGAGTGTTTTGACTTGCGTCCCGATCAAAGGCGTGGCTTCAGCGATCTGCAGCCCGACt                                                                                                  >  6:254084/1‑139 (MQ=255)
                                      aGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTttgtt                                                                                            <  1:94355/139‑1 (MQ=255)
                                       ggAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTttgttg                                                                                           <  8:128301/139‑1 (MQ=255)
                                             aaaaaTACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACa                                                                                     <  8:410791/136‑1 (MQ=255)
                                                aaTACATGGATCGCAAAAAACCTAGTCCGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACagag                                                                                  <  5:56977/139‑1 (MQ=255)
                                                         aTCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCCATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCAtat                                                                         >  2:341532/1‑137 (MQ=255)
                                                                 aaaCCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGa                                                                 <  7:76827/139‑1 (MQ=255)
                                                                         tCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCGTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTa                                                         <  4:385146/139‑1 (MQ=255)
                                                                           aGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGc                                                       <  4:381148/139‑1 (MQ=255)
                                                                               tttGGGCGCCTACAATCTCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt                                                   <  5:32451/139‑1 (MQ=255)
                                                                               tttGGGCGCCTACAATATCCGGATGCCCTTCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt                                                   <  2:141754/139‑1 (MQ=255)
                                                                               tttGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAATCAAGTAGCCTGt                                                   <  4:217985/139‑1 (MQ=255)
                                                                               tttGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt                                                   <  2:416097/139‑1 (MQ=255)
                                                                               tttGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt                                                   <  6:333233/139‑1 (MQ=255)
                                                                               tttGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTGTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt                                                   <  1:341532/139‑1 (MQ=255)
                                                                               caggggCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt                                                   <  1:405006/136‑1 (MQ=255)
                                                                                  gggCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt                                                   <  3:160072/136‑1 (MQ=255)
                                                                                  gggCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt                                                   >  4:160072/1‑136 (MQ=255)
                                                                                   ggCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAaag                                               >  4:105073/1‑139 (MQ=255)
                                                                                       ccTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAaagccaa                                           >  3:34132/1‑139 (MQ=255)
                                                                                       ccTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAaagcc                                             >  3:428567/1‑137 (MQ=255)
                                                                                       ccTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAaagcc                                             >  5:158221/1‑137 (MQ=255)
                                                                                           cAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAaagccaaga                                         <  1:184577/137‑2 (MQ=255)
                                                                                           cAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGacc                                       <  8:31012/139‑1 (MQ=255)
                                                                                             atatCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt                                                   >  4:283274/1‑125 (MQ=255)
                                                                                             atatCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt                                                   <  3:283274/125‑1 (MQ=255)
                                                                                              tatCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt                                                   <  7:287409/124‑1 (MQ=255)
                                                                                              tatCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt                                                   >  8:287409/1‑124 (MQ=255)
                                                                                                  cggatGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCGGACTGTTGTTGCAAACAGAGGCTTCAGATCCCATGGATGCAAGTCGCCTGTATAGCCATGACCGAGCaga                                >  2:389060/5‑137 (MQ=255)
                                                                                                  cggatGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCAt                                 <  3:445629/134‑1 (MQ=255)
                                                                                                    gatGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACagag                                                                                  <  8:271782/85‑1 (MQ=255)
                                                                                                    gatGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACagag                                                                                  >  7:271782/3‑87 (MQ=255)
                                                                                                     atGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAAc                              >  3:331009/2‑138 (MQ=255)
                                                                                                     atGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAAc                              >  6:280713/2‑138 (MQ=255)
                                                                                                     atGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAAc                              >  2:287788/2‑138 (MQ=255)
                                                                                                     atGCCCTGCAAGGAGATTTTAAACTTGCGTTCCGAAAAAAGCAGTGGCTTCGGCGATCTGCAGTACGACTTTTGTTGCAAACAGAGCCTTCATATCACAAGGAAGCAAGTAGCCTGTAAAGACAAGACCTAGCATAAc                              >  7:45134/2‑138 (MQ=255)
                                                                                                      cgCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAg                                    <  7:97539/130‑1 (MQ=255)
                                                                                                      agCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAg                                    >  8:97539/2‑131 (MQ=255)
                                                                                                        ccctGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACaaa                           >  3:416215/1‑138 (MQ=255)
                                                                                                        ccctGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACaaa                           >  6:92901/1‑138 (MQ=255)
                                                                                                             cAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGCCCGAGCATAACaaaa                          >  4:31929/1‑134 (MQ=255)
                                                                                                               acgAGCGTTTTAACTTGATTTCAGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAAAAGAGCCTTCATATCACAAGGAAGCAAGTAGCCTGTAAAGCAAAGACAGAGCATAACAAAAAAGCTCa                   >  8:98031/3‑139 (MQ=255)
                                                                                                                  agcgtTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAATCCAAGACCGAGCATAACAAAAAAGATCAgca                >  8:394523/4‑137 (MQ=255)
                                                                                                                    cgtTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCTACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGAAAAGCCAAGACCGAGCATAACAAAAAAGATCAGcaca              >  5:27674/2‑139 (MQ=255)
                                                                                                                    cgtTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGcaca              >  1:297778/2‑139 (MQ=255)
                                                                                                                    cgtTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCac               <  7:175446/137‑1 (MQ=255)
                                                                                                                         taaCTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGa          <  5:112663/136‑1 (MQ=255)
                                                                                                                         taaCTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGa          <  4:34132/136‑1 (MQ=255)
                                                                                                                         taaCTTGCGTTCCGAGCAAAGGCGGGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGa          <  2:126024/136‑1 (MQ=255)
                                                                                                                            cTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAAt      <  1:205034/139‑1 (MQ=255)
                                                                                                                                cGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACa                                                                                     <  2:245963/56‑1 (MQ=255)
                                                                                                                                cGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACa                                                                                     >  1:245963/1‑56 (MQ=255)
                                                                                                                                cGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAAt      <  7:222790/135‑1 (MQ=255)
                                                                                                                                cGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAAt      >  8:222790/1‑135 (MQ=255)
                                                                                                                                cGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAg  >  5:337123/1‑139 (MQ=255)
                                                                                                                                cGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGAATAACAAAAAAGATCAGCACACCGACAATCCAg  >  1:291754/1‑139 (MQ=255)
                                                                                                                                cGTTCCGATCAAAGGCGTCGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATACCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAg  >  7:144460/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                       
ATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAGAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCAGGGTGCCCCGCAAGGAGTGTTTTGACTTGCGTCCCGATCAAAGGCGTGGCTTCAGCGATCTGCAGCCCGACTTTTGTTGCAAACAGAGCCTTCATGTCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGCCCGAGCATGACAAAAAAGATCAGTACACCGACAATCCAG  >  CP006881/2181794‑2182060

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