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.5% 195.1 / NA 66I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (35/30);  minor base G (0/1);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132650‑2132925
                                                                                                                                         |                                                                                                                                          
tCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGACATGCCCGGATAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                                       <  2:205373/143‑1 (MQ=255)
    gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                                   >  2:1193075/1‑143 (MQ=255)
     tAGAGGGTAAACAGCTGATATTTTTTGCTGACGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAGGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCGCTGCGTCGc                                                                                                                                  >  1:1176564/1‑143 (MQ=255)
      agagGGTAAACAGCTGATATTTTTGGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGcc                                                                                                                                 <  2:220503/143‑1 (MQ=255)
          gggaaaCAGCTGATATTGTGTGCTCGCGCGGTCGCAGTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                             <  2:1122429/140‑1 (MQ=255)
            tagatatcTGATATTTTTGGCTGGGGCGGTCGGATTTGTTTTTGTGCTGTGGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttct                                                                                                                           <  1:993304/136‑1 (MQ=255)
                 tgCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAGACGTGGTCGTTCCCGGTGAATTTTTTCGGTGACGAGGTGATCTGCATTTCCCCGTGCTCTATATGCAGGCTCTGCGTCGCCTCTTCTACttt                                                                                                                      <  1:26492/142‑1 (MQ=255)
                  gCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTc                                                                                                                     >  2:631721/1‑143 (MQ=255)
                    tGATATTTTTTGCTGGCGCGGTCGACTTTGTCTTTGTGCTTGCGGTAAAGGTGGAAGTACCAGGAGCCTTTTATCGCGAACCGAGTGTTGTGGAGGTCATCGTGATTGATATGCCGGCCCCTCGTCGCCTCTTCTACTCTCCa                                                                                                                   >  1:1409943/1‑143 (MQ=255)
                    tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTCTCCa                                                                                                                   >  1:808829/1‑143 (MQ=255)
                     gATATTTTTTGCTGGCGCGGTCGGATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCAc                                                                                                                  >  2:1182241/1‑143 (MQ=255)
                     gATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTGGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCAc                                                                                                                  <  2:693676/143‑1 (MQ=255)
                         ttttttGCTGGCGCGGTTGAATTTGTCTTTGTTATTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGtt                                                                                                              <  1:1099027/143‑1 (MQ=255)
                              tGCGGGCGCGGTCGACGTGGTCGTTGTGCTTGCGGTAAAAGTGGAAGTTCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCCTTTCTCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                                                                                         <  1:171534/143‑1 (MQ=255)
                                cTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                                                                                       >  1:257404/1‑143 (MQ=255)
                                 tGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                                                                                      <  1:1049458/143‑1 (MQ=255)
                                   gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                                                                                    <  1:406894/143‑1 (MQ=255)
                                       ggTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacc                                                                                                <  2:598505/143‑1 (MQ=255)
                                        gTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACTCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGTCTCTGCGTCGCCTCTTCTCCTCTCCACGGTTAGccagactacccca                                                                                               >  1:1181110/1‑132 (MQ=255)
                                         tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccac                                                                                              >  2:1390356/1‑143 (MQ=255)
                                           gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCAGTGTGATGGAAAAGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacca                                                                                            >  2:1057105/1‑143 (MQ=255)
                                           gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCCCGGTTCGCCGTCCTccaccccca                                                                                            >  1:907368/1‑143 (MQ=255)
                                           gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccccca                                                                                            >  2:157384/1‑143 (MQ=255)
                                                   cTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAgg                                                                                    >  2:1066617/1‑143 (MQ=255)
                                                        tgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGCAGAt                                                                               >  1:452186/1‑143 (MQ=255)
                                                         gCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATg                                                                              <  2:1068036/143‑1 (MQ=255)
                                                           ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGCTCTGCATTTCACCGTGATTGATATCCAGGCTCTGCGTCGCCTCTTCTACTTTCAACGGTTCGCCGTACTCGACCACCATCCGCCggcaaatctc                                                                            >  1:154286/1‑135 (MQ=255)
                                                           ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                            <  2:121700/143‑1 (MQ=255)
                                                              cGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTa                                                                         >  2:1447431/1‑143 (MQ=255)
                                                               ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACCCCCCCACCATCCGCCGGAAGCTGTCGTaa                                                                        >  1:108507/1‑143 (MQ=255)
                                                                gtagaggTGGAAGTACCCGGAGAATTGTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAt                                                                       <  1:626921/137‑1 (MQ=255)
                                                                gTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCATCCGCCGGACGCTGTCGGAAt                                                                       >  2:1325805/1‑143 (MQ=255)
                                                                 tAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCAACGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATc                                                                      >  2:757065/1‑143 (MQ=255)
                                                                  aaaaGTGGAAGTACCCGGAGAATGTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATcc                                                                     <  2:1440927/143‑1 (MQ=255)
                                                                    aaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtg                                                                   <  2:334140/143‑1 (MQ=255)
                                                                     tgttgatgTACCCGGTCAATTTTTTCGGTGACGAGGTGATCTGCAGTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                  <  2:751655/136‑1 (MQ=255)
                                                                     aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGAACACGTCGTAATCCtgt                                                                  >  1:1076188/1‑143 (MQ=255)
                                                                     aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGACGTCGTAATCCtgt                                                                  >  1:1280360/1‑143 (MQ=255)
                                                                      gTGGAAGTCCCCGGAGAAGTGTTTCGGTGACGAGGTGATCTGCATTTCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                                 <  1:1131056/143‑1 (MQ=255)
                                                                        ggAAGTACCCGGAGAATTTTTTCGGTGAGGAGGTCCTCTGCAGTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCCTCTACTTTACACGGTCCGCCGTACTCCACCACCCTGCGCAGGAACACGTCCTAATGCTGTGcc                                                               >  1:106626/1‑143 (MQ=255)
                                                                         gtattcccgggAGAATTTTTTCGGTGAGGGGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCCCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCg                                                              <  1:1178703/134‑1 (MQ=255)
                                                                             tgcccGGAGAATTTGTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCTCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                          <  1:111246/141‑1 (MQ=255)
                                                                             tACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                          <  1:1193075/143‑1 (MQ=255)
                                                                             tACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                          <  2:842324/143‑1 (MQ=255)
                                                                                ccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTACACGTTTCTCCGTACTCCACCACCATCCGCAGGAAGCTGTCGTACCCCCATGCCGCTTTCAg                                                       >  1:506242/1‑143 (MQ=255)
                                                                                  ggAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCACGCACATGTCGTATTCCTGTGCCGcgtgcgttt                                                     >  2:1318185/1‑135 (MQ=255)
                                                                                        ttttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACGACCATCCGCCGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTAt                                               >  2:247504/1‑143 (MQ=255)
                                                                                          ttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTccgcaacc                                             >  1:481169/1‑136 (MQ=255)
                                                                                          tttaCGGTGACGGGCTGATCTGCATTTCACCATGATTGACATGCAGGCTCTGCGTCGCCTCTTCTACTTTCAAGGGTTCGCCGTACTCCGCCACCGTCCTCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCg                                             <  2:1268525/143‑1 (MQ=255)
                                                                                                ctgtcttggTGTTCTGCATTTCTCCGTGGTTGTTCTGCACGGTCTGCGTCGCCTCTTCCACTTTCCACGGTTCGCCGTACTCTACCACCCTCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                       <  1:68230/136‑1 (MQ=255)
                                                                                                        gTGATCTGTGTTTCCCCGTGATTGATATGCAGGCTCTCCGTCGCCTCGTCTACTTTCCTCGGTTCGCCGTACTCCACCCCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACTCTCttt                               <  1:1112947/143‑1 (MQ=255)
                                                                                                        gTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCATTGCCGCTTTCAGTTCGGTATCGGacacacacccctt                               >  1:832505/1‑138 (MQ=255)
                                                                                                          gATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGa                             <  2:49865/143‑1 (MQ=255)
                                                                                                           aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGaa                            <  2:920613/143‑1 (MQ=255)
                                                                                                           aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGaa                            >  2:1014007/1‑143 (MQ=255)
                                                                                                             cTGCATTTCACCGTGAGTGATATGCAGGCTCTGCGTCGACTCTTCTACTTTCCACGGTTAGCCGTACTCCACCAACAGCCGCAGGAAGGTGTCGTAATCGTGTGCAGCTTTCGGTTCGGTGTAGAACAAACACGCGTTGaaga                          >  2:364628/1‑141 (MQ=255)
                                                                                                                  tctacccGTGATTGATAGGCAGGCGCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCAt                     <  2:1433781/138‑1 (MQ=255)
                                                                                                                   ttcgctgTGATTTATATGCAGGCTCTGCGTCGCGTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATg                    <  1:737576/137‑1 (MQ=255)
                                                                                                                    tCACCGTGATTGATATGCAGGCCCTGCGTCGCATCCCCTACTTGCCGCGTCCCCCCGCCCTCCACCACCATCCGCACGACGCTGTCGTACTCCTGTGCCCCCTTCCGTCCGGGACCGACCAACCACTCTTTGAGACGCCaccc                   >  1:1033402/1‑140 (MQ=255)
                                                                                                                       ccGTGATTGGTTGGCAGGCCCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCa                <  1:631721/143‑1 (MQ=255)
                                                                                                                        cGTGATTGATATGCAGGCTCTGCGTCGTCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGGCCAg               >  1:384565/1‑143 (MQ=255)
                                                                                                                             tgtagatgcatgcTCTGCGCCGCCGCTCCTACTTTCCACGGTTCGCCGCATTCCTCCATCGTCCGCAGGAAGATGGCGTAATCCTGAGCCGCTCTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaaa          <  1:1252623/138‑1 (MQ=255)
                                                                                                                              tgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAAGGGTAGTAATCCTGTGCCGCCTTCAGTACGATCTCGAACCAACAACATTTGAAAGGCCAATCACAGGGAAAg         >  1:1176723/1‑143 (MQ=255)
                                                                                                                                atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTTCTCCACCTCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAc       <  2:1434947/143‑1 (MQ=255)
                                                                                                                                atatGCAGGCTCTGCGACACCTATTCGGCTATCCACGGTTAGCCGTACTCCACCACCATCCGCAGGAAGATGGCGTAATCCTGTGCCGCTTTCAGTTCGGCCTCGCACAAACCCTCTTTGAAACGCCATGCCCCGGTAAAGAc       >  1:883691/1‑143 (MQ=255)
                                                                                                                                  atGCAGGCTCTGCGTCGACTCTTATACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATACTGTGCCGCTTTCAGATCGGTATCGCACAGGCACGCGTTGAAACGCCAAGCTCAGGTAGACACCt     >  1:769465/1‑143 (MQ=255)
                                                                                                                                     cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGGCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCAGGCCCGGGTAAAGACCGGGt  >  2:38519/1‑143 (MQ=255)
                                                                                                                                         |                                                                                                                                          
TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132650‑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: