Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 908,737 T→C Q512Q (CAA→CAG poxB ← pyruvate oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913908,7370TC99.5% 698.8 / NA 214Q512Q (CAA→CAGpoxBpyruvate oxidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (100/113);  minor base A (0/1);  total (100/114)
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 220 total aligned reads displayed.

CTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTG  >  NC_000913/908647‑908826
                                                                                          |                                                                                         
cTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTg                                                                                         >  2:2511918/1‑93 (MQ=255)
cTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTg                                                                                         >  2:1686829/1‑93 (MQ=255)
 tGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGc                                                                                        >  2:5994010/1‑93 (MQ=255)
  gTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCa                                                                                       >  1:6108229/1‑93 (MQ=255)
   ttCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAg                                                                                      >  2:63671/1‑93 (MQ=255)
   ttCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAg                                                                                      >  2:5563373/1‑93 (MQ=255)
     cGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCCCCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAggg                                                                                    <  2:1784200/93‑1 (MQ=255)
       aGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCt                                                                                  <  1:1109740/93‑1 (MQ=255)
       aGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCt                                                                                  <  2:3606374/93‑1 (MQ=255)
       aGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCt                                                                                  <  1:2478458/93‑1 (MQ=255)
        gTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCtt                                                                                 >  1:4916056/1‑93 (MQ=255)
        gTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCtt                                                                                 >  1:6635232/1‑93 (MQ=255)
         tttGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtc                                                                                >  2:910560/1‑93 (MQ=255)
         tttGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtc                                                                                <  2:2705951/93‑1 (MQ=255)
         tttGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtc                                                                                >  1:4050564/1‑93 (MQ=255)
           tGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtcat                                                                              <  1:1440428/93‑1 (MQ=255)
           tGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtcat                                                                              >  2:6111703/1‑93 (MQ=255)
            gATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtcatc                                                                             <  2:1085806/93‑1 (MQ=255)
             aTCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtcatca                                                                            >  2:745045/1‑93 (MQ=255)
             aTCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTtcatca                                                                            >  2:5463958/1‑93 (MQ=255)
              tCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCaa                                                                           <  2:4745758/93‑1 (MQ=255)
               cTGCGTTGGAATGGCTAACTCTTCTTTGGCGCCCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAAc                                                                          <  2:6182802/93‑1 (MQ=255)
               cTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAAc                                                                          <  1:5287194/93‑1 (MQ=255)
               cTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAAc                                                                          >  2:4321481/1‑93 (MQ=255)
                tGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACt                                                                         <  2:3003102/93‑1 (MQ=255)
                  cGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTc                                                                       <  1:4632795/93‑1 (MQ=255)
                  cGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTc                                                                       >  1:6170048/1‑93 (MQ=255)
                   ggTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCa                                                                      <  1:4888415/93‑1 (MQ=255)
                   ggTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCa                                                                      >  1:2562487/1‑93 (MQ=255)
                    gTGGAATGGCTAACTCTCCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAg                                                                     >  2:4034701/1‑93 (MQ=255)
                     tGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGa                                                                    <  2:5064378/93‑1 (MQ=255)
                        aaTGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGc                                                                 <  1:5264098/93‑1 (MQ=255)
                        aaTGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGc                                                                 <  2:3069228/93‑1 (MQ=255)
                        aaTGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGc                                                                 <  2:4970128/93‑1 (MQ=255)
                        aaTGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGc                                                                 <  2:2360653/93‑1 (MQ=255)
                         aTGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCt                                                                <  1:3290656/93‑1 (MQ=255)
                         aTGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCt                                                                >  1:3025234/1‑93 (MQ=255)
                         aTGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCt                                                                <  1:6440854/93‑1 (MQ=255)
                         aTGGCTAACCCTTCTTTGGCGACCACCACATCCACCAACGCCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCt                                                                <  1:1213054/93‑1 (MQ=255)
                          tGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCtt                                                               >  2:2181285/1‑93 (MQ=255)
                           gacTATCTCTTCTTTGGCGACCACCACATCCACCACCACCGTACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCttt                                                              <  1:79372/91‑1 (MQ=255)
                             cTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCttttt                                                            <  2:2630115/93‑1 (MQ=255)
                              tAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTc                                                           <  1:2427410/93‑1 (MQ=255)
                              tAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTc                                                           <  1:5175261/93‑1 (MQ=255)
                               aaCTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCt                                                          <  2:2015548/93‑1 (MQ=255)
                                aCTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTa                                                         <  2:2886542/93‑1 (MQ=255)
                                 ctctTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTac                                                        >  1:700113/1‑93 (MQ=255)
                                 ctctTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTac                                                        >  2:3664457/1‑93 (MQ=255)
                                 ctctTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTac                                                        <  2:1425586/93‑1 (MQ=255)
                                 ctctTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTac                                                        <  1:3592417/93‑1 (MQ=255)
                                  tcttctTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTaca                                                       >  1:298672/1‑93 (MQ=255)
                                   cttcttTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTacac                                                      <  2:5709394/93‑1 (MQ=255)
                                   cttcttTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTacac                                                      >  2:5892892/1‑93 (MQ=255)
                                     tcttTGGCGACCGCCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACgg                                                    <  2:1570196/93‑1 (MQ=255)
                                      cttTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGa                                                   >  1:4142714/1‑93 (MQ=255)
                                          ggCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATAcc                                               >  1:5536333/1‑93 (MQ=255)
                                           gCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATAccc                                              <  2:928484/93‑1 (MQ=255)
                                           gCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATAccc                                              >  1:5805100/1‑93 (MQ=255)
                                           gCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATAccc                                              >  2:190178/1‑93 (MQ=255)
                                             gACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGt                                            <  2:6161206/93‑1 (MQ=255)
                                               ccaccaCATCCACCAACGCCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATCCCCGTaa                                          <  1:3661951/93‑1 (MQ=255)
                                                caccacATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAAt                                         <  1:6065626/93‑1 (MQ=255)
                                                caccacATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAAt                                         <  1:3724961/93‑1 (MQ=255)
                                                caccacATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAAt                                         <  2:3630321/93‑1 (MQ=255)
                                                  ccacATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGc                                       >  1:2335222/1‑93 (MQ=255)
                                                     caTCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGc                                    >  1:190876/1‑93 (MQ=255)
                                                      aTCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCa                                   <  2:5185528/93‑1 (MQ=255)
                                                      aTCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGCCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCa                                   <  1:3879353/93‑1 (MQ=255)
                                                       tCCACCAACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCAc                                  >  1:1007250/1‑93 (MQ=255)
                                                          accaACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCt                               <  2:4276269/93‑1 (MQ=255)
                                                           ccaACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCtt                              >  1:4434432/1‑93 (MQ=255)
                                                            caACACCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTc                             <  2:6243325/93‑1 (MQ=255)
                                                              acacCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCgg                           >  1:4838152/1‑93 (MQ=255)
                                                              acacCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCgg                           >  2:5007154/1‑93 (MQ=255)
                                                              acacCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCgg                           <  1:5410029/93‑1 (MQ=255)
                                                               cacCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGc                          >  1:6379922/1‑93 (MQ=255)
                                                               cacCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGc                          >  2:814849/1‑93 (MQ=255)
                                                               cacCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGc                          >  2:2534019/1‑93 (MQ=255)
                                                               cacCGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGc                          <  2:6087802/93‑1 (MQ=255)
                                                                 ccGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCaa                        <  2:4811549/93‑1 (MQ=255)
                                                                 ccGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCaa                        >  2:1413647/1‑93 (MQ=255)
                                                                  cGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAAt                       <  2:2802629/93‑1 (MQ=255)
                                                                  cGGACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAAt                       <  2:1218783/93‑1 (MQ=255)
                                                                    gACCGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGc                     >  1:1304043/1‑93 (MQ=255)
                                                                      ccGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCgg                   <  2:515207/93‑1 (MQ=255)
                                                                      ccGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCgg                   >  1:2906858/1‑93 (MQ=255)
                                                                      ccGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCgg                   >  2:391190/1‑93 (MQ=255)
                                                                      ccGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCgg                   >  1:2776721/1‑93 (MQ=255)
                                                                       cGTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCggg                  >  2:4552572/1‑93 (MQ=255)
                                                                        gTCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGc                 <  2:374818/93‑1 (MQ=255)
                                                                         tCGATGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCa                >  2:659807/1‑93 (MQ=255)
                                                                            aTGGAGAAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAg             >  1:2116683/1‑93 (MQ=255)
                                                                              ggAGAAGGCGCGCTGCAGGGCTTCATCAACTGCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGtt           <  2:4111159/93‑1 (MQ=255)
                                                                                agaAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTtg         <  2:1491116/93‑1 (MQ=255)
                                                                                agaAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTtg         >  2:3229297/1‑93 (MQ=255)
                                                                                agaAGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTtg         <  1:956228/93‑1 (MQ=255)
                                                                                  aaGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTCtgtg       <  2:6277044/93‑1 (MQ=255)
                                                                                   aGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTgtgt      >  2:3360907/1‑93 (MQ=255)
                                                                                   aGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTgtgt      <  1:3482075/93‑1 (MQ=255)
                                                                                   aGGCGCGCTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTgtgt      >  1:4280562/1‑93 (MQ=255)
                                                                                          |                                                                                         
CTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGCACGCTTCGGCAATGCGGGCAAAGTTTGTGTCGTG  >  NC_000913/908647‑908826

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