Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,019,162 A→G Y39H (TAC→CAC)  ompA ← outer membrane protein 3a

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,019,1620AG99.0% 316.7 / ‑4.7 100Y39H (TAC→CAC) ompAouter membrane protein 3a
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (64/35);  total (65/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 103 total aligned reads displayed.

ACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTG  >  NC_000913/1019072‑1019251
                                                                                          |                                                                                         
aCGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCt                                                                                         <  1:1325894/93‑1 (MQ=255)
  gggTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTgg                                                                                       >  1:799776/1‑93 (MQ=255)
    gTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGa                                                                                     >  1:4393983/1‑93 (MQ=255)
    gTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGa                                                                                     >  2:1057120/1‑93 (MQ=255)
     ttAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGAc                                                                                    <  1:3899150/93‑1 (MQ=255)
         ccTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGTGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGc                                                                                >  2:2664997/1‑93 (MQ=255)
          cTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGcc                                                                               >  2:635906/1‑93 (MQ=255)
          cTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGACGGGCCATTGTTGGTGATGAAACCAGTGTCATGGTGCTGGGACCAGcc                                                                               >  2:3747939/1‑93 (MQ=255)
           tGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGccc                                                                              >  1:3800548/1‑93 (MQ=255)
            ggTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGcccg                                                                             >  1:2899392/1‑92 (MQ=255)
            ggTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCa                                                                             >  1:2899393/1‑93 (MQ=255)
             gTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAg                                                                            <  1:2555372/93‑1 (MQ=255)
              tAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGt                                                                           >  1:479381/1‑93 (MQ=255)
              tAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGt                                                                           <  2:73911/93‑1 (MQ=255)
              tAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAAACAGTGTCATGGTGCTGGGACCAGCCCAGt                                                                           >  1:479382/1‑93 (MQ=255)
               aaCCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGtt                                                                          >  2:2207706/1‑93 (MQ=255)
               aaCCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGtt                                                                          >  1:1476752/1‑93 (MQ=255)
                  caccaAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAg                                                                       >  2:376413/1‑93 (MQ=255)
                  caccaAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAg                                                                       >  1:1623517/1‑93 (MQ=255)
                   accaAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGc                                                                      >  2:1621059/1‑93 (MQ=255)
                    ccaAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCa                                                                     >  2:3337284/1‑93 (MQ=255)
                     caAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCAc                                                                    >  1:2320214/1‑93 (MQ=255)
                     caAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCAc                                                                    >  1:931450/1‑93 (MQ=255)
                      aaaaGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCAcc                                                                   >  1:2210745/1‑93 (MQ=255)
                      aaaaGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCAcc                                                                   <  1:4239781/93‑1 (MQ=255)
                      aaaaGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGGGCTGGGACCATCCCAGTTTAGCAcc                                                                   >  2:2472043/1‑93 (MQ=255)
                       aaaGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCa                                                                  <  2:3440346/93‑1 (MQ=255)
                       aaaGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCa                                                                  >  2:2497932/1‑93 (MQ=255)
                            aCCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTa                                                             >  2:728594/1‑93 (MQ=255)
                            aCCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTa                                                             >  1:3109098/1‑93 (MQ=255)
                                gcgcCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAg                                                         >  1:2775061/1‑93 (MQ=255)
                                 cgcCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAgg                                                        <  1:4592330/93‑1 (MQ=255)
                                  gcCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAggg                                                       >  2:3367738/1‑92 (MQ=255)
                                  gcCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGgt                                                       >  1:718920/1‑93 (MQ=255)
                                  gcCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGgt                                                       >  2:1130993/1‑93 (MQ=255)
                                   cccAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGgtg                                                      >  2:3546452/1‑93 (MQ=255)
                                     cAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGtt                                                    <  2:211186/93‑1 (MQ=255)
                                        ttGGTTTTCATGGGTCGGGCCATTGGTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCAACAGTGTAACAGGTGTTAtt                                                 >  2:175163/1‑92 (MQ=255)
                                         tGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCt                                                >  1:766820/1‑93 (MQ=255)
                                          ggTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCtt                                               >  1:352611/1‑93 (MQ=255)
                                          ggTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCtt                                               <  1:4410364/93‑1 (MQ=255)
                                          ggTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCtt                                               >  1:2696813/1‑93 (MQ=255)
                                          ggTTTTCATGGGTCGGGCAATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCtt                                               <  1:3451267/93‑1 (MQ=255)
                                           gTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCttt                                              >  1:2099904/1‑93 (MQ=255)
                                           gTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCttt                                              >  1:117640/1‑93 (MQ=255)
                                           gTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCttt                                              >  2:2389897/1‑93 (MQ=255)
                                            ttttCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTc                                             >  1:100201/1‑93 (MQ=255)
                                            ttttCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTc                                             >  1:100200/1‑93 (MQ=255)
                                            ttttCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTc                                             <  2:4725385/93‑1 (MQ=255)
                                             tttCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCg                                            <  1:2261403/93‑1 (MQ=255)
                                               tCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGa                                          <  2:188065/93‑1 (MQ=255)
                                               tCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGa                                          <  2:2238327/93‑1 (MQ=255)
                                                   gggTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCgg                                      >  1:392577/1‑93 (MQ=255)
                                                    ggTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGc                                     <  1:3174873/93‑1 (MQ=255)
                                                    ggTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGc                                     <  1:3174872/93‑1 (MQ=255)
                                                     gTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGcc                                    <  1:1584314/93‑1 (MQ=255)
                                                     gTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGcc                                    <  1:2052852/93‑1 (MQ=255)
                                                         ggCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTgcg                                >  1:793355/1‑93 (MQ=255)
                                                          gCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTgcgc                               <  2:3683031/93‑1 (MQ=255)
                                                          gCCATTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTgcgc                               >  1:1320157/1‑93 (MQ=255)
                                                            cATTGTTTTTGATGAAACCAGTGTCATGGTGCTGGGACCAGGCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTa                             <  1:4099626/93‑1 (MQ=255)
                                                             aTTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTag                            >  2:1311368/1‑92 (MQ=255)
                                                             aTTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTAc                            >  1:4315975/1‑93 (MQ=255)
                                                             aTTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTAc                            >  2:1311366/1‑93 (MQ=255)
                                                             aTTGTTGTTGATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTAc                            >  1:2652742/1‑93 (MQ=255)
                                                                gttgttgATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGt                         <  2:2234159/93‑1 (MQ=255)
                                                                gttgttgATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGt                         <  2:2234160/93‑1 (MQ=255)
                                                                gttgttgATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGt                         <  1:3151071/93‑1 (MQ=255)
                                                                 ttgttgATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTa                        >  2:887565/1‑93 (MQ=255)
                                                                 ttgttgATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTa                        >  1:1704544/1‑93 (MQ=255)
                                                                 ttgttgATGAAACCAGTGCCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTa                        <  1:991453/93‑1 (MQ=255)
                                                                 tggttgATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTa                        >  2:3292618/3‑93 (MQ=255)
                                                                  tgttgATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAg                       >  1:3930649/1‑93 (MQ=255)
                                                                   gttgATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGc                      <  2:1498403/93‑1 (MQ=255)
                                                                   gttgATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGc                      >  1:2095009/1‑93 (MQ=255)
                                                                    ttgATGAAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCg                     >  1:1876360/1‑93 (MQ=255)
                                                                     tgatgaAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGa                    <  2:286032/93‑1 (MQ=255)
                                                                        tgaAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAc                 <  1:1858706/93‑1 (MQ=255)
                                                                        tgaAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAc                 >  1:3700881/1‑93 (MQ=255)
                                                                        tgaAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAc                 >  2:2215557/1‑93 (MQ=255)
                                                                        tgaAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAc                 >  1:2481643/1‑93 (MQ=255)
                                                                         gaAACCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAcc                >  1:3435241/1‑93 (MQ=255)
                                                                          aaaCCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAcca               <  2:2539692/93‑1 (MQ=255)
                                                                           aaCCAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAccag              >  2:4191392/1‑93 (MQ=255)
                                                                             ccAGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAccagcc            >  1:907225/1‑93 (MQ=255)
                                                                               aGTGTCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAccagccag          >  1:4011257/1‑93 (MQ=255)
                                                                                gtgtCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAccagccagg         >  1:1396377/1‑92 (MQ=255)
                                                                                gtgtCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGt         >  1:2585418/1‑93 (MQ=255)
                                                                                gtgtCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGt         <  1:1409037/93‑1 (MQ=255)
                                                                                 tgtCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTg        <  2:1996280/93‑1 (MQ=255)
                                                                                 tgtCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTg        <  1:3541818/93‑1 (MQ=255)
                                                                                  gtCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGc       <  1:4670853/93‑1 (MQ=255)
                                                                                  gtCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGc       <  2:2653408/93‑1 (MQ=255)
                                                                                  gtCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGc       <  2:2428082/93‑1 (MQ=255)
                                                                                   tCATGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGcc      >  1:219332/1‑93 (MQ=255)
                                                                                     aTGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCAc    >  2:4172095/1‑93 (MQ=255)
                                                                                     aTGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCAc    >  2:2549167/1‑93 (MQ=255)
                                                                                     aTGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCAc    >  2:1419431/1‑93 (MQ=255)
                                                                                      tGGTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACt   >  1:1147743/1‑93 (MQ=255)
                                                                                       ggTGCTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTg  <  1:1938499/93‑1 (MQ=255)
                                                                                          |                                                                                         
ACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTG  >  NC_000913/1019072‑1019251

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