Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,725,209 +T intergenic (‑41/‑488) U712_19205 ← / → U712_19210 Uncharacterized protein ywcI/Minor extracellular protease vpr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,725,2051.T81.5% 149.2 / 28.2 56intergenic (‑37/‑492)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base . (2/8);  new base T (31/13);  total (35/21)
Fisher's exact test for biased strand distribution p-value = 8.78e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATAT  >  CP006881/3725068‑3725337
                                                                                                                                            |                                                                                                                                         
tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCCCCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGAAGGCAGCCt                                                                                                                                          >  7:49583/1‑139 (MQ=255)
tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCt                                                                                                                                          >  7:49577/1‑139 (MQ=255)
tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCt                                                                                                                                          >  6:377912/1‑139 (MQ=255)
tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCt                                                                                                                                          <  6:324599/139‑1 (MQ=255)
tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCt                                                                                                                                          >  5:324599/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATATAAAAACATCCCTCCGCGTCTTTTTGGCAGGCCGCCtt                                                                                                                                        >  4:383286/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTGCGTTTTTGCAGGGAGCCtt                                                                                                                                        >  5:252234/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCGTTTTGGCAGGCAGCCtt                                                                                                                                        >  4:357733/1‑139 (MQ=255)
  aGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGCCCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCGTTTTGGCAGGCAGCCgtt                                                                                                                                       >  1:340555/1‑139 (MQ=255)
       aaTAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCGTTTTagc                                                                                                                                >  7:323905/1‑139 (MQ=255)
         tAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTagcag                                                                                                                              <  8:224879/139‑1 (MQ=255)
         tAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTagcag                                                                                                                              <  3:120999/139‑1 (MQ=255)
            accgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGccc                                                                                                                           >  8:68126/4‑139 (MQ=255)
            accgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTGTTTTGCAGGCAGCCGGTTTAGCAGccc                                                                                                                           >  4:325229/4‑139 (MQ=255)
              cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGt                                                                                                                         >  3:23586/2‑139 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGttt                                                                                                                       >  7:272392/1‑138 (MQ=255)
                               ttCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGcc                                                                                                         <  5:107579/138‑1 (MQ=255)
                                         cGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGaaaa                                                                                               >  6:28077/1‑138 (MQ=255)
                                                  gTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAt                                                                                    >  4:63699/1‑139 (MQ=255)
                                                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTGTTTTTGGCAGGCCGCCTGTTTAGCAGCCCGTTTTAGCAGCCGCAGCCCGCAAGAAAAGACGGCCGCtt                                                                                   >  2:383935/1‑139 (MQ=255)
                                                   ttCATCAAAAAGACAATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCGTTTTGGCCGGCAGCCTATTTAGCAGCCCGTATTCTCCCCCGCAGCCCCCAATAAAAGACGGcctcag                                                                                   >  4:112230/1‑135 (MQ=255)
                                                        caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCATCCCCTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttt                                                                              >  3:340290/1‑139 (MQ=255)
                                                           acacaCCATCCACACGCTCAAATATACAAGCAAAGGTTTAAAAGAAAAACATCCCTCCGTTTCTTTTTGGCAGGCAGCATTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCCAAATAAAACACGGCCGAtttcttttct                                                                           <  6:252234/135‑1 (MQ=255)
                                                           aaaGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttttct                                                                           >  7:216906/1‑139 (MQ=255)
                                                           aaaGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcgttta                                                                            >  5:347805/1‑137 (MQ=255)
                                                               aCCATCCACACCCGTAAAGATACACCAACACTTTTCTTTGAAAAACATCCGTCCTCTTTTTTTTGGAAGGCAGCCTTTTACGCGGCGCGTTTGCTCAGCCCCCGCCCCCACGAAAACTCGGCCGATTTCTTTTCTCGCa                                                                       <  8:49583/139‑1 (MQ=255)
                                                               aCCATCCACACCCGAAACGATACAACCAACGGTTTCTTACAAAAACACCCCTCCGCTTCTTTTGGGCAGCCAGCCTTTTTCCCCCCCCGTTTTCTCCCCCCCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCa                                                                       <  8:49577/139‑1 (MQ=255)
                                                                     cacaCCCTTAAATATAAAAGCAATGGTTTTATAAAAAAACATCCGTCGGTTTTTTTTTGGCGGGCTGGTTTTTTAGCCGCCGGTTTTCTCCCCCGCCGCGGCAATAAAAACACGGCCGATTTCTTTTCTCGCAAACTAg                                                                 <  4:340290/139‑1 (MQ=255)
                                                                     acagcccGAAAATAAACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTGTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCCCCGCCCGCAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAg                                                                 <  3:89543/135‑1 (MQ=255)
                                                                            gTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTATCAGCCGCAGCCCGCAAGAAAAGACGGCCTATTTATTTTGTCGCAAACTAGCggagat                                                          >  5:244016/1‑135 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAACTAGCGGctc                                                           >  4:378060/1‑137 (MQ=255)
                                                                               aaGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCATCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCGTTTCTCGCAAACTAGCGGCTCTat                                                        >  3:345310/1‑137 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACGTCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTGTGGCAAACAACCGGCGATAGCcaaa                                                   >  1:203756/1‑138 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCactgt                                                  >  4:89543/1‑135 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa                                                   >  3:400520/1‑138 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACAAGCGCCTCGAGCaaaa                                                   >  6:132349/1‑138 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTATCGCAAACAAGCCGCTCTAGCaaaa                                                   >  4:230857/1‑138 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGTTTTCTCCGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCGTTTGTCCGAAACTAGCGGCGAGAGCAAAAca                                                 >  1:122685/1‑137 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCCAGAAAAGACGGCCGATTTCGTTTCTCGCAAACGAGCGGCTCTAGCAAAACTc                                                >  8:47171/1‑139 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGc                                                                <  3:403585/123‑1 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGc                                                                >  4:403585/1‑123 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTACCAAACTACCGGCTCAAGCAAAACt                                                 >  2:41973/1‑138 (MQ=255)
                                                                                        gaaaaCGTATCAAACAAAACCCTCCCTCCGTTTATTTCGGCCACGCAGCCTTTCTATCCGCCGGTTCTCTCAGCCGCAGCCCCCATCAACAGACGCCCGATTTCTTTTCTCGCAAACAAGCGGCTCTAGCAAACCTCAt                                              <  8:323905/137‑1 (MQ=255)
                                                                                        gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCGCTAGCAAAACTCAt                                              >  6:399406/1‑139 (MQ=255)
                                                                                        gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCCCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                              <  7:68126/139‑1 (MQ=255)
                                                                                        gCAAACGTTTCATAGAAAAACATCCCTCCGCTTATTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGCCTCTAGCAAAACTCAt                                              >  5:172020/1‑139 (MQ=255)
                                                                                          aaaCTTTCCAAATAAAAACATCCCTCCGCTTCTTTTTGGCAGCCAGCCTTTTAAGCAGCCCGTTTTCTCAGCCGCCGCCCGAAACAAAAGACGGCCGATTTATTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                            <  8:216906/139‑1 (MQ=255)
                                                                                          aaaCGTTTCTTATAAAAACATCCCTCCGCTTCTTTTTGGCTGGCAGCCTTTTTAGCAGCCCTTTTTCTCCGCCGCAGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                            <  6:347805/139‑1 (MQ=255)
                                                                                            ccGTTTCATACAAAAACATCCCTCCGCTTTTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCGCAGCCGCAGCCCGAAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                          <  1:41973/138‑1 (MQ=255)
                                                                                            ccGTCTCATTGAAAAACAACCCTCGGCTGCTTTTTGGCCGGCAGCCTTTTTTGCCGCCCTTTTGCACCCCCGCCCCCCGAAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                          <  4:77369/138‑1 (MQ=255)
                                                                                            acttttCATACAAAAACACCCCTCCGCTTCTTTTTGGCTGGCAGCATTTTTAGCAGCCCTTTTCCTCACCCGCCCCCCCAAATAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                          <  5:377912/136‑1 (MQ=255)
                                                                                            aCGTTTCATATAAAAACACCCCTCCGTGTATTTTTGGCAGGCTGCCTTTTATGCCGCCCGTTTGCTCCGCCGCTCCCCGCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                          <  4:93119/139‑1 (MQ=255)
                                                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGCTTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAcg                                          >  5:115785/1‑138 (MQ=255)
                                                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCGCTAGCAAAACGCATGAct                                          >  7:368816/1‑139 (MQ=255)
                                                                                            aCGTTTATTACAAAAACACCCCCCCGCTTGTTTTTGCCCGCCCGCCTTTTTTGCCCCCGGTTTTCCCCGCCCCACCCCCAAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                          <  6:82942/139‑1 (MQ=255)
                                                                                             cGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGATTTCTCAGCCGCAGCCCGCAAGAAAAAACGGCCCAGTACTATAAACGCAAACTAACCCCTCTAGCAAAACTAAAGAct                                          >  8:393387/1‑138 (MQ=255)
                                                                                             cGTGTCATACAAAAACACCCCCCCGCGTCTTTTTGGCAGGCCCCCTTTTTCCCCGCCCGGTTGCACCCCCGCCGCCCGCAATAAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                          <  2:35259/139‑1 (MQ=255)
                                                                                                 tCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTtct                                                                                                                    <  2:10132/61‑1 (MQ=255)
                                                                                                 tCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTtct                                                                                                                    >  1:10132/1‑61 (MQ=255)
                                                                                                        aaaaCACCCCCCGGCTTCTTTTGGCCCGGCCGGTTTTTAAGCCCCCCGTTTACACCGCCCCCCCCCCCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGAcc                              <  1:121942/139‑1 (MQ=255)
                                                                                                                 cTCCGCTTCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGTTTTCTAAGCCGCAGCCCGCCAGAAAAGACGGCCGAGTTGTTTTGTCGCAAACATGCGGCGCGAGCAAAACTCATTCATCTATGAGATACCGGGGaaa                       >  6:262288/1‑137 (MQ=255)
                                                                                                                 cTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGAGTTCGTTTCTCGCCAACTAGCGGCTCTAGCAAAACTCAGGACTCTATGATAGACCGTGGaaa                       >  7:165254/1‑137 (MQ=255)
                                                                                                                 cTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGaaa                       >  1:358022/1‑137 (MQ=255)
                                                                                                                 cTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTAATTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTATATGATAGACCGTGGaaa                       >  3:318368/1‑137 (MQ=255)
                                                                                                                                    gcagCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAatat  >  3:57848/1‑139 (MQ=255)
                                                                                                                                            |                                                                                                                                         
TGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATAT  >  CP006881/3725068‑3725337

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