Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 254,298 T→C T473T (ACA→ACG pepD ← aminoacyl‑histidine dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913254,2980TC83.4% 584.7 / 113.5 265T473T (ACA→ACGpepDaminoacyl‑histidine dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (19/25);  new base C (123/98);  total (142/123)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

Only 100 of 270 total aligned reads displayed.

GGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCAT  >  NC_000913/254207‑254387
                                                                                           |                                                                                         
ggATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAgtg                                                                                          <  2:3855336/93‑1 (MQ=255)
 gATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAgcgt                                                                                         <  1:1847917/93‑1 (MQ=255)
  aTGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTAAGTCAGCAGCGTc                                                                                        <  1:1719057/93‑1 (MQ=255)
     aaaCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAg                                                                                     >  1:6982291/1‑93 (MQ=255)
      aaCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGt                                                                                    <  1:5879504/93‑1 (MQ=255)
      aaCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGt                                                                                    <  1:2685555/93‑1 (MQ=255)
       aCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTa                                                                                   >  2:4718991/1‑93 (MQ=255)
        cTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTaa                                                                                  <  1:4067320/93‑1 (MQ=255)
         tGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAAt                                                                                 >  2:6355155/1‑93 (MQ=255)
         tGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAAt                                                                                 >  2:450340/1‑93 (MQ=255)
         tGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAAt                                                                                 <  1:1428842/93‑1 (MQ=255)
          gCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATg                                                                                >  2:2894024/1‑93 (MQ=255)
            cGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGAc                                                                              <  2:3915850/93‑1 (MQ=255)
            cGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGGAATGAc                                                                              >  1:3671236/1‑93 (MQ=255)
             ggAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGAcc                                                                             <  1:5183489/93‑1 (MQ=255)
              gAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCt                                                                            <  1:7422431/93‑1 (MQ=255)
               aaGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTa                                                                           <  2:2738468/93‑1 (MQ=255)
               aaGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTa                                                                           <  2:3496870/93‑1 (MQ=255)
                aGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTAc                                                                          >  2:6911317/1‑93 (MQ=255)
                aGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTAc                                                                          >  2:4951360/1‑93 (MQ=255)
                aGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTAc                                                                          >  2:3437840/1‑93 (MQ=255)
                  gCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGc                                                                        <  2:7486838/93‑1 (MQ=255)
                  gCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGc                                                                        >  1:4107827/1‑93 (MQ=255)
                     aTTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCttt                                                                     <  1:6749813/93‑1 (MQ=255)
                       tAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCg                                                                   <  1:4176565/93‑1 (MQ=255)
                        aaaCGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGa                                                                  <  1:4990428/93‑1 (MQ=255)
                         aaCGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCGGTAATGACCTACGCTTTCGAt                                                                 >  2:5869467/1‑93 (MQ=255)
                          aCGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGAtg                                                                <  2:1088889/93‑1 (MQ=255)
                          aCGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGAtg                                                                >  1:2994018/1‑93 (MQ=255)
                           cGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtgt                                                               <  2:2282541/93‑1 (MQ=255)
                           cGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGAtgt                                                               <  1:4285158/93‑1 (MQ=255)
                            gCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGAtgtg                                                              >  2:1813525/1‑93 (MQ=255)
                             ccATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGa                                                             >  1:6584213/1‑93 (MQ=255)
                               aTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAAc                                                           <  2:2064952/93‑1 (MQ=255)
                               aTCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAAc                                                           >  1:620551/1‑93 (MQ=255)
                                tCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACt                                                          >  2:5970623/1‑93 (MQ=255)
                                 ccGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACtt                                                         <  1:1552135/93‑1 (MQ=255)
                                  cGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTg                                                        >  1:935093/1‑93 (MQ=255)
                                  cGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTg                                                        >  1:5389912/1‑93 (MQ=255)
                                  cGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTg                                                        >  1:6040555/1‑93 (MQ=255)
                                   ggCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACATACGCTTTCGATGTGAACTTGc                                                       <  2:3318969/93‑1 (MQ=255)
                                    gcagcaAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCt                                                      >  1:87168/1‑93 (MQ=255)
                                     cagcaAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTc                                                     <  2:253487/93‑1 (MQ=255)
                                      agcaAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCa                                                    >  2:1058171/1‑93 (MQ=255)
                                      agcaAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCa                                                    >  2:7298455/1‑93 (MQ=255)
                                      agcaAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCa                                                    >  1:206583/1‑93 (MQ=255)
                                      agcaAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCa                                                    <  2:31405/93‑1 (MQ=255)
                                       gcaAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCAt                                                   >  1:1106463/1‑93 (MQ=255)
                                       gcaAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCAt                                                   >  2:1920149/1‑93 (MQ=255)
                                         aaatcaaatAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATcc                                                 <  1:4920819/93‑1 (MQ=255)
                                         aaatcaaatAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATcc                                                 >  1:5440934/1‑93 (MQ=255)
                                          aatcaaatAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCg                                                <  2:6603928/93‑1 (MQ=255)
                                            tcaaatAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGga                                              <  1:3951262/93‑1 (MQ=255)
                                             caaatAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGgag                                             <  1:4415074/93‑1 (MQ=255)
                                              aaatAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGgaga                                            <  1:7307149/93‑1 (MQ=255)
                                              aaatAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGgaga                                            >  1:2570205/1‑93 (MQ=255)
                                              aaatAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGgaga                                            <  2:1493481/93‑1 (MQ=255)
                                                 taatTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAgtg                                         >  1:1788178/1‑93 (MQ=255)
                                                  aatTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAgtgt                                        >  2:216028/1‑93 (MQ=255)
                                                  aatTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAgtgt                                        >  2:5788730/1‑93 (MQ=255)
                                                  aatTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAgtgt                                        <  1:4605146/93‑1 (MQ=255)
                                                    ttACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTgg                                      >  2:3003262/1‑93 (MQ=255)
                                                    ttACTTCGCCGGAATTTCTCTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTgg                                      <  1:2740161/93‑1 (MQ=255)
                                                     tACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGa                                     >  1:2456781/1‑93 (MQ=255)
                                                      aCTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGAc                                    >  2:1190183/1‑93 (MQ=255)
                                                      aCTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGAc                                    <  2:5322297/93‑1 (MQ=255)
                                                            ccGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGTCCGGTGa                              >  2:1662609/1‑93 (MQ=255)
                                                            ccGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGa                              >  1:3460514/1‑93 (MQ=255)
                                                            ccGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGa                              <  1:4450174/93‑1 (MQ=255)
                                                             cGGAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGAt                             >  1:3294856/1‑93 (MQ=255)
                                                              ggAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATa                            >  1:5948950/1‑93 (MQ=255)
                                                              ggAATTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATa                            <  1:1260846/93‑1 (MQ=255)
                                                                 aTTTCTTTCAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGtt                         <  2:4071422/93‑1 (MQ=255)
                                                                   ttctttcAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTgg                       <  2:3241804/93‑1 (MQ=255)
                                                                     ctttcAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGcc                     <  1:316474/93‑1 (MQ=255)
                                                                     ctttcAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGcc                     >  1:6848219/1‑93 (MQ=255)
                                                                        tcAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGa                  >  1:4808965/1‑93 (MQ=255)
                                                                        tcAGCAGTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGa                  >  2:6660190/1‑93 (MQ=255)
                                                                          agcagTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATa                <  1:2352659/93‑1 (MQ=255)
                                                                           gcagTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAg               >  2:7187798/1‑93 (MQ=255)
                                                                            cagTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa              <  2:2803948/93‑1 (MQ=255)
                                                                            cagTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa              <  1:5707895/93‑1 (MQ=255)
                                                                             agTTCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGCGGACCGGTGATAGTTGGCCCGATAGaa             <  2:5164395/93‑1 (MQ=255)
                                                                               ttCAGTCAGCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc           >  2:6460715/1‑93 (MQ=255)
                                                                                tcagtcagCAGTGTCCAGTAATGACCTACACTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAcc          <  2:3661839/93‑1 (MQ=255)
                                                                                tcagtcagCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAcc          >  1:1382399/1‑93 (MQ=255)
                                                                                tcagtcagCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAcc          >  1:7212086/1‑93 (MQ=255)
                                                                                 cagtcagCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCa         <  2:3285135/93‑1 (MQ=255)
                                                                                  agtcagCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt        >  1:6370549/1‑93 (MQ=255)
                                                                                   gtcagCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATg       >  2:6944208/1‑93 (MQ=255)
                                                                                   gtcagCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATg       <  1:377968/93‑1 (MQ=255)
                                                                                   gtcagCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATg       <  1:6547916/93‑1 (MQ=255)
                                                                                   gtcagCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATg       >  1:7242632/1‑93 (MQ=255)
                                                                                    tcagCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGt      >  1:1602387/1‑93 (MQ=255)
                                                                                    tcagCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGt      >  2:2513335/1‑93 (MQ=255)
                                                                                    tcagCAGCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGt      >  2:861613/1‑93 (MQ=255)
                                                                                       gcagTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCa   >  2:4550073/1‑93 (MQ=255)
                                                                                       gcagCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCa   >  2:6141544/1‑93 (MQ=255)
                                                                                       gcagCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCa   >  2:5671351/1‑93 (MQ=255)
                                                                                       gcagCGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCa   >  2:1711813/1‑93 (MQ=255)
                                                                                           |                                                                                         
GGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCAT  >  NC_000913/254207‑254387

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