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,7870AC98.4% 174.4 / NA 64I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (33/29);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132649‑2132926
                                                                                                                                          |                                                                                                                                           
aTCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTg                                                                                                                                         >  1:1170060/1‑143 (MQ=255)
 tCTGGTAGAGGGTGAACAGCTTATATTTTTTGCTGGCGCGGTCGAATTTTTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                                        <  2:1362277/143‑1 (MQ=255)
  cTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCg                                                                                                                                       >  1:1110737/1‑143 (MQ=255)
     gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTGTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                                    <  1:570963/143‑1 (MQ=255)
     gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                                    <  1:552949/143‑1 (MQ=255)
      tAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                   >  2:901446/1‑143 (MQ=255)
      tAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                   >  1:164916/1‑143 (MQ=255)
      tAGAGGGTAAACAGCGGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                   <  2:1359762/143‑1 (MQ=255)
        gagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTAGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTAGCct                                                                                                                                 >  2:1221002/1‑143 (MQ=255)
          gggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCCCCGTGATTGATATGCCGGCTCTGCGTCGCctct                                                                                                                               >  2:1307401/1‑143 (MQ=255)
            gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGTGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtgttc                                                                                                                             >  1:914610/1‑143 (MQ=255)
                   gCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTc                                                                                                                      <  2:393174/143‑1 (MQ=255)
                    cTGATATCTTTTGCTGGCGCGGTCGACTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGCCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTcc                                                                                                                     <  1:1109082/143‑1 (MQ=255)
                     tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                    >  1:670066/1‑143 (MQ=255)
                       atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACg                                                                                                                  >  2:626537/1‑143 (MQ=255)
                       atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACg                                                                                                                  <  1:1193393/143‑1 (MQ=255)
                          ttttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGtt                                                                                                               <  1:371510/143‑1 (MQ=255)
                                  tGGCGCGGTCGAATTTGTCTTTGTACTTGCAGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCGTTTCACCGTGATTGATATGCAGACTCTGCGTCGCTTCTTCTAAGTTCCACGGTTCGCCGtcc                                                                                                       >  2:1037157/1‑141 (MQ=255)
                                    gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAACTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATCTCACCGGGATTGATATGCAGGCTCTGCGCCGCCTCTTCTACTGTCCACCGTTCGCCggcctc                                                                                                     >  2:299906/1‑138 (MQ=255)
                                     cgcgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCCCCTCTTCTCCTTTCCACGCTTCGCCGTACTcc                                                                                                    >  2:329426/1‑143 (MQ=255)
                                      gcgGTCTAATTTGTCTCTGTTCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTCTCGGTGCCGAGCTTTTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcca                                                                                                   <  1:51261/143‑1 (MQ=255)
                                      gcgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcca                                                                                                   >  1:959209/1‑143 (MQ=255)
                                      gcgGTCGAATCTGGCGTTGTGCTTGCGGTCAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcca                                                                                                   <  2:880872/143‑1 (MQ=255)
                                          tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccac                                                                                               >  1:1066190/1‑143 (MQ=255)
                                                  gTCTCTGTGCTTGCGGTAAAAGTGGAAGGACCTGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCa                                                                                       <  1:138409/143‑1 (MQ=255)
                                                    cttgctgCTGGCGGTAGAGGTGGTAGTACCCGGAGAGTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAgg                                                                                     <  2:424945/138‑1 (MQ=255)
                                                            ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCCCAAGACGAttcc                                                                             >  2:143894/1‑140 (MQ=255)
                                                            ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGCAGCCATTACCACCTCCCACGGTTCGCCTTACTCCACCACCATCCGCAGGCAGATGTc                                                                             >  1:625582/1‑143 (MQ=255)
                                                             tGCGTTAAACGTGGGAGTGCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTGCACGGTTCGCCGTAGTCCCCCACCATCCGCAGGAAGATGTCg                                                                            <  2:406519/143‑1 (MQ=255)
                                                             tGCGGTAAAAGTGGAAGTACTCGGAGACTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCg                                                                            <  1:845575/143‑1 (MQ=255)
                                                             tGCGGGAAAAGTGGGAGTACCCGGAGAATTGTGTCGGCGACGAGGTGATCTGCATTTCACCGTGATTGATGTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCg                                                                            <  1:116799/143‑1 (MQ=255)
                                                                ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTaa                                                                         <  1:1020891/143‑1 (MQ=255)
                                                                ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGAAATGCAGGCTCTGCGTCGCCTCTACAACTTTCAAATGTTATTTGTACCCCACCACCATCAGCAGGAAGATGTCGaaa                                                                         >  2:598681/1‑143 (MQ=255)
                                                                  tAAGAGTTGAAGTACCCGGCGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATc                                                                       <  2:738055/143‑1 (MQ=255)
                                                                    aaaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCTTCCTCTTGACGCTGTCGTCATCCt                                                                     >  1:917048/1‑143 (MQ=255)
                                                                         ggAAGTACCCGGAGAATTTGTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGCTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGcc                                                                <  1:1042473/143‑1 (MQ=255)
                                                                            aGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCCTCTACTGTCCACGGTTCGCCGTACTCCAACACCCGCCGCACGAAGATGCCGTAACCCTGTGCCGcc                                                             >  2:332411/1‑142 (MQ=255)
                                                                              tACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGTTTGACTGGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                           <  1:628074/143‑1 (MQ=255)
                                                                              gaCCCGGAGAAGTTTTTCGGTGACGAGGTGATATGCATTACACCGGGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                           <  1:249838/142‑1 (MQ=255)
                                                                                  cGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGt                                                       <  1:728004/143‑1 (MQ=255)
                                                                                    gagaAGTGTTGCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTc                                                     <  1:1276985/143‑1 (MQ=255)
                                                                                     agaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTGTACTTTCCACGGTTCGCCGTACTCCACCCACCTCCGCAGGAAGATGTCGTAATCCCGTGCCGCTGTCAGtcag                                                    >  2:717380/1‑140 (MQ=255)
                                                                                     agaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                    <  1:9751/143‑1 (MQ=255)
                                                                                      gaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCgg                                                   <  2:1110737/143‑1 (MQ=255)
                                                                                          tttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCAACCACCATCCGCAGGAAGATGTCGTAAACCTCTGCCCCTTCCACTTCggaaac                                               >  2:1338712/1‑139 (MQ=255)
                                                                                            tttCGGTGCCGGTGTGATCTGCTTTTCACCGTGATTGATATGCCGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCTTCCGCAGGAAGATGTAATATTCTGTTGCCGCTTTCAGTTCGGTATCGa                                             <  2:828289/143‑1 (MQ=255)
                                                                                                  tGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaac                                       >  1:492039/1‑143 (MQ=255)
                                                                                                   gACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCCGTATCGaacaaaca                                      >  1:1053709/1‑143 (MQ=255)
                                                                                                   gACCAGGTGAGCTGCACTGCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGCTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGATACCGaactgaaa                                      <  2:1089504/143‑6 (MQ=255)
                                                                                                    aCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAacac                                     >  1:946861/1‑143 (MQ=255)
                                                                                                           gATCTGCATTTCACCGTGATTGAGATGCAGGCCCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGCGATCGCCCACACACTCTTTGa                              >  2:714572/1‑143 (MQ=255)
                                                                                                              cTGCATTTCGCCGTGATTGATATGCAGGTTCTGCGTCGCCGCTTTTACTTTCCGCGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTTTCGAACAAACACTCTTTGAAAc                           <  2:1299129/143‑1 (MQ=255)
                                                                                                                gCAGTTCAGCGGGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTGCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGc                         <  2:1112311/143‑1 (MQ=255)
                                                                                                                 cATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACGCTGTGACACAcc                        >  1:318547/1‑143 (MQ=255)
                                                                                                                    ttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCAACATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAGCGCCATg                     >  1:462801/1‑143 (MQ=255)
                                                                                                                       aCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACCCTCTTTGAAACGCCATGccc                  >  1:1074103/1‑143 (MQ=255)
                                                                                                                         cGTGATTGATCTGCAGGCCCTGCGTCGCGCCTTCTACTCTCCGCAGTTCTCCGTACTCCACCTCCAGCCGCAGCCAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg                <  2:151355/143‑1 (MQ=255)
                                                                                                                         cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGCTGTCGTACGCCCGTGCCGCTTTCAGCACGATATCGAACAGACACCCGTTCAAACGCCATGCGCAg                >  2:1115740/1‑143 (MQ=255)
                                                                                                                             attgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaa            <  1:707754/143‑1 (MQ=255)
                                                                                                                              tcggtgtGCAGGCTCTGCGTCTCGTCTGCTACTTTCCATGTTTTGACGTTCTTCCCCACACTCCTCAGTAAGACGTAGTAATCCTGTGCTGCTTTCAGTTCGCTATCGAACAAACGCTCTTTGAAACGCCATGCCCAGGTaaa           <  2:996701/137‑1 (MQ=255)
                                                                                                                               tgatATGCAGGCTCTGCGACGCCTCTTCTACTTTCCACGGTTCGACGTACTCCACCACCATCCGCAGGAAGAGGTCGTAATCCTGTGCCGCTTTCAGTTCGCTATCGAACAAACACTCTACGAGACGCCACGCCCAGGTGAAg          >  2:1310091/1‑143 (MQ=255)
                                                                                                                                   aagCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTCTCGAGCAAACACTCTTTGAACCGCCATGCCCCGCTCAAGACCt      >  1:863622/3‑143 (MQ=255)
                                                                                                                                     gCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCCTTGGAACGCCATGCCCAGGTAAAGACCTgg    >  1:405492/1‑143 (MQ=255)
                                                                                                                                      cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTGGTAATCCCGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCGGGTCAGACCTgggt  >  2:163048/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                           
ATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132649‑2132926

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: