Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,182,021 C→A G138G (GGG→GGT 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,182,0210CA98.8% 251.3 / NA 85G138G (GGG→GGTU712_10925Menaquinol‑cytochrome c reductase cytochrome b subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (48/36);  total (49/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAATATCAGGGTGCCCCGCAAGGAGTGTTTTGACTTGCGTCCCGATCAAAGGCGTGGCTTCAGCGATCTGCAGCCCGACTTTTGTTGCAAACAGAGCCTTCATGTCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGCCCGAGCATGACAAAAAAGATCAGTACACCGACAATCCAGTTCAGCTCGCGCGGCTTTTTATACGCCCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCACCCCAGTGGTGCA  >  CP006881/2181885‑2182157
                                                                                                                                        |                                                                                                                                        
cAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGacc                                                                                                                                        <  8:31012/139‑1 (MQ=255)
       cggatGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCGGACTGTTGTTGCAAACAGAGGCTTCAGATCCCATGGATGCAAGTCGCCTGTATAGCCATGACCGAGCaga                                                                                                                                 >  2:389060/5‑137 (MQ=255)
       cggatGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCAt                                                                                                                                  <  3:445629/134‑1 (MQ=255)
          atGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAAc                                                                                                                               >  6:280713/2‑138 (MQ=255)
          atGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAAc                                                                                                                               >  2:287788/2‑138 (MQ=255)
          atGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAAc                                                                                                                               >  3:331009/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                                                                                                                            >  6:92901/1‑138 (MQ=255)
             ccctGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACaaa                                                                                                                            >  3:416215/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                                                                                                           <  4:34132/136‑1 (MQ=255)
                              taaCTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGa                                                                                                           <  5:112663/136‑1 (MQ=255)
                              taaCTTGCGTTCCGAGCAAAGGCGGGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGa                                                                                                           <  2:126024/136‑1 (MQ=255)
                                 cTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAAt                                                                                                       <  1:205034/139‑1 (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)
                                      gttccGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGt                                                                                                  >  5:115781/4‑139 (MQ=255)
                                            aTCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCt                                                                                            >  1:268589/1‑139 (MQ=255)
                                                aaGGCGTGGCTGCGGCTATCTGCAGTCCTACTTTTGTTGCAAACAGATCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTcgcg                                                                                        <  2:291754/139‑1 (MQ=255)
                                                 aGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTAGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTcgcgc                                                                                       <  3:172371/139‑1 (MQ=255)
                                                      tGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGttt                                                                                  >  5:48002/1‑139 (MQ=255)
                                                      tGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGttt                                                                                  >  5:206196/1‑139 (MQ=255)
                                                      tGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGttt                                                                                  >  5:117808/1‑139 (MQ=255)
                                                      tGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGttt                                                                                  >  2:111950/1‑139 (MQ=255)
                                                      tGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGttt                                                                                  >  7:390949/1‑139 (MQ=255)
                                                              gCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACg                                                                          <  8:390949/139‑1 (MQ=255)
                                                              gCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACg                                                                          <  5:280713/139‑1 (MQ=255)
                                                              gCGATCTGCAGTCCGACGTTTGTTGCATACAGAGCCTTAATATAGCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACg                                                                          <  1:291888/139‑1 (MQ=255)
                                                                  tCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACg                                                                          <  1:296234/135‑1 (MQ=255)
                                                                  tCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACg                                                                          >  2:296234/1‑135 (MQ=255)
                                                                         tccGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGaaa                                                               <  4:331009/138‑1 (MQ=255)
                                                                           cGACTTTTGTTGCAAACAGAGCCTTCATATCCCAATGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGaaaaa                                                             <  1:389060/139‑1 (MQ=255)
                                                                           cGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCGGTAAAGCCAAGACCGCGCATAACAAAACAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGaaaac                                                             >  6:54308/1‑138 (MQ=255)
                                                                           cGACTTTTGTAGCAAACAGAGCCTACATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGATCAGCTCGCGCGGTTTTTTATACGCGCCTTGGaaaaa                                                             <  5:20722/139‑1 (MQ=255)
                                                                              cTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGAc                                                          >  1:195179/1‑139 (MQ=255)
                                                                                  tgttgCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCt                                                       <  1:287788/138‑1 (MQ=255)
                                                                                       cAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCg                                                  >  6:198253/1‑138 (MQ=255)
                                                                                       cAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCg                                                  <  1:274586/138‑1 (MQ=255)
                                                                                       cAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCg                                                  <  7:115332/138‑1 (MQ=255)
                                                                                              agCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTaaa                                          <  4:428567/139‑1 (MQ=255)
                                                                                              agCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTaaa                                          <  4:416215/139‑1 (MQ=255)
                                                                                              agCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTaa                                           >  8:434404/1‑138 (MQ=255)
                                                                                              agCCTTCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTaa                                           >  1:59413/1‑138 (MQ=255)
                                                                                                   tCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTcgcg                                                                                        >  2:170689/1‑88 (MQ=255)
                                                                                                   tCATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTcgcg                                                                                        <  1:170689/88‑1 (MQ=255)
                                                                                                    cATATCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTGTTTTATACGCGCGTTGGAAAAAGACCCGAAGCGTATGTAAAAAcat                                     >  2:279471/1‑138 (MQ=255)
                                                                                                     atatCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAAcatc                                    <  7:394523/135‑1 (MQ=255)
                                                                                                          ccAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATcaca                              <  6:337123/139‑1 (MQ=255)
                                                                                                          ccAAGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATcaca                              <  2:195179/139‑1 (MQ=255)
                                                                                                            aaGGAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAAt                            >  4:258937/1‑139 (MQ=255)
                                                                                                              ggAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGa                          >  2:60983/1‑139 (MQ=255)
                                                                                                              ggAAGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGa                          >  2:54305/1‑139 (MQ=255)
                                                                                                                 aGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGa                                                                                                           >  2:227772/1‑55 (MQ=255)
                                                                                                                 aGCAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGa                                                                                                           <  1:227772/55‑1 (MQ=255)
                                                                                                                   cAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAATACCCTCAGCGTATGTAAAAACATCATCACAATGACCAgg                     >  2:401211/1‑139 (MQ=255)
                                                                                                                   cAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATcaca                              >  4:436526/1‑130 (MQ=255)
                                                                                                                   cAAGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATcaca                              <  3:436526/130‑1 (MQ=255)
                                                                                                                    aaGTAGCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTGTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGc                    >  4:256149/1‑139 (MQ=255)
                                                                                                                         gCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCTCACCGACAATCCAGTTCAGCTCGCTCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCTTATGTAAAAACATCATCACAATGACCAGGCTCGCa               <  8:144460/139‑1 (MQ=255)
                                                                                                                         gCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCa               <  3:105073/139‑1 (MQ=255)
                                                                                                                         gCCTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCa               >  5:378051/1‑139 (MQ=255)
                                                                                                                          ccTGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCAc              >  2:30635/1‑139 (MQ=255)
                                                                                                                          ccTGTAAACCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCAc              <  2:297778/139‑1 (MQ=255)
                                                                                                                            tGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCa                                                     <  4:27699/98‑1 (MQ=255)
                                                                                                                            tGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCa                                                     >  1:321077/1‑98 (MQ=255)
                                                                                                                            tGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCa                                                     <  2:321077/98‑1 (MQ=255)
                                                                                                                            tGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCa                                                     >  3:27699/1‑98 (MQ=255)
                                                                                                                            tGTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCAccc            >  8:439024/1‑139 (MQ=255)
                                                                                                                             gTAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCAccc            >  3:86762/1‑138 (MQ=255)
                                                                                                                              tAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGCCCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCAccc            >  4:53710/1‑137 (MQ=255)
                                                                                                                              tAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTa                                            >  4:54016/1‑105 (MQ=255)
                                                                                                                              tAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTa                                            <  3:54016/105‑1 (MQ=255)
                                                                                                                              tAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCAccc            >  8:81114/1‑137 (MQ=255)
                                                                                                                              tAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCAccc            >  3:104970/1‑137 (MQ=255)
                                                                                                                              tAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCACCCCa          <  1:111950/139‑1 (MQ=255)
                                                                                                                              tAAAGCCAAGACCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCACCCCa          <  7:434404/139‑1 (MQ=255)
                                                                                                                                aagccaagtcCGAGCTTAAGAAAAAAGCTCAGGACACCGACTCTCGACTTCAGCTCGGGCTGTTTTTTATACGCGGCTTTGAAAAAGACCCTCTGCGTATGTAATAAAATCATCACAATGACCAGGCTTTCACCCCAgt        <  8:45134/139‑1 (MQ=255)
                                                                                                                                aagccaagacCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCg                 <  1:110236/130‑1 (MQ=255)
                                                                                                                                aagccaagacCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCg                 >  2:110236/1‑130 (MQ=255)
                                                                                                                                aagccaagacCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCACCCCAgt        <  1:30635/139‑1 (MQ=255)
                                                                                                                                 agccaagccCGAGCATGACAAAAAAGATCAGTACACCGACAATCCATTTCAGCTCGCGCGGCTTTTTATACGCCCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCACCCCAgtg       <  5:254084/139‑1 (MQ=255)
                                                                                                                                 agccaagacCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCACCCCAgtg       >  4:405477/1‑139 (MQ=255)
                                                                                                                                   ccaagacCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCg                                                  >  8:203792/1‑94 (MQ=255)
                                                                                                                                   ccaagacCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCg                                                  <  7:203792/94‑1 (MQ=255)
                                                                                                                                     aagacCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCACCCCAgtggtg    >  6:323503/1‑138 (MQ=255)
                                                                                                                                     aagacCGAGCATAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCACCCCAGTGGTGc   <  3:31929/139‑1 (MQ=255)
                                                                                                                                      agccCGAGCAAAACAAAAAAGATCAGCACACCGACAATCCAGTTCAGCTCGCGCGGTTTTTTATACGCGAGTTGGCAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCACCCCAGTGGTGCa  >  2:398348/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                        
CAATATCAGGGTGCCCCGCAAGGAGTGTTTTGACTTGCGTCCCGATCAAAGGCGTGGCTTCAGCGATCTGCAGCCCGACTTTTGTTGCAAACAGAGCCTTCATGTCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGCCCGAGCATGACAAAAAAGATCAGTACACCGACAATCCAGTTCAGCTCGCGCGGCTTTTTATACGCCCCTTGGAAAAAGACCCTCAGCGTATGTAAAAACATCATCACAATGACCAGGCTCGCACCCCAGTGGTGCA  >  CP006881/2181885‑2182157

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: