Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 460,114 T→C V668A (GTA→GCA)  lon → ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913460,1140TC86.2% 522.6 / 71.7 203V668A (GTA→GCA) lonATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base T (15/13);  new base C (97/78);  total (112/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

Only 100 of 257 total aligned reads displayed.

CATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTG  >  NC_000913/460022‑460206
                                                                                            |                                                                                            
cATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATccacgtccacgt                                                                                              <  2:1758783/93‑1 (MQ=255)
 aTTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTa                                                                                             <  2:425488/93‑1 (MQ=255)
  ttCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCAc                                                                                            <  2:3350155/93‑1 (MQ=255)
  ttCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCAc                                                                                            <  2:511999/93‑1 (MQ=255)
    cAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCg                                                                                          >  2:302925/1‑93 (MQ=255)
    cAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCg                                                                                          >  1:5669788/1‑93 (MQ=255)
     aGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCgg                                                                                         >  1:2237231/1‑93 (MQ=255)
       gcggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGaa                                                                                       <  2:2406374/93‑1 (MQ=255)
       gcggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGaa                                                                                       <  2:1303225/93‑1 (MQ=255)
       gcggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGaa                                                                                       <  1:5921392/93‑1 (MQ=255)
         ggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAgg                                                                                     >  2:1729151/1‑93 (MQ=255)
         ggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAgg                                                                                     <  1:5167515/93‑1 (MQ=255)
         ggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAgg                                                                                     >  1:3636679/1‑93 (MQ=255)
           ggTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTg                                                                                   <  1:4362188/92‑1 (MQ=255)
            gTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGc                                                                                  >  2:2172728/1‑93 (MQ=255)
             ttAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCg                                                                                 >  1:5356243/1‑93 (MQ=255)
             ttAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCg                                                                                 >  1:2129394/1‑93 (MQ=255)
              tAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGa                                                                                >  2:6623183/1‑93 (MQ=255)
              tAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGa                                                                                <  1:671329/93‑1 (MQ=255)
               aaCGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGaa                                                                               >  2:657691/1‑92 (MQ=255)
               aaCGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGAc                                                                               >  2:4047982/1‑93 (MQ=255)
               aaCGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGAc                                                                               <  1:5856139/93‑1 (MQ=255)
                 cGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGc                                                                             >  2:4897244/1‑93 (MQ=255)
                  ggtggtTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGcc                                                                            >  1:852522/1‑93 (MQ=255)
                   gtggtTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCg                                                                           <  2:2242382/93‑1 (MQ=255)
                   gtggtTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCg                                                                           >  1:6632449/1‑93 (MQ=255)
                      gtTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGaaa                                                                        >  1:296665/1‑93 (MQ=255)
                      gtTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGaaa                                                                        <  2:5083941/93‑1 (MQ=255)
                      gtTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGaaa                                                                        <  2:4031584/93‑1 (MQ=255)
                        tCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGa                                                                      <  2:1234013/93‑1 (MQ=255)
                          gtgcgcgtgcgGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATg                                                                    >  1:4312078/1‑93 (MQ=255)
                          gtgcgcgtgcgGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATg                                                                    <  2:2635701/93‑1 (MQ=255)
                           tgcgcgtgcgGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATgg                                                                   >  2:5272608/1‑93 (MQ=255)
                           tgcgcgtgcgGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATgg                                                                   >  2:1706657/1‑93 (MQ=255)
                           tgcgcgtgcgGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATgg                                                                   >  1:1954061/1‑93 (MQ=255)
                           tgcgcgtgcgGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATgg                                                                   >  2:3408054/1‑93 (MQ=255)
                            gcgcgtgcgGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGt                                                                  >  1:4096544/1‑93 (MQ=255)
                             cgcgtgcgGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTGCACGCACCGGAAGGTGCGACGCCGAAAGATGGTc                                                                 <  2:451671/93‑1 (MQ=255)
                              gcgtgcgGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTcc                                                                >  1:2337595/1‑93 (MQ=255)
                              gcgtgcgGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTcc                                                                >  1:3982394/1‑93 (MQ=255)
                                gtgcgGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGCAAGATGGTCCGa                                                              >  1:1157474/1‑93 (MQ=255)
                                 tgcgGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAg                                                             >  2:4440260/1‑93 (MQ=255)
                                  gcgGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGt                                                            <  2:4190505/93‑1 (MQ=255)
                                   cgGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTg                                                           >  2:1952059/1‑93 (MQ=255)
                                    ggAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGc                                                          >  2:1076805/1‑93 (MQ=255)
                                     gAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGcc                                                         >  2:2923519/1‑93 (MQ=255)
                                     gAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGcc                                                         >  1:931844/1‑93 (MQ=255)
                                     gAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGcc                                                         >  2:2731092/1‑93 (MQ=255)
                                      aaaaaCTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCg                                                        >  2:3201956/1‑93 (MQ=255)
                                         aaCTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTa                                                     <  2:388939/93‑1 (MQ=255)
                                         aaCTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTa                                                     >  2:2217171/1‑93 (MQ=255)
                                           cTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTAtt                                                   <  2:353098/93‑1 (MQ=255)
                                             ggggATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGc                                                 >  1:1783475/1‑93 (MQ=255)
                                             ggggATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGc                                                 <  1:323882/93‑1 (MQ=255)
                                             ggggATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGc                                                 <  1:3692448/93‑1 (MQ=255)
                                              gggATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCt                                                >  2:6473316/1‑93 (MQ=255)
                                                gATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTAt                                              <  1:1234719/93‑1 (MQ=255)
                                                gATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTAt                                              >  2:3690334/1‑93 (MQ=255)
                                                 aTCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTAtg                                             >  2:3629258/1‑93 (MQ=255)
                                                    aaCCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGc                                          >  2:6074337/1‑93 (MQ=255)
                                                      cccTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCAc                                        >  1:1029768/1‑93 (MQ=255)
                                                      cccTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCAc                                        >  2:1373270/1‑93 (MQ=255)
                                                      cccTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCAc                                        >  1:525503/1‑93 (MQ=255)
                                                       ccTGATTTTTACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCAcc                                       <  2:5302561/93‑1 (MQ=255)
                                                        cTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACcg                                      >  1:4737435/1‑93 (MQ=255)
                                                           aTTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCgcgc                                   <  1:1440159/93‑1 (MQ=255)
                                                            tttttACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCt                                  <  2:3061034/93‑1 (MQ=255)
                                                            tttttACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCt                                  <  1:3163125/93‑1 (MQ=255)
                                                             ttttACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTg                                 <  1:1528391/93‑1 (MQ=255)
                                                              tttACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTgg                                >  2:158653/1‑93 (MQ=255)
                                                               ttACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGt                               >  1:1651601/1‑93 (MQ=255)
                                                               ttACGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGAGGGTCGGAGTGCGGGTATTCCTATGTGCACCGCGCTGGt                               <  2:2092004/93‑1 (MQ=255)
                                                                 aCGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGttt                             <  2:8281/93‑1 (MQ=255)
                                                                 aCGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGttt                             >  1:125797/1‑93 (MQ=255)
                                                                  cGAAAAACGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTc                            <  1:6430794/93‑1 (MQ=255)
                                                                    aaaaaCGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCtt                          >  1:2509479/1‑93 (MQ=255)
                                                                    aaaaaCGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCtt                          <  1:4743375/93‑1 (MQ=255)
                                                                     aaaaCGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTg                         <  2:5923207/93‑1 (MQ=255)
                                                                     aaaaCGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTg                         >  2:3109462/1‑93 (MQ=255)
                                                                      aaaCGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGc                        >  1:3481770/1‑93 (MQ=255)
                                                                       aaCGTGACATCTACGTCCACGCCCCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGcc                       <  1:421132/93‑1 (MQ=255)
                                                                       aaCGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGcc                       >  1:2672582/1‑93 (MQ=255)
                                                                         cGTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTg                     <  2:5625072/93‑1 (MQ=255)
                                                                          gTGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGa                    <  1:4479205/93‑1 (MQ=255)
                                                                           tGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGAc                   <  1:5346551/93‑1 (MQ=255)
                                                                           tGACATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGAc                   >  1:6354250/1‑93 (MQ=255)
                                                                             aCATCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCTGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCg                 <  2:4790299/93‑1 (MQ=255)
                                                                               aTCCACGTCCACGCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGt               >  1:374726/1‑93 (MQ=255)
                                                                                tccacgtccacgCACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTa              >  1:6563839/1‑93 (MQ=255)
                                                                                   acgtccacgcACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAAcc           >  2:3849455/1‑93 (MQ=255)
                                                                                     gtccacgcACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCg         <  2:4547940/93‑1 (MQ=255)
                                                                                     gtccacgcACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCg         >  2:6012392/1‑93 (MQ=255)
                                                                                      tccacgtACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCgg        >  1:1591822/1‑93 (MQ=255)
                                                                                        cacgtACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGtt      >  1:4901761/1‑93 (MQ=255)
                                                                                        cacgcACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGtt      >  2:4840668/1‑93 (MQ=255)
                                                                                        cacgcACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGtt      <  1:4928216/93‑1 (MQ=255)
                                                                                        cacgcACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGtt      >  2:5734325/1‑93 (MQ=255)
                                                                                        cacgcACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGtt      >  2:1814147/1‑93 (MQ=255)
                                                                                        cacgcACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGtt      >  2:2765002/1‑93 (MQ=255)
                                                                                         acgcACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTc     >  2:6355342/1‑93 (MQ=255)
                                                                                            |                                                                                            
CATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGAAAGATGGTCCGAGTGCCGGTATTGCTATGTGCACCGCGCTGGTTTCTTGCCTGACCGGTAACCCGGTTCGTG  >  NC_000913/460022‑460206

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