Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,972,637 2 bp→GC coding (430‑431/546 nt) wecD → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,972,6370CG99.6% 772.6 / ‑5.0 235R144A (CGT→GCT) wecDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (108/126);  minor base T (0/1);  total (108/127)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01
*NC_0009133,972,6380GC99.6% 771.4 / ‑4.7 236R144A (CGT→GCT) wecDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (108/127);  minor base A (0/1);  total (108/128)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 242 total aligned reads displayed.

AATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAGAAAGCA  >  NC_000913/3972553‑3972730
                                                                                    ||                                                                                            
 aTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTa                                                                                      >  1:5237367/1‑93 (MQ=255)
  tGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTaa                                                                                     >  2:5352873/1‑93 (MQ=255)
    cGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTaaaa                                                                                   >  1:4311891/1‑93 (MQ=255)
     gACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaac                                                                                  >  1:1543994/1‑93 (MQ=255)
      aCAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaaca                                                                                 >  1:3671828/1‑93 (MQ=255)
       cAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaacaa                                                                                <  2:1015244/93‑1 (MQ=255)
         gATGCGCTAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACt                                                                              <  1:2714382/93‑1 (MQ=255)
         gATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACt                                                                              <  2:5731659/93‑1 (MQ=255)
         gATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACt                                                                              >  1:4082329/1‑93 (MQ=255)
         gATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACt                                                                              >  2:2982416/1‑93 (MQ=255)
         gATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACt                                                                              <  1:1599652/93‑1 (MQ=255)
          aTGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACtt                                                                             <  2:5477873/93‑1 (MQ=255)
          aTGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACtt                                                                             <  1:2856628/93‑1 (MQ=255)
          aTGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACtt                                                                             <  2:1730231/93‑1 (MQ=255)
          aTGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACtt                                                                             <  2:2034890/93‑1 (MQ=255)
           tGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACttt                                                                            >  1:620563/1‑93 (MQ=255)
           tGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACttt                                                                            >  2:2369670/1‑93 (MQ=255)
             cgcgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTgg                                                                          >  2:6358752/1‑92 (MQ=255)
                  aTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTg                                                                     <  2:5302087/93‑1 (MQ=255)
                   ttGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTgg                                                                    >  2:5972411/1‑93 (MQ=255)
                   ttGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTgg                                                                    <  2:6247971/93‑1 (MQ=255)
                     ggCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAACCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCg                                                                  <  1:5965163/93‑1 (MQ=255)
                       ccTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGAc                                                                <  2:3736024/93‑1 (MQ=255)
                         tgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGAccc                                                              >  2:2554107/1‑93 (MQ=255)
                            tggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCaaaa                                                           >  1:1909072/1‑93 (MQ=255)
                             ggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAAt                                                          <  2:6632623/93‑1 (MQ=255)
                             ggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAAt                                                          >  2:3826886/1‑93 (MQ=255)
                             ggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAAt                                                          <  1:2045813/93‑1 (MQ=255)
                              gctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATg                                                         <  2:5218516/93‑1 (MQ=255)
                              gctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATg                                                         <  1:6048015/93‑1 (MQ=255)
                              gctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATg                                                         <  1:5514339/93‑1 (MQ=255)
                              gctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATg                                                         <  1:1329238/93‑1 (MQ=255)
                               ctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATgg                                                        <  1:3477100/93‑1 (MQ=255)
                                tggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGGGGGGGGCGCCCCAAATggg                                                       >  2:2790234/1‑93 (MQ=255)
                                 ggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGc                                                      >  2:3494779/1‑93 (MQ=255)
                                   aCGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCaa                                                    <  1:6533522/93‑1 (MQ=255)
                                    cgcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAac                                                   >  1:421509/1‑93 (MQ=255)
                                       ggTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAcc                                                <  2:4061561/93‑1 (MQ=255)
                                       ggTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAcc                                                >  2:2503752/1‑93 (MQ=255)
                                        gTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCCACAcca                                               >  1:2982393/1‑92 (MQ=255)
                                           cAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgccg                                            <  2:1535412/93‑1 (MQ=255)
                                            aGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACcgccgc                                           >  1:108532/1‑93 (MQ=255)
                                            aGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACcgccgc                                           >  1:3817890/1‑93 (MQ=255)
                                            aGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACcgccgc                                           <  2:4100237/93‑1 (MQ=255)
                                            aGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACcgccgc                                           >  1:4974399/1‑93 (MQ=255)
                                            aGGCGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACcgccgc                                           >  2:6076767/1‑93 (MQ=255)
                                               tGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCt                                        <  1:5624907/93‑1 (MQ=255)
                                               tGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCt                                        <  2:4932806/93‑1 (MQ=255)
                                               tGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCt                                        <  2:5850778/93‑1 (MQ=255)
                                                gctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCtt                                       <  2:2089761/93‑1 (MQ=255)
                                                gctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCtt                                       <  2:2163835/93‑1 (MQ=255)
                                                gctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCtt                                       >  1:3167362/1‑93 (MQ=255)
                                                 ctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTa                                      <  2:1451741/93‑1 (MQ=255)
                                                   gagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaaa                                    <  1:4459687/93‑1 (MQ=255)
                                                   gagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaaa                                    <  2:4801434/93‑1 (MQ=255)
                                                    agctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAAc                                   <  2:1932095/93‑1 (MQ=255)
                                                     gctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACg                                  >  2:5250399/1‑93 (MQ=255)
                                                       tgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGat                                <  2:5114308/93‑1 (MQ=255)
                                                        gaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGata                               >  1:6584282/1‑93 (MQ=255)
                                                         aTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac                              <  1:1630402/93‑1 (MQ=255)
                                                            cAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacata                           >  2:3693394/1‑93 (MQ=255)
                                                            cAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacata                           <  2:1171444/93‑1 (MQ=255)
                                                            cAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacata                           <  1:4688526/93‑1 (MQ=255)
                                                             aaaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacatac                          >  1:72674/1‑93 (MQ=255)
                                                              aaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                         >  2:3859056/1‑93 (MQ=255)
                                                              aaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                         >  2:4313368/1‑93 (MQ=255)
                                                              aaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                         >  1:503275/1‑93 (MQ=255)
                                                              aaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                         >  2:2568782/1‑93 (MQ=255)
                                                              aaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                         <  2:1873208/93‑1 (MQ=255)
                                                              aaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                         >  2:1419219/1‑93 (MQ=255)
                                                              aaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                         >  1:839044/1‑93 (MQ=255)
                                                               aCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaa                        >  1:2473287/1‑93 (MQ=255)
                                                               aCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaa                        >  1:2413825/1‑93 (MQ=255)
                                                               aCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaa                        >  2:149130/1‑93 (MQ=255)
                                                               aCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaa                        >  1:3962875/1‑93 (MQ=255)
                                                                cGGCGCTAAACTGGGCGTATGCTCACGGTAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAg                      <  2:4855183/93‑1 (MQ=255)
                                                                  gcgcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgt                     <  1:1965432/93‑1 (MQ=255)
                                                                  gcgcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgt                     >  2:676640/1‑93 (MQ=255)
                                                                   cgcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtg                    <  1:6226813/93‑1 (MQ=255)
                                                                   cgcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtg                    >  1:5403764/1‑93 (MQ=255)
                                                                   cgcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtg                    >  2:6528290/1‑93 (MQ=255)
                                                                    gcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtgg                   >  1:661455/1‑93 (MQ=255)
                                                                    gcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtgg                   >  1:1703066/1‑93 (MQ=255)
                                                                     cTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggt                  <  2:621440/93‑1 (MQ=255)
                                                                      tAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggtg                 >  1:3100711/1‑93 (MQ=255)
                                                                       aaaCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGc                >  2:3604274/1‑93 (MQ=255)
                                                                       aaaCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGc                >  2:3301220/1‑93 (MQ=255)
                                                                        aaCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCg               >  1:2227226/1‑93 (MQ=255)
                                                                         aCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGa              >  1:4934937/1‑93 (MQ=255)
                                                                          cTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGaa             <  2:1367211/93‑1 (MQ=255)
                                                                          cTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGaa             <  2:447157/93‑1 (MQ=255)
                                                                          cTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGaa             >  1:1507543/1‑93 (MQ=255)
                                                                           tGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAAt            >  2:5064996/1‑93 (MQ=255)
                                                                           tGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAAt            <  1:2258436/93‑1 (MQ=255)
                                                                            gggCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATg           >  1:4020137/1‑93 (MQ=255)
                                                                            gggCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATg           >  1:1754126/1‑93 (MQ=255)
                                                                            gggCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATg           <  2:782267/93‑1 (MQ=255)
                                                                             ggCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGt          <  1:5752359/93‑1 (MQ=255)
                                                                              gCGTATTCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACTCCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTa         <  1:3214690/93‑1 (MQ=255)
                                                                              gCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTa         >  2:4591248/1‑93 (MQ=255)
                                                                                    ||                                                                                            
AATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAGAAAGCA  >  NC_000913/3972553‑3972730

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