Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,724,542 A→G S95S (TCT→TCC)  U712_19200 ← SacPA operon antiterminator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,724,5420AG98.8% 298.5 / NA 83S95S (TCT→TCC) U712_19200SacPA operon antiterminator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (34/48);  minor base T (1/0);  total (35/48)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTTTACATGTCCGATCGCCCACAGGCCGATCTCGTATTCTTTTTTATAGAGCGCTTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCAGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCCCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCAATATGCTCCTCCGGCAGTGTTTGCAGGATTTGT  >  CP006881/3724411‑3724656
                                                                                                                                   |                                                                                                                  
tctTTTACATGTCCGATCGCCCACAGGCCGATCTCATATTCTTTTTTATAGAGCGCATTAATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTATATCGCAAAGGACAAATCGTCGGAAAgcg                                                                                                             >  5:130067/1‑139 (MQ=255)
  ttttACATGTCCGATCGCCCACAGGCCGATCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAgcg                                                                                                             <  8:35334/137‑1 (MQ=255)
  ttttACATGTCCGATCGCCCACAGGCCGATCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAgcg                                                                                                             <  4:267319/137‑1 (MQ=255)
  ttttACATGTCCGATCGCCCACAGGCCGATCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAgcg                                                                                                             <  1:188420/137‑1 (MQ=255)
  ttttACATGTCCGATCGCCCACAGGCCGATCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAgcg                                                                                                             <  5:79737/137‑1 (MQ=255)
  ttttACATGTCCGATCGCCCACAGGCCGATCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAgcg                                                                                                             <  7:186288/137‑1 (MQ=255)
  ttttACATGTCCGATCGCCCACAGGCCGATCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAgcg                                                                                                             <  5:203562/137‑1 (MQ=255)
     tACATGTCCGATCGCCCACAGGCCGATCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGat                                                                                                         <  4:4749/138‑1 (MQ=255)
      aCATGTCCGATCGCCCACAGGCCGATCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGat                                                                                                         >  4:21401/1‑137 (MQ=255)
      aCATGTCCGATCGCCCACAGGCCGATCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGat                                                                                                         <  3:21401/137‑1 (MQ=255)
          gTCCGATCGCCCACAGGCCGATCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAAt                                                                                                   >  5:192988/1‑139 (MQ=255)
                       caGGCCGATCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAgc                                                                                      >  8:281696/1‑139 (MQ=255)
                          gCCGATCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCTATTTATTCTGAACAAGCAGCCCATTTTGAATCCTTTCCATCGCAATGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAgcggc                                                                                   >  3:336833/1‑139 (MQ=255)
                              aTCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCATAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGcgccg                                                                               >  3:32222/1‑139 (MQ=255)
                              aTCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGcgccg                                                                               >  8:168099/1‑139 (MQ=255)
                              aTCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGcgccg                                                                               >  5:136221/1‑139 (MQ=255)
                              aTCTCATATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGcgccg                                                                               >  2:306433/1‑139 (MQ=255)
                                   atATTCTTTTTTCTAGAGCGCCTTTATCTCATGCCGCAATTTATGCTGCACAAGCAGCCCATTTAGAATCCTTTCGATCGCAACGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCGCAGCGGCGCCGCGAGc                                                                          <  7:194641/138‑1 (MQ=255)
                                   atATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTCAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAg                                                                           <  7:281696/137‑1 (MQ=255)
                                    tATTCTTTTTTATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCt                                                                         <  3:281845/139‑1 (MQ=255)
                                           ttttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAgcg                                                                                                             <  4:207468/96‑1 (MQ=255)
                                           ttttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAgcg                                                                                                             >  3:207468/1‑96 (MQ=255)
                                           ttttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGc                                                                          <  4:4827/131‑1 (MQ=255)
                                           ttttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGc                                                                          >  3:4827/1‑131 (MQ=255)
                                           ttttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTctcc                                                                     <  5:2287/136‑1 (MQ=255)
                                           ttttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTctcc                                                                     <  7:168099/136‑1 (MQ=255)
                                           ttttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTctcc                                                                     <  4:80699/136‑1 (MQ=255)
                                           ttttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTctcc                                                                     <  7:268429/136‑1 (MQ=255)
                                           ttttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTctcc                                                                     <  4:75883/136‑1 (MQ=255)
                                           ttttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTcgcagtc                                                                  >  1:62175/1‑133 (MQ=255)
                                              tataGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGACTCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCg                                                                <  4:231614/138‑1 (MQ=255)
                                                taGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCa                                                              <  1:280967/138‑1 (MQ=255)
                                                taGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCa                                                              <  4:113674/138‑1 (MQ=255)
                                                   agCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCTTTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGc                                                          <  7:112132/139‑1 (MQ=255)
                                                   agCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGc                                                          <  7:24929/139‑1 (MQ=255)
                                                   agCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGc                                                          <  8:327227/139‑1 (MQ=255)
                                                     cgcCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCt                                                         >  7:206326/1‑138 (MQ=255)
                                                        ctttATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCtgatg                                                     <  3:116325/138‑1 (MQ=255)
                                                        ctttATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCtgatg                                                     <  2:181381/138‑1 (MQ=255)
                                                                 aTGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATc                                                <  1:38534/135‑1 (MQ=255)
                                                                 aTGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCGGATGATATCCt                                              >  5:46004/1‑137 (MQ=255)
                                                                  tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAgcg                                                                                                             >  1:100146/1‑73 (MQ=255)
                                                                  tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAgcg                                                                                                             >  5:74594/1‑73 (MQ=255)
                                                                  tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAgcg                                                                                                             <  6:74594/73‑1 (MQ=255)
                                                                  tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAgcg                                                                                                             <  2:100146/73‑1 (MQ=255)
                                                                  tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGat                                                                                                         <  3:313637/77‑1 (MQ=255)
                                                                  tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTccg                                           >  1:79579/1‑137 (MQ=255)
                                                                  tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAg                                           >  4:338844/1‑139 (MQ=255)
                                                                  tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAg                                           >  5:107176/1‑139 (MQ=255)
                                                                  tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAg                                           >  1:254261/1‑139 (MQ=255)
                                                                  tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCGCTCCTTCCGCATAGCTGATGATATCCTCAg                                           >  1:284375/1‑139 (MQ=255)
                                                                  tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCGCTCCGTCCTCATATCTGAGGATATCCTCAg                                           >  7:49259/1‑139 (MQ=255)
                                                                  tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCACTCCTTCCGCAAAGCTGATGATATCCACAg                                           >  4:275524/1‑139 (MQ=255)
                                                                  tGCAGCAATTTATTCTGCACAAGCAGCACATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGat                                                                                                         >  4:313637/1‑77 (MQ=255)
                                                                  tGCAGCAATTTATTCTGCACAAGAAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGAGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGATCTGCTTCCGCATAGCGGATGATATCCTCAg                                           >  8:10954/1‑139 (MQ=255)
                                                                     agcaATTTATTCTGCACACGCAGCCCATTATGAATCATTTCGATCGCAAAGGACAAATGGTCGGAAAGCGGGTTATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTa                                        <  3:275524/139‑1 (MQ=255)
                                                                     agcaATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTa                                        <  3:224600/139‑1 (MQ=255)
                                                                     agcaATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTa                                        <  3:338844/139‑1 (MQ=255)
                                                                     agcaATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCACAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTa                                        <  8:285677/139‑1 (MQ=255)
                                                                              ttCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATTAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTatat                               <  7:10954/139‑1 (MQ=255)
                                                                              ttCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTatat                               <  1:306433/139‑1 (MQ=255)
                                                                                   cacaAGCAGCCCATTTTGAATCGTTTCGATCCCTAAGTACAAATGGTCGGAACGCGCGATATGAATTTGTTCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTTATGATATCCTCAGCTATTTCTATATGctc                           <  5:312579/138‑1 (MQ=255)
                                                                                   cacaAGCAGCCCATTTTGAATCCTTTCTATCGCAAAGGACAAATGGTCTGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGAGAGCTATACTTAAGCATAGCTGATGATATACTCAGCTAATTCAatattcaaa                          >  6:237389/1‑134 (MQ=255)
                                                                                   cacaAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGctcc                          >  7:283529/1‑139 (MQ=255)
                                                                                   cacaAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGctc                           <  1:83537/138‑1 (MQ=255)
                                                                                   cacaAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGc                             <  2:284375/136‑1 (MQ=255)
                                                                                   cacaAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTAGGAAAGCGCGATATGAAAGTCATCGCTCAGAGGCGCCGCGAGCAATCCATACGCATAGCTGATGATATACACAGCAATTTCTATATGatc                           >  6:312579/1‑138 (MQ=255)
                                                                                       agcagcCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCg                      <  6:136221/139‑1 (MQ=255)
                                                                                        gcagcCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCTTAGCTGATGATATCCTCAGCTAt                                       <  3:200262/121‑1 (MQ=255)
                                                                                        gcagcCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCTTAGCTGATGATATCCTCAGCTAt                                       >  4:200262/1‑121 (MQ=255)
                                                                                         cagcCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGc                    <  3:90795/139‑1 (MQ=255)
                                                                                          agcCCATTTTTAATCCTATCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCa                   <  4:173912/139‑1 (MQ=255)
                                                                                            cccATTTTGAATCCTTTCGATCGCAAAGTAAAAATGGTCGTAAAGAGATATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATCATATCCTCAGCTATTTCTATATGCTCCTCCGGCa                   <  5:237389/137‑1 (MQ=255)
                                                                                            cccATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAgt                 <  6:110431/139‑1 (MQ=255)
                                                                                             ccATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAg                                                                           <  2:117933/80‑1 (MQ=255)
                                                                                             ccATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAg                                                                           >  1:117933/1‑80 (MQ=255)
                                                                                                      aTCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGa       <  6:107176/139‑1 (MQ=255)
                                                                                                          tttCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCAGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTg   >  7:239171/1‑139 (MQ=255)
                                                                                                           ttCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCa                                                              >  2:72071/1‑79 (MQ=255)
                                                                                                           ttCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCa                                                              <  1:72071/79‑1 (MQ=255)
                                                                                                           ttCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATTATATACTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGt  >  5:281384/1‑139 (MQ=255)
                                                                                                           ttCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGt  >  1:175339/1‑139 (MQ=255)
                                                                                                           ttCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATAAGAATGTGATCGCTCAGCGGCGCCGCGCGCTCTCCTTCCCCATAGCTGATGATATACTCATCTATTTCTATATACTCCTCAGGAAGTGTTAGCACGATATGt  >  1:334631/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                  
TCTTTTACATGTCCGATCGCCCACAGGCCGATCTCGTATTCTTTTTTATAGAGCGCTTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCAGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCCCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCAATATGCTCCTCCGGCAGTGTTTGCAGGATTTGT  >  CP006881/3724411‑3724656

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