Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 459,573 A→G I488V (ATT→GTT)  lon → ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913459,5730AG86.7% 714.8 / 94.4 293I488V (ATT→GTT) lonATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base A (22/17);  new base G (110/144);  total (132/161)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.70e-01

Only 100 of 298 total aligned reads displayed.

TGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGATTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGCACTGAAAAAA  >  NC_000913/459481‑459665
                                                                                            |                                                                                            
   aaGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTc                                                                                           <  1:7247944/93‑1 (MQ=255)
      tGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGcc                                                                                        >  1:4616242/1‑93 (MQ=255)
      tGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGcc                                                                                        >  2:1239281/1‑93 (MQ=255)
        gATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCctc                                                                                      <  2:2187684/93‑1 (MQ=255)
          ttACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCtctc                                                                                    <  1:6139571/93‑1 (MQ=255)
          ttACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCtctc                                                                                    >  2:6519792/1‑93 (MQ=255)
           tACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTcc                                                                                   >  2:4565104/1‑93 (MQ=255)
             nnaTCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCgg                                                                                 <  2:1474892/91‑1 (MQ=255)
             cGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCgg                                                                                 <  1:7132385/93‑1 (MQ=255)
              gATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGATTCGCCTCTCCGGt                                                                                >  2:6586374/1‑93 (MQ=255)
               aTCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGATTCGCCTCTCCGGtt                                                                               >  1:2887817/1‑93 (MQ=255)
               aTCTCAGCGACGCGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGtt                                                                               <  1:5579599/93‑1 (MQ=255)
                tctcAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTta                                                                              >  2:6015498/1‑93 (MQ=255)
                tctcAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTta                                                                              >  2:2570398/1‑93 (MQ=255)
                tctcAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTta                                                                              <  2:5051662/93‑1 (MQ=255)
                tctcAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGATTCGCCTCTCCGGTta                                                                              <  1:3121419/93‑1 (MQ=255)
                 ctcAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTtat                                                                             <  1:5010466/93‑1 (MQ=255)
                 ctcAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTtat                                                                             <  1:5490027/93‑1 (MQ=255)
                  tcAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTtata                                                                            <  2:5542605/93‑1 (MQ=255)
                    aGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATAcc                                                                          <  2:5656890/93‑1 (MQ=255)
                      cGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGa                                                                        <  2:7369079/93‑1 (MQ=255)
                      cGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGa                                                                        <  1:1053560/93‑1 (MQ=255)
                      cGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGa                                                                        <  2:3364230/93‑1 (MQ=255)
                       gACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGaa                                                                       >  2:2823641/1‑93 (MQ=255)
                       gACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGaa                                                                       >  2:4070208/1‑93 (MQ=255)
                       gACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGATTCGCCTCTCCGGTTATACCGaa                                                                       <  1:6956628/93‑1 (MQ=255)
                         cGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGATTCGCCTCTCCGGTTATACCGAAGa                                                                     <  2:3452720/93‑1 (MQ=255)
                          gTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGAt                                                                    >  1:5684314/1‑93 (MQ=255)
                          gTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGAt                                                                    >  1:7513975/1‑93 (MQ=255)
                           tGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATg                                                                   >  1:4673647/1‑93 (MQ=255)
                           tGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATg                                                                   >  2:4267650/1‑93 (MQ=255)
                           tGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATg                                                                   >  1:1021994/1‑93 (MQ=255)
                           tGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATg                                                                   >  1:2294042/1‑93 (MQ=255)
                           tGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATg                                                                   >  1:6405878/1‑93 (MQ=255)
                               gTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGaaaa                                                               >  2:7025587/1‑93 (MQ=255)
                               gTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGaaaa                                                               <  2:4503291/93‑1 (MQ=255)
                               gTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGATTCGCCTCTCCGGTTATACCGAAGATGaaaa                                                               >  1:2397057/1‑93 (MQ=255)
                                tttGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGaaaaa                                                              >  1:4441869/1‑93 (MQ=255)
                                     cgcgACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGATTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGaa                                                         <  2:1915997/93‑1 (MQ=255)
                                      gcgACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAAc                                                        <  2:544568/93‑1 (MQ=255)
                                      gcgACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAAc                                                        <  1:4069654/93‑1 (MQ=255)
                                      gcgACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAAc                                                        >  1:2064110/1‑93 (MQ=255)
                                        gACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACAt                                                      <  2:5988252/93‑1 (MQ=255)
                                        gACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACAt                                                      >  1:6771513/1‑93 (MQ=255)
                                        gACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACAt                                                      >  1:7350642/1‑93 (MQ=255)
                                         aCGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATc                                                     >  2:3099597/1‑93 (MQ=255)
                                         aCGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATc                                                     >  1:4409158/1‑93 (MQ=255)
                                         aCGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGATTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATc                                                     <  2:68169/93‑1 (MQ=255)
                                           gTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGc                                                   >  2:3962346/1‑93 (MQ=255)
                                           gTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGc                                                   <  1:974858/93‑1 (MQ=255)
                                               aaCTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGATTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCaaa                                               >  1:1515109/1‑93 (MQ=255)
                                                aCTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAAc                                              <  1:2493636/93‑1 (MQ=255)
                                                aCTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAAc                                              >  1:332403/1‑93 (MQ=255)
                                                  tCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGt                                            >  2:1208201/1‑93 (MQ=255)
                                                  tCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGt                                            >  1:2713005/1‑93 (MQ=255)
                                                  tCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGATTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGt                                            >  2:1615508/1‑93 (MQ=255)
                                                     aTGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCAc                                         <  1:487141/93‑1 (MQ=255)
                                                     aTGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCAc                                         <  1:6007183/93‑1 (MQ=255)
                                                      tGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTta                                                                              <  1:2955262/55‑1 (MQ=255)
                                                      tGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTta                                                                              >  2:2955262/1‑55 (MQ=255)
                                                       gAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACct                                       <  1:4190192/93‑1 (MQ=255)
                                                       gAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACct                                       <  2:2126509/93‑1 (MQ=255)
                                                         nnaTTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACctgc                                     <  1:2198256/91‑1 (MQ=255)
                                                         aCATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACctgc                                     <  1:4732883/93‑1 (MQ=255)
                                                          taTTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTGACctgct                                    <  2:179075/92‑1 (MQ=255)
                                                          cATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACctgct                                    >  1:6499027/1‑93 (MQ=255)
                                                              ccGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCg                                <  2:4435861/93‑1 (MQ=255)
                                                                ggCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGaa                              <  1:5181417/93‑1 (MQ=255)
                                                                ggCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGaa                              <  1:6821260/93‑1 (MQ=255)
                                                                ggCACCGCTGCTGGATCGTATGGAAGTGATTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGaa                              <  2:2017963/93‑1 (MQ=255)
                                                                 gCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAg                             <  1:1300482/93‑1 (MQ=255)
                                                                  cACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGc                            >  2:4781279/1‑93 (MQ=255)
                                                                    ccGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAg                          <  2:5285361/93‑1 (MQ=255)
                                                                     cGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTTAACATCGCCAAACGTCACCTGCTGCCGAAGCAGa                         <  2:6756689/93‑1 (MQ=255)
                                                                       ctgctgGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGAtt                       >  1:3912894/1‑93 (MQ=255)
                                                                        tgctgGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTg                      >  1:3794242/1‑93 (MQ=255)
                                                                        tgctgGATCGTATGGAAGTGGTTCGCCTCTCCCGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTg                      >  1:5973432/1‑93 (MQ=255)
                                                                         gctgGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGa                     <  1:3598041/93‑1 (MQ=255)
                                                                         gctgGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGa                     <  2:640707/93‑1 (MQ=255)
                                                                          ctgGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGaa                    >  2:6894169/1‑93 (MQ=255)
                                                                          ctgGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGaa                    <  2:3668669/93‑1 (MQ=255)
                                                                          ctgGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGaa                    <  1:6560234/93‑1 (MQ=255)
                                                                           tgGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAAc                   <  1:5924781/93‑1 (MQ=255)
                                                                           tgGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAAc                   >  1:6250691/1‑93 (MQ=255)
                                                                           tgGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAAc                   >  2:795075/1‑93 (MQ=255)
                                                                             gATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGt                 <  2:2446566/93‑1 (MQ=255)
                                                                               tCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTaa               <  2:7174758/93‑1 (MQ=255)
                                                                               tCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTaa               <  2:249240/93‑1 (MQ=255)
                                                                                cGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAAt              <  2:2382738/93‑1 (MQ=255)
                                                                                cGTATGGAAGTGATTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAAt              <  1:399063/93‑1 (MQ=255)
                                                                                 gTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATg             >  2:4145279/1‑93 (MQ=255)
                                                                                  tATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGc            >  2:6224242/1‑93 (MQ=255)
                                                                                   aTGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGCa           <  1:5743226/93‑1 (MQ=255)
                                                                                   aTGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGCa           <  2:854629/93‑1 (MQ=255)
                                                                                    tGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGCAc          >  2:2701126/1‑93 (MQ=255)
                                                                                      gAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGCACTg        >  1:1831787/1‑93 (MQ=255)
                                                                                      gAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGCACTg        >  2:7425965/1‑93 (MQ=255)
                                                                                         gTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGCACTGaaa     <  1:5428716/93‑1 (MQ=255)
                                                                                         gTGATTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGCACTGaaa     <  1:1892095/93‑1 (MQ=255)
                                                                                            aTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGCACTGaaaaaa  <  2:3793125/93‑1 (MQ=255)
                                                                                            |                                                                                            
TGGAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGATTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGCACTGAAAAAA  >  NC_000913/459481‑459665

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