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.T79.1% 158.1 / 43.8 67intergenic (‑37/‑492)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base . (4/10);  new base T (36/15);  total (42/25)
Fisher's exact test for biased strand distribution p-value = 6.07e-03
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.

TGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCA  >  CP006881/3725068‑3725331
                                                                                                                                           |                                                                                                                                  
tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGgccgccg                                                                                                                                    >  1:339178/1‑138 (MQ=255)
tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCt                                                                                                                                   >  2:1625/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCCCACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCCGGCAGCCtt                                                                                                                                 >  1:140874/1‑139 (MQ=255)
 gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTGTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGgccgccgt                                                                                                                                   >  6:365729/1‑137 (MQ=255)
       aaTAGACCGCGGCAGCCGCTGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGCCCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATACCTCCGCTTCGTTTTGGCAGGCAGCCtttttag                                                                                                                           >  3:27401/1‑135 (MQ=255)
       aaTAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCCCACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCtttttag                                                                                                                           >  3:222290/1‑135 (MQ=255)
         tAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTagcag                                                                                                                       <  5:133854/139‑1 (MQ=255)
            accgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGccc                                                                                                                    >  7:75303/4‑139 (MQ=255)
            accgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGccc                                                                                                                    >  1:131621/4‑139 (MQ=255)
              cgCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGt                                                                                                                  >  7:320780/2‑139 (MQ=255)
                cGGCAGCCGCAGGTTTTCTGTCCGGCGTCACCCAGTTCATCAAACAGCCCCGCCACACCCGGAAAGATACAAGCAAACGTGTCATAGAAACACCTCCCTCCGCGGCGTTTTGGCAGGCAGCCTGTTTAGCAGCCCGgtt                                                                                                                >  7:123502/1‑139 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGttt                                                                                                                >  4:202184/1‑138 (MQ=255)
                 ggCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTGTTTGGCAGGCAGCCTGTTTAGCAGCCCGtttt                                                                                                               >  7:168846/1‑139 (MQ=255)
                   cAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTtct                                                                                                             <  1:162536/139‑1 (MQ=255)
                           gTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCg                                                                                                       >  5:112546/1‑137 (MQ=255)
                                 cTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGcccg                                                                                                <  5:310645/138‑2 (MQ=255)
                                      cGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGcccgcaaga                                                                                           >  2:266343/1‑132 (MQ=255)
                                         cGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTTGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCCCAAGaaaa                                                                                        >  5:122742/1‑138 (MQ=255)
                                              cccAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTATAAGCCGCCGCCCGCAACAAAAGACGGc                                                                                  >  2:80637/1‑139 (MQ=255)
                                                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                             >  4:378708/1‑139 (MQ=255)
                                                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCAGCAAGAAAAGACGGCCGaat                                                                             >  8:429851/1‑137 (MQ=255)
                                                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCGTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                             >  1:29366/1‑139 (MQ=255)
                                                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGCCGTTTTATCCGCCCGGTTTCTCAGCCGCACCCCCCAAGAAAAGAAGGCCgcgg                                                                             >  1:160090/1‑136 (MQ=255)
                                                        caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCCGGCAGCCTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttttc                                                                      >  3:100804/1‑139 (MQ=255)
                                                        caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttt                                                                        >  8:126288/1‑139 (MQ=255)
                                                         aaaaaGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttctttt                                                                       >  4:382248/1‑139 (MQ=255)
                                                             aGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTATTTTTGGCAGGCAGCCTTTTTAACAGCCCGTTTTCTCATCCCCACCCCGCAACAAAAGACCGCCGATTTATTTTCACg                                                                   >  7:160644/1‑139 (MQ=255)
                                                                         cccGAAAAGAAACAAGCAACCGTTTCAAAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGc                                                       <  1:408660/139‑1 (MQ=255)
                                                                         agcttAAAGAAACAAGCAATGGTTTCATAGAAAAACATCCCAGCGATTATTTTTGGCAGGCAGCTTTTTATGCCCCCCTTTTTCCCCGCCGCCGCCCGCACGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGc                                                       <  8:123947/135‑1 (MQ=255)
                                                                          ccGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGccc                                                                                                 <  3:367072/96‑1 (MQ=255)
                                                                          ccGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGccc                                                                                                 >  4:367072/1‑96 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCTGGCCGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGGCGAGTTCTTTTGTCGGAAACTAGGGGctc                                                     >  1:279595/1‑137 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                     >  7:136224/1‑137 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTAGCAAACTAGCggggc                                                     >  7:397391/1‑134 (MQ=255)
                                                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                     >  8:242218/1‑137 (MQ=255)
                                                                                acaTACAAGCAAACTTTTCATACAAAAACATCCCTCCGCTTCTTTTTGGCCGGCTGCCTTTTTAGCAGCCCGTTTTCTCCGCCGCCGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                <  7:429851/137‑1 (MQ=255)
                                                                                acaTACAACCAACCTTTGCATAAAAAAACACCCCACCGATTATTTTGGCCCGGCTCCATTTTAAGCAGCCCGTTTCCTCCCCCGCACCCCCAAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                <  5:140466/137‑1 (MQ=255)
                                                                                aGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                <  7:126288/139‑1 (MQ=255)
                                                                                aGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                <  1:80637/139‑1 (MQ=255)
                                                                                  aTACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCGCTAGCaa                                               >  1:49152/1‑138 (MQ=255)
                                                                                   tACAAGCAAACGTTTCATAGAAAAAAATCCCTCCCCTTCTTTTTGGCAGGCAGCCTTTTTAGAAGCCCGTTTTCTCAGCAGCAGCCCGCAAGAAAAGACGGCCGATTTCTTATCTAGCAAACAATCGGCTCTAGCaac                                              >  7:22241/1‑137 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                            >  7:345503/1‑139 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                            >  3:124657/1‑139 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                            >  3:299563/1‑139 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                            >  3:118354/1‑139 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTAGCAAACTAGCCGCTCAAGCaaaa                                             >  8:282559/1‑138 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTAATTTATAGCAAACAAGCGGCTCTAGCaaaa                                             >  8:175919/1‑138 (MQ=255)
                                                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCGTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                            >  6:120380/1‑139 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTGTTTTTTGCAGGCATCGTTTTTAGCAGCCCGTTTTATAAGCCGCAGCCCGCAAGAAAAGACGGACGATTTCTTGTGTGTCAAACATGCGGTTCGAACaaaaat                                           >  6:353776/1‑136 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCa                                                                                                     >  6:331303/1‑80 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCa                                                                                                     <  5:331303/80‑1 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTAGCAAACTAGCGGCTCAAGCCAAACt                                           >  4:103100/1‑138 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTGTAGCCAACTATCGGCGCTAGCAAAAca                                           >  1:230933/1‑137 (MQ=255)
                                                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCCCAGCCCGCAAGAAAATACGGCCGAGTTCTTTTCTCGCCAACTAGCGGCTCTAGCCAAAcct                                          >  7:123947/1‑137 (MQ=255)
                                                                                        gCAACCGTTTCATACAAAAAAAGCCCCCCGCTTCTTTTTGGCCGCCACCCTTTTTAGCAGCCCGTTTTCTCAGCCCCAGCCCCAAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                        <  3:398258/139‑1 (MQ=255)
                                                                                        gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTCGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                        <  1:1625/139‑1 (MQ=255)
                                                                                        gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                        <  8:320780/139‑1 (MQ=255)
                                                                                          aaaCTTTTCATATAAAAACATCCCTCCGCTTCTTTTTGCCAGGCACCATTTTTCCCAGCCCTTTTTCACCGCCGCAGCCCGAAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                      <  8:345503/139‑1 (MQ=255)
                                                                                           aaCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCGAGACAAACTCATGa                                      >  6:216563/1‑138 (MQ=255)
                                                                                            ccttttCATCCAAAAACATCCCGCCGCTTCTTTTTGGCAGGACGCCTTTTACGCACCCCGTTTTCGCAGCCCCACCCCCAAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                    <  7:197823/136‑1 (MQ=255)
                                                                                            aCGTTTCATATAAAAAGATCCCGCCGCTGCTTTTTGGCAGGCAGCGTTTTTAGCCGCCCGTTTTCTCAGCCCCAGCCCGCAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                    <  1:266343/139‑1 (MQ=255)
                                                                                            aCGTTTCATATAAAAACATCCCTCCGTTTCTTTTTGGCAGGCAGCTTTTTTTCCCGCCGGTTTTCTCAGCCGCAGCGCGAAAAAAATGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                    <  7:282559/139‑1 (MQ=255)
                                                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                    <  3:202184/139‑1 (MQ=255)
                                                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                    >  3:109979/1‑139 (MQ=255)
                                                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTGTTTTCTCGCCAACTAGCGGCTCTAGCAAAACTCATGAct                                    >  2:90411/1‑139 (MQ=255)
                                                                                             cGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                    >  5:347465/1‑138 (MQ=255)
                                                                                                          aacaacCCTCCGGTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCTTTTTCGCCGCCCCCCCCCCAAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGt                      <  8:136224/139‑1 (MQ=255)
                                                                                                                    cGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGaa            >  5:398337/1‑139 (MQ=255)
                                                                                                                         cttTTTGGCAGGCCCCCTTTATACCACCCCGTTTACTAACCCGCCCCCCGCAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCt        <  1:23235/138‑1 (MQ=255)
                                                                                                                             ttgcctggccgCCTTTTACGCCGCCGGTTTTCGCCGCCCCACCCCCCAACAAACCACGGCCGATTTCTTTTCTTGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGt    <  3:378708/132‑1 (MQ=255)
                                                                                                                             ttGCCAGGCCGCCTTTTAATCAGCCCGTTTTCTCCGCGTCCCCCCGAAATAACCCACGGCCGATTTCTTTTCTCGAAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGt    <  4:108444/138‑1 (MQ=255)
                                                                                                                             tgggcagccaGCGTTTTTAGCCGCCCGTTTTCTACGCCCCACCCCCCAACAAAACCCGCCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGt    <  8:345811/136‑1 (MQ=255)
                                                                                                                               ggcaggcaGCATTTTTAGCAGCCCGTTTTATCCGCCGCACCCCGAAATAAAAGACGGCCTATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTca  <  4:222290/138‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                  
TGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCA  >  CP006881/3725068‑3725331

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: