Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 908,654 A→G L540P (CTC→CCC)  poxB ← pyruvate oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913908,6540AG87.7% 701.8 / 76.7 277L540P (CTC→CCC) poxBpyruvate oxidase
Reads supporting (aligned to +/- strand):  ref base A (19/15);  new base G (122/121);  total (141/136)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 280 total aligned reads displayed.

GCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCT  >  NC_000913/908562‑908746
                                                                                            |                                                                                            
   aGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGtt                                                                                           >  1:6061054/1‑93 (MQ=255)
   aGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGtt                                                                                           >  2:5276122/1‑93 (MQ=255)
     tttgtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTg                                                                                         <  1:2846047/93‑1 (MQ=255)
     tttgtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTg                                                                                         >  2:4745820/1‑93 (MQ=255)
     tttgtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTg                                                                                         >  1:6141582/1‑93 (MQ=255)
      ttgtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGa                                                                                        <  2:6732246/93‑1 (MQ=255)
       tgtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGAt                                                                                       <  1:5681568/93‑1 (MQ=255)
       tgtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGAt                                                                                       <  1:6531209/93‑1 (MQ=255)
       tgtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGAt                                                                                       <  1:7210546/93‑1 (MQ=255)
        gtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATc                                                                                      <  2:5525619/93‑1 (MQ=255)
        gtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATc                                                                                      <  1:2922734/93‑1 (MQ=255)
           ttCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGc                                                                                   <  1:3249874/93‑1 (MQ=255)
             cGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCgg                                                                                 <  1:4593769/93‑1 (MQ=255)
             cGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCgg                                                                                 <  2:1029688/93‑1 (MQ=255)
              gCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGt                                                                                <  1:3415958/93‑1 (MQ=255)
               ccAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTg                                                                               >  2:4455096/1‑93 (MQ=255)
                cAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTgg                                                                              >  1:3944378/1‑93 (MQ=255)
                   ttCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAAt                                                                           >  1:1890236/1‑93 (MQ=255)
                   ttCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAAt                                                                           <  2:2925499/93‑1 (MQ=255)
                     cGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATgg                                                                         <  1:6960595/93‑1 (MQ=255)
                      gATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGc                                                                        <  2:6655520/93‑1 (MQ=255)
                      gATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGc                                                                        >  1:7328312/1‑93 (MQ=255)
                       aTCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCt                                                                       >  1:5841484/1‑93 (MQ=255)
                         cACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTaa                                                                     <  1:124158/93‑1 (MQ=255)
                         cACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTaa                                                                     >  1:5997302/1‑93 (MQ=255)
                         cACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTaa                                                                     <  2:6493241/93‑1 (MQ=255)
                           cTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAAct                                                                   >  2:1592905/1‑93 (MQ=255)
                             tcatcaCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAActct                                                                 >  2:4198093/1‑93 (MQ=255)
                               atcaCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACtcttc                                                               <  2:5761246/93‑1 (MQ=255)
                                tcaCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACtcttct                                                              >  2:2278361/1‑93 (MQ=255)
                                tcaCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACtcttct                                                              >  1:3504039/1‑93 (MQ=255)
                                tcaCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACtcttct                                                              <  1:2535874/93‑1 (MQ=255)
                                tcaCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACtcttct                                                              <  1:5063358/93‑1 (MQ=255)
                                tcaCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACtcttct                                                              <  1:5398752/93‑1 (MQ=255)
                                 caCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTcttctt                                                             <  2:842403/93‑1 (MQ=255)
                                  aCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCttt                                                            <  2:5664397/93‑1 (MQ=255)
                                   ccGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTg                                                           <  1:3725427/93‑1 (MQ=255)
                                   ccGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTg                                                           >  1:6738992/1‑93 (MQ=255)
                                    cgcgTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTgg                                                          >  1:434898/1‑93 (MQ=255)
                                     gcgTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGc                                                         <  1:5482670/93‑1 (MQ=255)
                                       gTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGa                                                       >  2:3486220/1‑93 (MQ=255)
                                        tCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGac                                                      >  1:3781532/1‑93 (MQ=255)
                                          cGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGacca                                                    >  1:4083562/1‑93 (MQ=255)
                                          cGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGacca                                                    >  1:4493396/1‑93 (MQ=255)
                                           gCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGaccac                                                   >  1:5318558/1‑93 (MQ=255)
                                            cTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGaccacc                                                  >  1:201452/1‑93 (MQ=255)
                                             tgatgaTTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGAccacca                                                 >  2:1879213/1‑93 (MQ=255)
                                             tgatgaTTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGAccacca                                                 >  1:917613/1‑93 (MQ=255)
                                              gatgatTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACcaccac                                                >  2:6829842/1‑93 (MQ=255)
                                              gatgatTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACcaccac                                                >  2:6161244/1‑93 (MQ=255)
                                               atgatTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACcaca                                               <  2:578323/93‑1 (MQ=255)
                                               atgatTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACcaca                                               <  1:2775283/93‑1 (MQ=255)
                                                tgatTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACAt                                              >  1:7353466/1‑93 (MQ=255)
                                                  atTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATcc                                            <  2:2807731/93‑1 (MQ=255)
                                                  atTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATcc                                            >  1:4261416/1‑93 (MQ=255)
                                                  atTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATcc                                            <  1:2724145/93‑1 (MQ=255)
                                                  atTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATcc                                            >  1:4434195/1‑93 (MQ=255)
                                                    tGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATccac                                          <  2:5191215/93‑1 (MQ=255)
                                                      cgcgcAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATccacca                                        <  2:6316612/93‑1 (MQ=255)
                                                       gcgcAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCaa                                       >  2:3218666/1‑93 (MQ=255)
                                                       gcgcAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCaa                                       >  2:2850901/1‑93 (MQ=255)
                                                       gcgcAGCATATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCaa                                       <  2:1995098/93‑1 (MQ=255)
                                                         gcagcaTATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAaca                                     <  1:7467466/93‑1 (MQ=255)
                                                          cagcaTATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCGCCACATCCACCAacac                                    <  2:2045124/93‑1 (MQ=255)
                                                            gcaTATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCg                                  >  1:3639226/1‑93 (MQ=255)
                                                            gcaTATACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCg                                  >  1:7120885/1‑93 (MQ=255)
                                                              atatACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGa                                >  2:596785/1‑93 (MQ=255)
                                                              atatACAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGa                                <  2:7049843/93‑1 (MQ=255)
                                                               tataCAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGAc                               >  2:1329448/1‑93 (MQ=255)
                                                               tataCAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGAc                               <  2:3749719/93‑1 (MQ=255)
                                                                 taCAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCg                             <  1:2998085/93‑1 (MQ=255)
                                                                   cAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTc                           <  1:1093639/93‑1 (MQ=255)
                                                                   cAGGCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTc                           >  1:3327304/1‑93 (MQ=255)
                                                                     ggCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGa                         <  2:254646/93‑1 (MQ=255)
                                                                     ggCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGa                         <  2:5414992/93‑1 (MQ=255)
                                                                      gCTGAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGAt                        <  2:459107/93‑1 (MQ=255)
                                                                       cTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATg                       <  2:7449689/93‑1 (MQ=255)
                                                                         gAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGga                     >  2:3304957/1‑93 (MQ=255)
                                                                         gAAACCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGga                     <  1:6780372/93‑1 (MQ=255)
                                                                          aaaCCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGgag                    >  1:4662615/1‑93 (MQ=255)
                                                                           aaCCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGgaga                   <  1:2273991/93‑1 (MQ=255)
                                                                           aaCCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGgaga                   >  2:643642/1‑93 (MQ=255)
                                                                           aaCCTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGgaga                   >  1:5724144/1‑93 (MQ=255)
                                                                              cTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAgg                <  1:6444294/93‑1 (MQ=255)
                                                                              cTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAgg                <  2:5861869/93‑1 (MQ=255)
                                                                              cTTTGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAgg                <  2:6190167/93‑1 (MQ=255)
                                                                               tttGGCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGgc               >  1:6038549/1‑93 (MQ=255)
                                                                                  ggCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGcgcg            <  1:1353414/93‑1 (MQ=255)
                                                                                  ggCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGcgcg            <  1:1920466/93‑1 (MQ=255)
                                                                                   gCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGt           <  1:2429818/93‑1 (MQ=255)
                                                                                   gCCTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGt           <  1:4191486/93‑1 (MQ=255)
                                                                                    ccTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGtt          >  2:236687/1‑93 (MQ=255)
                                                                                    ccTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGtt          >  2:4917616/1‑93 (MQ=255)
                                                                                     cTGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTg         >  2:3962921/1‑93 (MQ=255)
                                                                                      tGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGc        <  1:6334597/93‑1 (MQ=255)
                                                                                      tGTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGc        >  1:4913778/1‑93 (MQ=255)
                                                                                       gTTCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCa       <  1:4745820/93‑1 (MQ=255)
                                                                                        ttCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAg      >  2:4281606/1‑93 (MQ=255)
                                                                                         tCGGGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAgg     <  2:3127845/93‑1 (MQ=255)
                                                                                            aGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCt  <  1:1483067/93‑1 (MQ=255)
                                                                                            |                                                                                            
GCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCT  >  NC_000913/908562‑908746

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