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,5730AG99.7% 1030.3 / ‑6.2 300I488V (ATT→GTT) lonATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (162/137);  total (162/138)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 306 total aligned reads displayed.

GAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGATTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGCACTGAAA  >  NC_000913/459483‑459662
                                                                                          |                                                                                         
gAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGtt                                                                                         >  2:1391045/1‑93 (MQ=255)
gAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGtt                                                                                         <  2:5925339/93‑1 (MQ=255)
  aGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCg                                                                                       >  2:5252479/1‑93 (MQ=255)
   gTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGc                                                                                      >  2:3246350/1‑93 (MQ=255)
   gTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGc                                                                                      <  2:6289294/93‑1 (MQ=255)
       aTTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCctct                                                                                  >  2:8382901/1‑93 (MQ=255)
       aTTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCctct                                                                                  >  2:4589236/1‑93 (MQ=255)
        ttACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCtctc                                                                                 <  1:1158948/93‑1 (MQ=255)
        ttACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCtctc                                                                                 <  1:4813599/93‑1 (MQ=255)
        ttACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCtctc                                                                                 >  1:379773/1‑93 (MQ=255)
           cGATCTTAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCgg                                                                              <  1:873983/93‑1 (MQ=255)
            gATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGt                                                                             <  2:6470303/93‑1 (MQ=255)
             aTCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGtt                                                                            <  1:7172090/93‑1 (MQ=255)
              tctcAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTta                                                                           <  2:2698744/93‑1 (MQ=255)
              tctcAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTta                                                                           >  2:5981379/1‑93 (MQ=255)
               ctcAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTtat                                                                          >  1:3093869/1‑93 (MQ=255)
                tcAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTtata                                                                         >  1:625011/1‑93 (MQ=255)
                tcAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTtata                                                                         <  1:188592/93‑1 (MQ=255)
                 cAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATAc                                                                        >  1:5374675/1‑93 (MQ=255)
                   gCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCg                                                                      >  1:4718851/1‑93 (MQ=255)
                      aCGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAg                                                                   >  2:7487068/1‑93 (MQ=255)
                      aCGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAg                                                                   >  1:6024118/1‑93 (MQ=255)
                      aCGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGGTATACCGAAg                                                                   >  2:4253/1‑93 (MQ=255)
                       cGTTATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGa                                                                  <  2:2552141/93‑1 (MQ=255)
                       cGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGa                                                                  <  2:2283128/93‑1 (MQ=255)
                       cGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGa                                                                  <  2:2868838/93‑1 (MQ=255)
                         tGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATg                                                                <  1:4541643/93‑1 (MQ=255)
                         tGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATg                                                                >  1:987555/1‑93 (MQ=255)
                         tGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATg                                                                >  1:3610096/1‑93 (MQ=255)
                         tGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATg                                                                <  2:3654626/93‑1 (MQ=255)
                           aTGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGaa                                                              >  2:1837732/1‑93 (MQ=255)
                             gTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGaaaa                                                            >  2:3474608/1‑93 (MQ=255)
                               ttGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAAc                                                          <  1:5571272/93‑1 (MQ=255)
                                   cgcgACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGaa                                                      <  1:5019998/93‑1 (MQ=255)
                                      gACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACAt                                                   >  1:1020309/1‑93 (MQ=255)
                                       aCGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTTTACCGAAGATGAAAAACTGAACATc                                                  >  1:8229749/1‑93 (MQ=255)
                                       aCGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATc                                                  >  1:6486686/1‑93 (MQ=255)
                                       aCGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATc                                                  <  2:2491160/93‑1 (MQ=255)
                                       aCGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATc                                                  <  1:6096067/93‑1 (MQ=255)
                                       aCGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATc                                                  >  1:2878018/1‑93 (MQ=255)
                                       aCGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATc                                                  >  2:5584249/1‑93 (MQ=255)
                                       aCGTCGAACTCCATGAACATTCAGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGACAAACTGAACATc                                                  >  2:5178853/1‑93 (MQ=255)
                                        cGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATcc                                                 >  1:6078800/1‑92 (MQ=255)
                                        cGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCg                                                 >  1:6925499/1‑93 (MQ=255)
                                          tCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGcc                                               >  1:7429905/1‑93 (MQ=255)
                                           cGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTTGACGGTTATACCGAAGATGCAAACATGAACATCGCCa                                              >  1:4230425/1‑93 (MQ=255)
                                           cGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCa                                              >  1:8659735/1‑93 (MQ=255)
                                            gAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCaa                                             >  2:1452974/1‑93 (MQ=255)
                                             aaCTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCaaa                                            >  2:8403481/1‑93 (MQ=255)
                                              aCTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAAc                                           <  1:5144409/93‑1 (MQ=255)
                                              aCTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAAc                                           >  1:2580084/1‑93 (MQ=255)
                                                tCCATGAACATTCCTGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATAACGAAGATGAACAACTGAACATCGCCAAAAGt                                         >  2:6905615/1‑93 (MQ=255)
                                                tCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACgg                                         >  2:3113627/1‑92 (MQ=255)
                                                 ccATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGta                                        >  2:8660599/1‑92 (MQ=255)
                                                 ccATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTc                                        >  1:8809766/1‑93 (MQ=255)
                                                 ccATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTc                                        <  2:7647076/93‑1 (MQ=255)
                                                  cATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCa                                       <  2:4942279/93‑1 (MQ=255)
                                                     gAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACct                                    >  1:2034891/1‑93 (MQ=255)
                                                     gAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACct                                    <  1:4412364/93‑1 (MQ=255)
                                                     gAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACct                                    <  2:205080/93‑1 (MQ=255)
                                                      aaCATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACctg                                   >  1:1329469/1‑93 (MQ=255)
                                                       aCATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACctgc                                  >  2:7801812/1‑93 (MQ=255)
                                                        cATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACctgct                                 <  2:2978341/93‑1 (MQ=255)
                                                        cATTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACctgct                                 <  1:7988725/93‑1 (MQ=255)
                                                         aTTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACctgctg                                >  2:8217363/1‑93 (MQ=255)
                                                         aTTCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACctgctg                                >  2:7057257/1‑93 (MQ=255)
                                                          ttCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCtgctgc                               <  1:270056/93‑1 (MQ=255)
                                                           tCCGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGcc                              <  2:2898985/93‑1 (MQ=255)
                                                             cGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGa                            >  1:1745380/1‑93 (MQ=255)
                                                             cGGCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGa                            >  2:1182777/1‑93 (MQ=255)
                                                              ggCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGaa                           >  2:8034336/1‑93 (MQ=255)
                                                              ggCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGaa                           <  2:2795264/93‑1 (MQ=255)
                                                              ggCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGaa                           <  1:2376156/93‑1 (MQ=255)
                                                              ggCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGaa                           >  2:6168862/1‑93 (MQ=255)
                                                               gCACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAg                          <  2:721746/93‑1 (MQ=255)
                                                                cACCGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGc                         <  1:4184137/93‑1 (MQ=255)
                                                                   cGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGa                      <  2:8013050/93‑1 (MQ=255)
                                                                   cGCTGCTGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGa                      >  2:6592476/1‑93 (MQ=255)
                                                                    gctgctGGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGAt                     >  2:4984415/1‑93 (MQ=255)
                                                                     ctgctgGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGAtt                    <  2:4469332/93‑1 (MQ=255)
                                                                        gtgGATCGTATGGAAGTGATTCGCCTCGCCGGTTATACCGAAGATGAAAAACTGAACTTCGCCAAACGTCACCTGCTGCCGAAGCAGATTGaa                 <  2:3569483/92‑1 (MQ=255)
                                                                        ctgGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGaa                 <  2:4410254/93‑1 (MQ=255)
                                                                        ctgGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGaa                 <  2:3805934/93‑1 (MQ=255)
                                                                         tgGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAAc                >  1:6894941/1‑93 (MQ=255)
                                                                         tgGATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAAc                <  2:1572272/93‑1 (MQ=255)
                                                                          ggATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGACCATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACg               >  1:7662111/1‑93 (MQ=255)
                                                                          ggAGCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACg               <  2:3432801/93‑1 (MQ=255)
                                                                           gATCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGt              >  2:1981300/1‑93 (MQ=255)
                                                                             tCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTaa            <  1:6657117/93‑1 (MQ=255)
                                                                             tCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTaa            <  2:6972865/93‑1 (MQ=255)
                                                                             tCGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTaa            <  1:5724158/93‑1 (MQ=255)
                                                                              cGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAAt           <  1:3307116/93‑1 (MQ=255)
                                                                              cGTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAAt           <  1:8375574/93‑1 (MQ=255)
                                                                               gTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATg          >  1:1626351/1‑93 (MQ=255)
                                                                               gTATGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATg          <  2:4098995/93‑1 (MQ=255)
                                                                                  tGGAAGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGCAc       >  1:5213456/1‑93 (MQ=255)
                                                                                    gAAGGGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGCACTg     >  1:2272613/1‑93 (MQ=255)
                                                                                     aaGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGCACTGa    >  2:4817367/1‑93 (MQ=255)
                                                                                     aaGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGCACTGa    >  2:8577008/1‑93 (MQ=255)
                                                                                      aGTGGTTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGCACTGaa   >  2:97970/1‑93 (MQ=255)
                                                                                          |                                                                                         
GAAGTGGATTACGATCTCAGCGACGTGATGTTTGTCGCGACGTCGAACTCCATGAACATTCCGGCACCGCTGCTGGATCGTATGGAAGTGATTCGCCTCTCCGGTTATACCGAAGATGAAAAACTGAACATCGCCAAACGTCACCTGCTGCCGAAGCAGATTGAACGTAATGCACTGAAA  >  NC_000913/459483‑459662

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: