Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,181,895 G→A H180H (CAC→CAT)  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,8950GA98.6% 211.3 / ‑4.6 69H180H (CAC→CAT) U712_10925Menaquinol‑cytochrome c reductase cytochrome b subunit
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (32/36);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

TAAAGCGGTCCTGAAATCCCCTGCTTCCGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAGAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCAGGGTGCCCCGCAAGGAGTGTTTTGACTTGCGTCCCGATCAAAGGCGTGGCTTCAGCGATCTGCAGCCCGACTTTTGTTGCAAACAGAGCCTTCATGTCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGCCC  >  CP006881/2181761‑2182023
                                                                                                                                      |                                                                                                                                
tAAAGCGGTCCTGAAATCCCCTGCTTCCGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGcc                                                                                                                              >  7:410791/1‑139 (MQ=255)
       gTCCTGAAATCCCCTGCTTCCGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAg                                                                                                                       <  7:398490/139‑1 (MQ=255)
       gTCCTGAAATCCCCTGCTTCCGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAg                                                                                                                       >  2:205034/1‑139 (MQ=255)
        tCCTGAAATCCCCTGCTTCCGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAgg                                                                                                                      >  4:172371/1‑139 (MQ=255)
        tCCTGAAATCCCCTGCTTCCGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAgg                                                                                                                      >  3:217985/1‑139 (MQ=255)
        tCCTGAAATCCCCTGCTTCCGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAg                                                                                                                       >  1:141754/1‑138 (MQ=255)
        tCCTGAAATCCCCTGCTTCCGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAg                                                                                                                       >  6:112663/1‑138 (MQ=255)
          cTGAAATCCCCTGCTTCCGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGTGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGa                                                                                                                     <  4:23306/138‑1 (MQ=255)
          cTGAAATCCCCTGCTTCCGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGa                                                                                                                     <  5:312359/138‑1 (MQ=255)
          cTGAAATCCCCTGCTTCCGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGa                                                                                                                     <  8:141596/138‑1 (MQ=255)
                tCCCCTGCTTCCGGATCATAATAAACTCTGCAGCCATCAGGCCGAACAATGAGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGgttt                                                                                                              <  4:321878/139‑1 (MQ=255)
                tCCCCTGCTTCCGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGtttt                                                                                                              <  7:53542/139‑1 (MQ=255)
                  cccTGCTTCCGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGttttaa                                                                                                            >  1:126024/1‑137 (MQ=255)
                  cccTGCTTCCGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGttttaa                                                                                                            >  6:56977/1‑137 (MQ=255)
                  cccTGCTTCCGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGttttaa                                                                                                            >  8:231773/1‑137 (MQ=255)
                  cccTGCTTCCGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGgtttaa                                                                                                            >  3:321878/1‑137 (MQ=255)
                       cttactgATCATAATAAAATGTGCAGCCATCAGCCCGAACAAAGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACtt                                                                                                         <  6:122855/131‑1 (MQ=255)
                           cGGATCATAATTAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGtt                                                                                                    <  7:395482/138‑2 (MQ=255)
                           cGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGttc                                                                                                   <  5:207488/139‑3 (MQ=255)
                           cGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGtt                                                                                                    <  3:403812/138‑2 (MQ=255)
                           cGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGtt                                                                                                    <  7:231773/138‑2 (MQ=255)
                            ggATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGa                                                                                                                     >  6:229491/1‑120 (MQ=255)
                            ggATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGa                                                                                                                     <  5:229491/120‑1 (MQ=255)
                                 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)
                                            tgCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGtt                                                                                                                >  3:40333/1‑109 (MQ=255)
                                            tgCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGtt                                                                                                                <  4:40333/109‑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                                                                               >  4:440753/1‑131 (MQ=255)
                                                      aGCCCGAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTcg                                                                               <  3:440753/132‑2 (MQ=255)
                                                           gAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAgt                                                                   >  1:416097/1‑138 (MQ=255)
                                                           gAACAATGCGGCAGGAAGAAAAAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAgt                                                                   >  7:128301/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              <  6:333233/139‑1 (MQ=255)
                                                                                                                tttGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt              <  2:416097/139‑1 (MQ=255)
                                                                                                                tttGGGCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTGTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt              <  1:341532/139‑1 (MQ=255)
                                                                                                                caggggCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt              <  1:405006/136‑1 (MQ=255)
                                                                                                                   gggCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt              >  4:160072/1‑136 (MQ=255)
                                                                                                                   gggCGCCTACAATATCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt              <  3:160072/136‑1 (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              <  3:283274/125‑1 (MQ=255)
                                                                                                                              atatCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt              >  4:283274/1‑125 (MQ=255)
                                                                                                                               tatCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt              <  7:287409/124‑1 (MQ=255)
                                                                                                                               tatCCGGATGCCCTGCAAGGAGCGTTTTAACTTGCGTTCCGATCAAAGGCGTGGCTTCGGCGATCTGCAGTCCGACTTTTGTTGCAAACAGAGCCTTCATATCCCAAGGAAGCAAGTAGCCTGt              >  8:287409/1‑124 (MQ=255)
                                                                                                                                      |                                                                                                                                
TAAAGCGGTCCTGAAATCCCCTGCTTCCGGATCATAATAAAATGTGCAGCCATCAGCCCGAACAATGCGGCAGGAAGAAAGAATACATGGATCGCAAAAAACCTAGTCAGCGTTTGGGCGCCTACAATATCAGGGTGCCCCGCAAGGAGTGTTTTGACTTGCGTCCCGATCAAAGGCGTGGCTTCAGCGATCTGCAGCCCGACTTTTGTTGCAAACAGAGCCTTCATGTCCCAAGGAAGCAAGTAGCCTGTAAAGCCAAGCCC  >  CP006881/2181761‑2182023

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: ‑