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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,132,7870AC96.9% 185.4 / ‑5.4 65I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base C (38/24);  minor base T (2/0);  total (41/24)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132646‑2132925
                                                                                                                                             |                                                                                                                                           
cGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGGGACGAGGTGATCTGCATGTCGCCGTGATTGATAGGCAGGAt                                                                                                                                            >  2:344512/1‑143 (MQ=255)
 cgATCTGGTAGGGGGTAAACAGCTGATATTTTTTGCCGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTc                                                                                                                                           <  1:384814/142‑1 (MQ=255)
   aTCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTg                                                                                                                                         >  1:188187/1‑143 (MQ=255)
       ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAAGTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAGTTTTTTCGGGGACGAGTGTATCTGCATTTGACCGTGATTGATATGCAGGCGCTGCggc                                                                                                                                     >  2:489190/1‑141 (MQ=255)
            agGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctc                                                                                                                                >  2:915909/1‑143 (MQ=255)
            agGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCATTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGCTGACGAGATGATCTGCATTGCACCGTGATTGATATGCAGGCTCAGAGTCGCctc                                                                                                                                >  2:184023/1‑143 (MQ=255)
             gggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGCTTGATATGCACGCTCTGCGTCGCctct                                                                                                                               >  1:439761/1‑143 (MQ=255)
              ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                              >  2:786902/1‑143 (MQ=255)
              ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTATTCGTGCTTGCGGTAAAAGTGGAGGTACCCAGAGAATTTTTTCGGTGACGAGGTGAACTGCATTTCACCGTGAGTGATATGCAGGCTCTGAGTCGCctctg                                                                                                                              >  2:459141/1‑142 (MQ=255)
                 aagcAGCGGACATTGTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGGAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTTATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                           <  2:78407/140‑1 (MQ=255)
                  aaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTCTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTAc                                                                                                                          <  1:402004/143‑1 (MQ=255)
                        tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGGCGCCTCTTCTCCTTTCCa                                                                                                                    >  1:924689/1‑143 (MQ=255)
                          atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGATGTACCCGGAGAATTTTTTCGGTGACGAGGTGATTTGCATTTCACCGTGATTGATCTGCATGTTCTGCGTGGCCTCTTCTACCTTcctcc                                                                                                                  >  2:661946/1‑140 (MQ=255)
                          atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACg                                                                                                                  >  1:16597/1‑143 (MQ=255)
                          atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCCGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTGTCCACg                                                                                                                  >  1:996216/1‑143 (MQ=255)
                                    ctgcggcgGTCGAATTTGTCTTTGTGCTTGCGGTGAGAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                                                                                       <  2:633218/138‑1 (MQ=255)
                                    cTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                                                                                       <  1:1144814/143‑1 (MQ=255)
                                       gcgcGGTCGAACTTGTCGTTGTGCGTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                                                                                    <  2:264342/143‑1 (MQ=255)
                                         gcgGTCGAATCTGTGTTTGTTCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcca                                                                                                  <  2:1129681/143‑1 (MQ=255)
                                          cgGGCGAAGTTGTCTTGGTGCTTGCGGTAAAAGTGGAAGTGCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccac                                                                                                 <  1:467703/143‑1 (MQ=255)
                                            gTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGATTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacca                                                                                               <  1:856970/143‑1 (MQ=255)
                                               gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGCCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTGCTcccccacac                                                                                            >  2:794917/1‑141 (MQ=255)
                                                 aTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGACGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATc                                                                                          >  1:214380/1‑143 (MQ=255)
                                                     gTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATGTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACCCCACCACCATCCGCa                                                                                      >  2:244609/1‑143 (MQ=255)
                                                        tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACAGTGATTGATATCAAGGCTCTGCGTCTCCTCCTCTACTTTCCACCGTTCGCACGTGGCCACCACCATcacaacga                                                                                   >  1:1075028/1‑136 (MQ=255)
                                                         ttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                  <  1:383837/143‑1 (MQ=255)
                                                         ttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATACGCAGCTTCGGCAACTGCTCTTCTACTTTTCACGGTTCGCCGTACTCCACCACCATCAGCAGGaa                                                                                  >  1:1138383/1‑143 (MQ=255)
                                                             gCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCTCCGTACTCAACCAGCATCCGCAGGACGATg                                                                              >  1:370280/1‑143 (MQ=255)
                                                              cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGCAAGATGt                                                                             >  2:1020168/1‑143 (MQ=255)
                                                               ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTGTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                            <  2:268218/143‑1 (MQ=255)
                                                                  cGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCGGGAAGATGTCGTa                                                                         <  1:915909/143‑1 (MQ=255)
                                                                      aaaaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCTGCAGGAAGATGTCGTAATcc                                                                     >  1:1069102/1‑143 (MQ=255)
                                                                       ggtggggAAGTCCCTGGAGAACTTTGTCGGTGACGAGGTGATCTGCATTTCACCGTGGTTGATATGCAGGCTCCGCGTCGCCTCTTCGACTTGCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCt                                                                    <  2:58672/138‑1 (MQ=255)
                                                                         aGTGGAAGTACCCGGAGAATTTTTTTTTTTTTTTTGTGATCTGGATTTCACCGTACTTGATATGGAGGCTCGGCGTCGCCTCTTTTACGTTCTTCGGTTCGCCGTACTGCAACACCCTCTGCAGGCATTTTTCGTAATACtgt                                                                  >  2:104650/1‑143 (MQ=255)
                                                                          gTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCATCCGCAGGAAGATGTCGTAATCCtgtg                                                                 >  1:47374/1‑143 (MQ=255)
                                                                              aaGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCACCCGCAGGCCGATCTCGTACTCCTGTGCCGc                                                             >  1:34030/1‑143 (MQ=255)
                                                                                 tACCCGGAGAATTGTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                          <  1:682164/143‑1 (MQ=255)
                                                                                    ccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTACTCCTGTGCCGCTTTCAg                                                       >  2:444660/1‑143 (MQ=255)
                                                                                       gagaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCCGTTc                                                    >  1:861842/1‑143 (MQ=255)
                                                                                       gagaATTTTTTCGGTGACGAGGTCATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTc                                                    <  1:1060013/143‑1 (MQ=255)
                                                                                          attgctttCGGTGACGAGGTGGTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTGCCGCGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGt                                                 <  2:823470/138‑1 (MQ=255)
                                                                                            ttttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGCAt                                               >  1:201031/1‑143 (MQ=255)
                                                                                             tttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATc                                              >  2:795498/1‑143 (MQ=255)
                                                                                               tttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATGTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCAATCCGCCGCAAGCTGTAGTAACCCTGTCCCCCCTCCAGTTCGGTATCgc                                            >  2:988364/1‑142 (MQ=255)
                                                                                                    gTGGCGGGGTGGTGTGCATTTCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTGGCCGTACTCCACCCCCATCCGCAGGAAGATGTCGTAATCCTGGGCCGCTTTCAGTTCGGTATCGaacaaa                                       <  2:225048/143‑1 (MQ=255)
                                                                                                    gTGAGGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                       <  1:1056652/143‑1 (MQ=255)
                                                                                                    gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGGAGCGaacaaa                                       >  1:352326/1‑143 (MQ=255)
                                                                                                         gAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTCCCGCTTTCAGTTCGGTATCGACCAAACActc                                  >  1:570749/1‑143 (MQ=255)
                                                                                                               agcTGCAGTTCACTGGGATTGAGATGCAGGCTGTGCGTCGCCGCTTCTAGTTTCCACGGTTCGCCGTCCTCCCCCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGaa                            <  2:1096650/141‑1 (MQ=255)
                                                                                                               aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTACTCCTGTGCCGCTTTCAGTTCGGTACCGCACAAACACCCTTTGac                            >  2:590432/1‑142 (MQ=255)
                                                                                                                  tGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACCCTCTGaaagg                         >  2:673215/1‑141 (MQ=255)
                                                                                                                       ttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTTCACTTTCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTGTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATg                    <  1:654650/143‑1 (MQ=255)
                                                                                                                       ttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCCGTTCGGTATCGAACAAACACTCTCTGAACCGCCATg                    >  2:103518/1‑143 (MQ=255)
                                                                                                                       ttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATg                    >  1:407465/1‑143 (MQ=255)
                                                                                                                          aCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTCCTTTCCACGGTACGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGAATCGAACAAACACTACTTGAAACCCCATGccc                 >  2:627956/1‑143 (MQ=255)
                                                                                                                          aCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCGGGAAAATGTCGTAACACTGTACCGCTTTCAGCTCGGTATCAGCCAAACACTCTTTGAGACGCCATGccc                 >  2:48674/1‑143 (MQ=255)
                                                                                                                          aCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCGGTTGGGTATCGAACAAACACTCTTTGAAACGCCATGccc                 >  1:257669/1‑143 (MQ=255)
                                                                                                                            cGTGATGGGTATGCAGGCGCTGCGTCGCCGCTTCGACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg               <  2:463329/143‑1 (MQ=255)
                                                                                                                            cGTGATGGCTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg               <  2:947950/143‑1 (MQ=255)
                                                                                                                                 ttgttttGCAGGCCCTGCGTCGCCTCTTCTACTTTCCTCGGTTCGTCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaaa          <  2:330518/137‑1 (MQ=255)
                                                                                                                                  tgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCTCCACCGCCATCCGCAGGAAGATGTCCTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACCCTTTGACACGCCATGCGCAGGTaacg         >  2:525523/1‑141 (MQ=255)
                                                                                                                                     tatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAcc      >  2:798882/1‑143 (MQ=255)
                                                                                                                                       tGCAGGCTCTGCGTCCCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTg    <  1:711127/143‑1 (MQ=255)
                                                                                                                                        gCAGGCTCTGCATCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTACTGGAAACGCCATGCCCAGGTAAAGACCTgg   >  2:315334/1‑143 (MQ=255)
                                                                                                                                         cAGGCTCTGCGTCGCCTCTTCTACTTTCCCCGGTTCGCCGTACTCCCCCACCATCCGCTGGACGATGTCGTCATCCTGTGCCGCCTTCAGTTCGGTATCGAACCAACACTCTTTGATACGCCATGCCCAGGTAGAGAcatagt  >  2:421921/1‑138 (MQ=255)
                                                                                                                                             |                                                                                                                                           
CGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132646‑2132925

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: