Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,132,787 A→C I204S (ATC→AGC)  wcaA ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,132,7870AC97.0% 183.5 / NA 66I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (29/35);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGG  >  NC_000913/2132647‑2132924
                                                                                                                                             |                                                                                                                                         
ggATGTGGTCGCGGGCAAACAGCTGAGATTTTTTGCGGGCGCGGGCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTc                                                                                                                                         <  2:1008762/143‑1 (MQ=255)
 gATCTGGTGGAGGGTAGGCGGCTGATATTTTTGGCTGGCGCGGTCGAATTTGTCTTTGTGCTGGCGGAAAAAGTGGAAGGACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGctct                                                                                                                                        <  2:194369/143‑1 (MQ=255)
   tCTGGTAGAGGGTAAACAGCTGATAGGTGTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                                      <  1:562111/143‑1 (MQ=255)
     tGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATGTGCAGGCTCTGCGt                                                                                                                                    >  2:748844/1‑143 (MQ=255)
      ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTc                                                                                                                                   >  1:1117297/1‑143 (MQ=255)
        tAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                 >  2:186565/1‑143 (MQ=255)
        tAGAGGGTAAACAGCTGATATTTTTTGCTGGCGAGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCAGTTCACCGCGATTGGTATGCAGCATCttgggcgc                                                                                                                                 >  2:125364/1‑136 (MQ=255)
         agagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGGATTTGTCGTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATGGTTTCGGTGACGAGGGGATCTGCATGTCACCGTGATTGATATGCACGATCTGCGTCGcc                                                                                                                                >  2:1040163/1‑143 (MQ=255)
         agagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGcc                                                                                                                                >  1:295528/1‑143 (MQ=255)
          gagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCct                                                                                                                               >  2:536350/1‑143 (MQ=255)
           acgctAAACAGCTGATCTTTTTTGCTGGCGCAGTCGAAGCTCTCCTTGTGCTTCCGGTAAAAGTGGCAGTACCCGGCGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctc                                                                                                                              <  1:315414/139‑1 (MQ=255)
             ggcaaaCAGCAGATGTTTGTTGCGGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                            <  2:455412/140‑1 (MQ=255)
                 aaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTAc                                                                                                                        <  1:944102/143‑1 (MQ=255)
                      cTGATATTTTTTGCTGGCGCGGGCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGGACACGGAGAATTTTTTCGGTGACGAGGTGATCTGCAGTTCACCGTGAGTGATATGCAGGCGCGGCGTCGCGTCTTGTAGTGTcc                                                                                                                   >  2:1009879/1‑143 (MQ=255)
                       tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTAGAGGAAGTGCAAGTGCAGGCCCTGCGTCGCCTCTTCTCCTTTccc                                                                                                                  >  2:1139527/1‑142 (MQ=255)
                       tGATATTGTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGGATCCGGAGAATGTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                  <  1:1169801/143‑1 (MQ=255)
                         atatTTTTTGCTGGTTCGTGCGAATTTGTCTTTGTGCTTGCGGTAAAAGGGGAAGTACCCGGAGGATTTTTTGGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGTCTGTGCGTCTCCTCTTCTACTGTCCACg                                                                                                                >  2:87889/1‑143 (MQ=255)
                         atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTGACCGTGATTGATACGCAGCCTCTGCGTGGCCCCTACTACTTTCCagg                                                                                                                >  1:1143359/1‑141 (MQ=255)
                          tatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTGTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTCCTTTCCACgg                                                                                                               >  1:381477/1‑143 (MQ=255)
                           atgttttGCTGGCGCGGTCGAGTTTGTCTTTGTGCTGGCGGTAAAGGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGt                                                                                                              <  2:781996/140‑1 (MQ=255)
                                ttGCTGGCGCGGTGGAATTGGGCTTTGTGCTTGCGGTAAAAGGGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGcc                                                                                                         <  2:678199/143‑1 (MQ=255)
                                ttGCTGGCGCGGTCGACTTTGTGTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGcc                                                                                                         <  2:973307/143‑1 (MQ=255)
                                ttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGcc                                                                                                         <  1:570048/143‑1 (MQ=255)
                                      gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                                                                                   <  1:36159/143‑1 (MQ=255)
                                      gcgcGGTCGAATTTGTCTGGGGGTGTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                                                                                   <  2:286253/143‑1 (MQ=255)
                                       cgcgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcc                                                                                                  <  2:1117297/143‑1 (MQ=255)
                                              gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCCCCGTAATccaccacca                                                                                           >  1:1154835/1‑143 (MQ=255)
                                              gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTCCGCCGTACTccaccgcaa                                                                                           >  1:174025/1‑139 (MQ=255)
                                              gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTACAACTTGCCACGGTTCGCCGTACTccaccacca                                                                                           >  1:1062140/1‑143 (MQ=255)
                                              gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGACTCTGCGTCGCCTCTTCTACTTTCCACGGTACGCCGTCCTccacaccca                                                                                           >  2:668129/1‑138 (MQ=255)
                                                 tttGTCTTTGTGCTTGCTGTATAAGTGGAAGTACCCGGAGAATGTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATcc                                                                                        <  1:948183/143‑1 (MQ=255)
                                                  ttGTCTTGGGTCTTGTGGTTGAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATCTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCg                                                                                       <  2:368948/143‑1 (MQ=255)
                                                      cTTTGTGCTTGCGGTTCAAGTGGAAGTACCCGGAGAATTTTTTCCGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGGAgg                                                                                   <  1:114342/143‑1 (MQ=255)
                                                             cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACCTTCCACCGTTCGCCGTACTCCACCACCACCCGCAGCAAGCTGt                                                                            >  2:668974/1‑143 (MQ=255)
                                                              ttGCGGTAAAAGTGGAAGTTCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                           <  1:1016214/143‑1 (MQ=255)
                                                              ttGCGGTAAAAGTGGAAGTGCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                           <  1:1066977/143‑1 (MQ=255)
                                                              ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGCTGTc                                                                           >  2:1077979/1‑143 (MQ=255)
                                                              ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                           <  1:901531/143‑1 (MQ=255)
                                                               tGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCg                                                                          <  2:736207/143‑1 (MQ=255)
                                                               cgCGGTAAAAGTGGAAGTACCCGGAGAATTTTTGCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCg                                                                          <  1:455211/142‑1 (MQ=255)
                                                                gccgtgaaaggtgAACTACCGGGAGACTTTTTTCGGTCACGATGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGt                                                                         <  2:837204/131‑1 (MQ=255)
                                                                  ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCGGCGTCGCCTCTTCTACTTTCCACGGATCACCGTACTCCACCCCCATCCGCAAGAAGATGTCGCaa                                                                       >  2:426732/1‑143 (MQ=255)
                                                                      aaaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTaagcat                                                                   >  1:587432/1‑139 (MQ=255)
                                                                        atgcgAAGTTCCCGGAGAACTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                 <  1:659849/139‑1 (MQ=255)
                                                                            gAAGTCCCCGGAGAATTTTTTCGGTGACGGGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTGCTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCg                                                             <  2:772942/143‑1 (MQ=255)
                                                                              aGTACCCGGAGAATGTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGTAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCt                                                           <  1:423840/143‑1 (MQ=255)
                                                                                tatccGGAGAATTTTTTCGGCGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                         <  1:587719/140‑1 (MQ=255)
                                                                                tACCCGGAGAATCTTTTCGGTGACGAGGCTATCTGCATTTCACCGTGATCTATACTCAGTCTCCTCGTCGCCTCCTCTACTTGCCACGGTTCGCCCAACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                         <  1:779295/143‑1 (MQ=255)
                                                                                 aCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTc                                                        >  1:1118420/1‑143 (MQ=255)
                                                                                      gagaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCCGCGTAGCCTCTTCTACTTTCCACGGTTCGCCGTACTGCACCACCATCCGCAGGAAGATGTCGTAAGCCAGCGCAGCTATCAgaac                                                   >  2:604532/1‑140 (MQ=255)
                                                                                       agaATTTTTTCGGTGTCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                  <  1:1056568/143‑1 (MQ=255)
                                                                                                     gACGGGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaaca                                    <  1:284178/143‑1 (MQ=255)
                                                                                                      aCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTACTCCTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTAGCGagcaaccag                                   >  1:717586/1‑142 (MQ=255)
                                                                                                        gAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTCCTTTCTACGGTTCGCCGTACTCCCCCACCATCCGCAGGAAAGTGTCCTAATACCGTGCCGCTCTCAGCTCGGAATCGactcaaagagc                                 >  2:692635/1‑133 (MQ=255)
                                                                                                        gAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGGACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCACTTCGATATCGAacacacacac                                 >  2:66020/1‑141 (MQ=255)
                                                                                                             gATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGa                            <  1:745908/143‑1 (MQ=255)
                                                                                                               tCTGCATTTCACCGTGATTGATGTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGaaa                          <  2:1037979/143‑1 (MQ=255)
                                                                                                                    aTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACCCTCTTTGAAACGCCa                     >  1:451544/1‑143 (MQ=255)
                                                                                                                       tCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGc                  >  1:1018754/1‑143 (MQ=255)
                                                                                                                         aCCGTGATTGCTCTGCAGGCTCGGCGTCTCCTCTTCTACTGTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGccc                <  2:879103/143‑1 (MQ=255)
                                                                                                                            gTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACACCACCACCATACACAGGCAGCAGTCGTAATCCTGTGCGCCATTCAGTTCAGTATCGAACAAACACACTTTGAAACGACATGCACAgg             >  2:209560/1‑143 (MQ=255)
                                                                                                                                 tgatATGCAGGCGCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAg        <  1:186565/143‑1 (MQ=255)
                                                                                                                                  gatATGCACGCTCTGCGTCGCCTCGTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGa       <  1:17968/143‑1 (MQ=255)
                                                                                                                                     tatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGGAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTGTCGAAGAAACACGTTTTGAAACGCCTTGCCCCGGTACAGAcc     >  1:569811/1‑143 (MQ=255)
                                                                                                                                      atGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCt    >  1:659340/1‑143 (MQ=255)
                                                                                                                                       tGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTg   <  1:95623/143‑1 (MQ=255)
                                                                                                                                        gCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCAGCGGCAGGAAAATGTCGTAATCCTGTGCCGCTTTCCGCTCGGCAGCGAGCAAGAGCTCTGCGAAACACCATGCCCAGGAAAacacagcg  >  2:171501/1‑139 (MQ=255)
                                                                                                                                             |                                                                                                                                         
GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGG  >  NC_000913/2132647‑2132924

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