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,6370CG100.0% 1221.9 / NA 377R144A (CGT→GCT) wecDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (163/214);  total (163/214)
*NC_0009133,972,6380GC99.7% 1214.4 / NA 377R144A (CGT→GCT) wecDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (162/213);  minor base T (0/1);  total (163/214)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 383 total aligned reads displayed.

AATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAG  >  NC_000913/3972553‑3972724
                                                                                    ||                                                                                      
aaTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGTTCACGGt                                                                                 <  1:2809079/93‑1 (MQ=255)
    cGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTaaaa                                                                             >  2:4494388/1‑93 (MQ=255)
       cAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaacaa                                                                          <  2:9566160/93‑1 (MQ=255)
       cAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaacaa                                                                          >  2:3816470/1‑93 (MQ=255)
       cAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaacaa                                                                          >  2:6176753/1‑93 (MQ=255)
         gATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACt                                                                        >  2:5754841/1‑93 (MQ=255)
         gATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACt                                                                        >  1:5881282/1‑93 (MQ=255)
           tGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACttt                                                                      >  1:8144201/1‑93 (MQ=255)
           tGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACttt                                                                      <  2:9246100/93‑1 (MQ=255)
               cgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCgg                                                                  >  1:4085398/1‑93 (MQ=255)
                gAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCggg                                                                 <  1:2020143/93‑1 (MQ=255)
                 aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGt                                                                <  1:2307634/93‑1 (MQ=255)
                 aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGt                                                                <  1:8210149/93‑1 (MQ=255)
                 aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGt                                                                >  1:6759537/1‑93 (MQ=255)
                 aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGt                                                                <  2:3491068/93‑1 (MQ=255)
                 aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGt                                                                <  2:8601217/93‑1 (MQ=255)
                 aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGt                                                                <  2:6405935/93‑1 (MQ=255)
                  aTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTg                                                               <  2:1555212/93‑1 (MQ=255)
                  aTTGGCCTGATGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCATGGTAAAACAACTTTGCGGGTg                                                               <  2:2032175/93‑1 (MQ=255)
                   ttGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTgg                                                              >  1:4314791/1‑93 (MQ=255)
                   atGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTgg                                                              <  1:5168411/92‑1 (MQ=255)
                    tGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGc                                                             <  1:9458622/93‑1 (MQ=255)
                      gCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGa                                                           >  2:501718/1‑93 (MQ=255)
                         tgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGAccc                                                        >  1:1099395/1‑93 (MQ=255)
                         tgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGAccc                                                        >  2:7202754/1‑93 (MQ=255)
                          gctggctgGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCa                                                       >  2:1320859/1‑93 (MQ=255)
                           ctggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCaa                                                      <  1:6042506/93‑1 (MQ=255)
                            tggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCaaa                                                     >  1:7702242/1‑93 (MQ=255)
                             ggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCaaag                                                    >  1:5398729/1‑92 (MQ=255)
                              gctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATg                                                   <  2:7376704/93‑1 (MQ=255)
                               ctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATgg                                                  <  1:9983163/93‑1 (MQ=255)
                                tggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATggg                                                 <  2:4419372/93‑1 (MQ=255)
                                 ggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGc                                                <  1:5594564/93‑1 (MQ=255)
                                 ggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGc                                                <  2:7996418/93‑1 (MQ=255)
                                 ggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGc                                                <  2:2510034/93‑1 (MQ=255)
                                  gACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCa                                               <  2:8109669/93‑1 (MQ=255)
                                    cgcgGTGGAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTTAAGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAac                                             <  1:9629579/93‑1 (MQ=255)
                                    cgcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAac                                             >  1:460208/1‑93 (MQ=255)
                                     gcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAaca                                            <  2:7570696/93‑1 (MQ=255)
                                     gcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAaca                                            <  1:3713221/93‑1 (MQ=255)
                                      cgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAacac                                           <  2:7006067/93‑1 (MQ=255)
                                      cgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAacac                                           <  1:1015665/93‑1 (MQ=255)
                                       ggTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAcc                                          <  2:1774756/93‑1 (MQ=255)
                                        gTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccg                                         >  2:2824743/1‑93 (MQ=255)
                                        gTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccg                                         >  2:6771038/1‑93 (MQ=255)
                                         tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgc                                        >  1:9540395/1‑93 (MQ=255)
                                         tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgc                                        >  1:4873163/1‑93 (MQ=255)
                                         tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgc                                        <  1:6527386/93‑1 (MQ=255)
                                         tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgc                                        >  1:8769169/1‑93 (MQ=255)
                                           cAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgccg                                      >  1:6899561/1‑93 (MQ=255)
                                            aGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACcgccgc                                     >  2:2672946/1‑93 (MQ=255)
                                             ggTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCcgcg                                    >  2:2123270/1‑93 (MQ=255)
                                               tGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCt                                  >  2:6591719/1‑93 (MQ=255)
                                                 ctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTa                                >  1:3978788/1‑93 (MQ=255)
                                                  tgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaa                               >  2:7470479/1‑93 (MQ=255)
                                                  tgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaa                               <  1:4500574/93‑1 (MQ=255)
                                                   gagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCGAAATGGGCAACACGGCCGCGCTTaaa                              >  1:3118178/1‑93 (MQ=255)
                                                   gagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaaa                              <  1:5010076/93‑1 (MQ=255)
                                                    agctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAAc                             <  1:2249858/93‑1 (MQ=255)
                                                     gctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAAct                            >  1:7576833/1‑92 (MQ=255)
                                                     gctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACg                            >  2:4588472/1‑93 (MQ=255)
                                                     gctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGCCGACCCAAATGGGCAACACCGCCGCGCTTAAACg                            >  1:1612143/1‑93 (MQ=255)
                                                        gaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGata                         >  2:1885376/1‑93 (MQ=255)
                                                         aTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTCAACGatac                        >  2:414399/1‑93 (MQ=255)
                                                         aTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac                        <  1:8015430/93‑1 (MQ=255)
                                                         aTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac                        >  1:172197/1‑93 (MQ=255)
                                                         aTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac                        <  1:1692789/93‑1 (MQ=255)
                                                         aTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac                        <  1:7282756/93‑1 (MQ=255)
                                                         aTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac                        <  1:7920544/93‑1 (MQ=255)
                                                          tGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGataca                       >  2:9356646/1‑93 (MQ=255)
                                                          tGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGataca                       <  1:7090865/93‑1 (MQ=255)
                                                          tGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGataca                       <  1:914475/93‑1 (MQ=255)
                                                            cAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacata                     <  1:7857531/93‑1 (MQ=255)
                                                            cAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacata                     >  1:1398637/1‑93 (MQ=255)
                                                             aaaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacatac                    <  1:9269690/93‑1 (MQ=255)
                                                             aaaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacatac                    >  2:5457642/1‑93 (MQ=255)
                                                              aaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCACCACCGCCGCGCTTCACCGAtacataca                   >  1:6447446/1‑93 (MQ=255)
                                                              aaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                   >  1:6320514/1‑93 (MQ=255)
                                                                cGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaaa                 <  2:7344369/93‑1 (MQ=255)
                                                                 ggCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAg                <  1:8587671/93‑1 (MQ=255)
                                                                  gcgcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgt               >  1:9716134/1‑93 (MQ=255)
                                                                   cgcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtg              <  2:9855078/93‑1 (MQ=255)
                                                                    gcTAAACTGGGCGTATGCTCACGGTGAAACAACTTTGGGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGGTACATACAAAgcgg             >  2:5893560/1‑93 (MQ=255)
                                                                     cTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggt            <  2:3391234/93‑1 (MQ=255)
                                                                      tAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAGTGGGCAACACCGCCGTGCTTAAACGATACATACAAAgtggtg           <  1:5279465/93‑1 (MQ=255)
                                                                      tAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggtg           >  2:884585/1‑93 (MQ=255)
                                                                       aaaCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGc          >  1:8279559/1‑93 (MQ=255)
                                                                       aaaCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGc          >  2:2303704/1‑93 (MQ=255)
                                                                         aCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCGGGGCTTAAAGGAAACATACAAAGGGGGGCGa        >  1:9844606/1‑93 (MQ=255)
                                                                         aCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGa        >  1:3377044/1‑93 (MQ=255)
                                                                         aCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGAAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGa        >  1:1194782/1‑93 (MQ=255)
                                                                          cTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGaa       <  1:5529121/93‑1 (MQ=255)
                                                                          cTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGaa       <  2:9675617/93‑1 (MQ=255)
                                                                          cTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGaa       >  2:6258713/1‑93 (MQ=255)
                                                                           tGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAAt      <  1:9809958/93‑1 (MQ=255)
                                                                             ggCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGt    >  2:4724041/1‑93 (MQ=255)
                                                                             ggCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGt    <  1:7456934/93‑1 (MQ=255)
                                                                             ggCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCACTTAAACGATACATACAAAGTGGCGCGTAtct    >  2:8433519/1‑91 (MQ=255)
                                                                              gCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTa   <  1:5997104/93‑1 (MQ=255)
                                                                              gCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTa   <  2:92620/93‑1 (MQ=255)
                                                                                    ||                                                                                      
AATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAG  >  NC_000913/3972553‑3972724

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