Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,725,205:1 +T 100% intergenic (‑37/‑492) 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.T83.3% 168.0 / 32.0 66intergenic (‑37/‑492)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base . (2/9);  new base T (40/12);  total (45/21)
Fisher's exact test for biased strand distribution p-value = 4.36e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

TGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAAT  >  CP006881/3725068‑3725335
                                                                                                                                             |                                                                                                                                      
tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCGTTTTTGCAGGGAGccg                                                                                                                                        >  1:193373/1‑138 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAATAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCtt                                                                                                                                     >  2:349905/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCtt                                                                                                                                     >  3:400564/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCtt                                                                                                                                     >  7:199403/1‑139 (MQ=255)
   gCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGATAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCtttt                                                                                                                                   >  3:140908/1‑139 (MQ=255)
   gCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGTCAGGCAGCCtttg                                                                                                                                   >  3:146771/1‑138 (MQ=255)
         tAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGATCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTagcag                                                                                                                            <  1:4283/139‑1 (MQ=255)
            accgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTTGCAGGCAGCCTTTTTAGCAGccc                                                                                                                         >  1:96139/4‑139 (MQ=255)
             ccgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCGTTTTTGCAGCCCg                                                                                                                        >  1:28893/3‑139 (MQ=255)
              cgCGGCAGCCGCGGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACCCCCGTAAAGATACAAGCAAACGTTTGATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGt                                                                                                                       >  5:72402/2‑139 (MQ=255)
              cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGt                                                                                                                       >  1:127420/2‑139 (MQ=255)
              cgCGGCAGCCGCAGTTTTTATTTCCGGCGTCACCCAGTTCATCAAAAATACCATCCACACCCCTAAAGATAAAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCATTTTTAGCAGcccct                                                                                                                       >  5:124056/2‑137 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGtttt                                                                                                                    >  6:170300/1‑139 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGtttt                                                                                                                    >  3:200187/1‑139 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGttt                                                                                                                     >  8:66908/1‑138 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCATTTTGGCAGTCAGCCTGTTTAGCAGCCCGttt                                                                                                                     >  4:148249/1‑138 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTATCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGttt                                                                                                                     >  8:291731/1‑138 (MQ=255)
                                    tCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGcccgc                                                                                                   >  3:158151/1‑134 (MQ=255)
                                    tCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGcccgc                                                                                                   >  7:262785/1‑134 (MQ=255)
                                                cAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTTGCAGGCCGCCTTTTTAGCAGCCCGTTTTATAAGCCGCCGCCCGCAAGAAAAGACGGCCg                                                                                    >  5:61889/1‑139 (MQ=255)
                                                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                                 >  8:257671/1‑139 (MQ=255)
                                                       tcaAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCGTTTTATCAGCCCGGTTTCTCATCCGCAGCCCGCAAGAAAAGACGGCAGAgttctt                                                                             >  3:216001/1‑139 (MQ=255)
                                                        caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttagtt                                                                            >  5:75339/1‑135 (MQ=255)
                                                           aaaGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttttct                                                                         >  7:251906/1‑139 (MQ=255)
                                                               aCCATCAACACCCGAAACGATACAAGAAAACTTTTCATAGAAAAACATCCCTCCGCTTGTTTTTGGCAGGCAGCCTTTTTCGCAGCCGGTTTTCACAGCCGCCCCCCGCAACAAAAGACGGCCGATTTCTTTTCTCGCa                                                                     <  7:257671/139‑1 (MQ=255)
                                                                  aTCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCACGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCACCCCCCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAAc                                                                  >  4:143474/1‑139 (MQ=255)
                                                                     caccCCCGTAACGATACACGCAAACGTTTCATACAAAAACAGCCCGCCGCTTCTTTTGGGCCGGCCGCCTTTTTAGCCGCCCGTTTTCGCAGCCGCCGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAg                                                               <  5:170300/139‑1 (MQ=255)
                                                                       cacCCGTAAAGATACAAGCAAACGTGTCATAGAAAAACCTCCCTCCGCTTCGTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGCAACGACGGCCGATGTATTTTCGCGCAAACTAGCg                                                             >  4:353614/1‑139 (MQ=255)
                                                                            gTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGGAAACTAGCGGctct                                                        >  6:391037/1‑139 (MQ=255)
                                                                            gTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTAAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                         >  8:47837/1‑138 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                         >  1:215999/1‑137 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                         >  8:141297/1‑137 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                         >  1:54496/1‑137 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCCCCGCGTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGGTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGAGTTCATTTCTCGCAAACAAGCGGctc                                                         >  1:79504/1‑137 (MQ=255)
                                                                               aactaaCAAGAAATCGTTTCATAGAAAAACTTCCCTCGGTGTGTTTTTGCCCGGCCGCTTTTTATGCAGCCGGTTTGCTCAGCCCCCGCCCGAAAAAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAg                                                      <  2:107167/133‑1 (MQ=255)
                                                                                aGATACAACCAAACTTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCCGGCAGCTTTTTTAGCCGCCGGTTTTCTCCCCCGCCGCCGGAAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                    <  7:88928/139‑1 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGtt                                                                                                                      <  1:28423/70‑1 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGtt                                                                                                                      >  2:28423/1‑70 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGATTTCGTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                >  7:10708/1‑139 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                >  3:141809/1‑139 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCAAGCaaaa                                                 >  8:88928/1‑138 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAAATAGCGGCTCAAGCAAAAc                                                >  3:126785/1‑139 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAGCTAGCGGCGCGAGCACAAc                                                >  6:150847/1‑139 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACCTCCCAACGCGTCGTTTTGGCAGGCCCCCATTTTTTCAGCCCCGTTTATACGCCGCCGCCCCCCAGAAAAACCCGCCCAATTATTTTAACCCAAACGAGCCGCACTAGCcaaaaa                                               >  2:18345/1‑136 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACt                                               >  2:256650/1‑138 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTAGCAAACTAGCGGCTCTAGCAAAACt                                               >  3:147980/1‑138 (MQ=255)
                                                                                       cccAAACGTTTAATATAAAAACATCACTCCGGTTTTTTTTGGCCGCCCCTTTTTTTAGCCCCCCTTTTTCACCCCCCCCCCCCCAAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCa                                             <  4:275910/137‑1 (MQ=255)
                                                                                        gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCCGCCCGTTTTCGCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                            <  7:224776/139‑1 (MQ=255)
                                                                                        ccAACCGTTTCATACAAAAACACCCCTCCGCTTCTTTTGGCCCGGCAGCCTTTTTACCACCCCGTTTTCCCCGCCCCCGCCCCCAACAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                            <  8:199403/138‑1 (MQ=255)
                                                                                        ccAAACTTTTCATCTAAAAAAACCCCTCCGATTCTTTTTGCCAGGCCGCCTTTTTCCCCGCCCGTTTTCTCAGCCGCAGCCCCAAACAAACCACGGCCGATTTATTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                            <  3:360237/138‑1 (MQ=255)
                                                                                        ccAAACTTTTAATATAAAAACAACCCGCCGCTTATTTTTGCCGGGCCGCGTTTTACCCAGCCGGTTTTCTCCGCCCCCCCCCCCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                            <  4:141809/138‑1 (MQ=255)
                                                                                        ccAAACGTTTCATACAAAAACAGCCCGCCGCTTCTTTTTGGCCGGCAGCCTTTTTAGCAGCCCGTTTGCTCACCCGCCCCCCCCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                            <  7:256564/138‑1 (MQ=255)
                                                                                        ccAAACGTTTAATACAAAAACACCCCTCCGCTTTTTTTTGGCCGCCAGCGTTTTTAGCCCCCCGTTTTCTCCCCCGCCCCCCCAAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                            <  5:370053/138‑1 (MQ=255)
                                                                                        ccAAACGTTCCAAAGAAAAACATCCCTCCGCTTCTTTTTGGCTGGATGCTTTTTTAGCCGCCGGTTTTCTCAGCCGCCCCCCGAAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                            <  7:291731/138‑1 (MQ=255)
                                                                                          aaccGTTTCATATAAAAACATCCCTCCGCTTATTTTTGGCAGCCCGATTTTTTAGCACCCGTTTTTCTCACCCGCAGCCCTAAACAAACCACGGCCGATTTCTTTTGTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                          <  4:147980/136‑1 (MQ=255)
                                                                                           aaCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTAGCAAACTAGCGGCTCTAGAAAAACTCATTAc                                         >  1:134789/1‑139 (MQ=255)
                                                                                            aggTTTAATAAAAAAACATCCCTCGTTTTTTTTTTGGCCGGCAGCTTTTTTCGCACCCGTTTTACTCCCCCCCCGCCCCAAAAAAACCACGGCCGATTTCTTTTCTCGAAAACTAGCGGCTCTAGCAAAACTCATGAct                                        <  4:158151/137‑1 (MQ=255)
                                                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTATTTTCTCGCAAACTAGCGGCTCGAGCAAAAAGCATGAcg                                        >  3:11818/1‑138 (MQ=255)
                                                                                            aCGTTTCATAGAAAAACATCCCTCCGCTGCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCGCCGCCCCAGCCCGCAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                        <  7:66908/139‑1 (MQ=255)
                                                                                            aCGTTTCATACAAAAACATCCCTCCGCTTATTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCACACCCGCAGCCCGCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                        <  1:256650/139‑1 (MQ=255)
                                                                                            aCGTTTCATACAAAAACACCCCGCCGCTGCTTTTTGGCCGCCAGCCTTTTTCGCAGCCCGTTTTCTACCCCGCACCCCGCAATAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                        <  7:47837/139‑1 (MQ=255)
                                                                                            aCGTTTAATACAAAAACACCCCTCCGTTTGTTTTTGCCCGCACGCCTTTTTCGCAGCCGTTTTTCTCTGCCGCCGCCCGAAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                        <  8:251906/139‑1 (MQ=255)
                                                                                                 tCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCGAGCAAAACTCATGCCTCGAt                                    >  5:22612/1‑138 (MQ=255)
                                                                                                  cATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGCTTTATTTTATCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGa                                  >  5:329831/1‑139 (MQ=255)
                                                                                                   caccaaaaaCATCCCTCCGCTTCTTTTTGGCCGCCAGCCTTTTAAGCCCCCGTTTTTCACACCCCCACCCCGCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGAt                                 <  6:305581/135‑1 (MQ=255)
                                                                                                       aaaaaCATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGAc                             >  2:290193/1‑139 (MQ=255)
                                                                                                       aaaaaCATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGAc                             >  5:333614/1‑139 (MQ=255)
                                                                                                       aaaaaCATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAACACGGCCGCTTTCTTTTATCGCAAACAAGCGGCTCTAGCAAAACTCATGACGCGATGATAGAc                             >  7:315088/1‑139 (MQ=255)
                                                                                                         aaaCATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAACCTCATGACTGTATGATAGACCg                           >  8:8881/1‑139 (MQ=255)
                                                                                                                 cTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGaaa                     >  7:241909/1‑137 (MQ=255)
                                                                                                                 cTCCGCGTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTAAGCCGCAGCCCCCAAGAAAAGACGCCCGATTTATTTTATCGCAAAATAGCGGCTATAGGAAAACTAATTACTCTAAGATATACCGTGGAAag                    >  6:134157/1‑138 (MQ=255)
                                                                                                                                   ggcagCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTAGCAAACTAGCGGCTCTAGCAAAACTCATTACTATATGATTGCCACTGGAAAGAGAATTCTTTGTCACAat  >  4:290715/1‑138 (MQ=255)
                                                                                                                                             |                                                                                                                                      
TGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAAT  >  CP006881/3725068‑3725335

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: