Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,519,341 A→G noncoding (69/76 nt) lysV → tRNA‑Lys

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,519,3410AG99.1% 549.5 / ‑5.6 221noncoding (69/76 nt)lysVtRNA‑Lys
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (121/96);  total (122/99)
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 229 total aligned reads displayed.

TCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCG  >  NC_000913/2519251‑2519430
                                                                                           |                                                                                         
tCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGcc                                                                                         >  2:9937997/1‑93 (MQ=37)
 cccGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGccc                                                                                        >  2:2239781/1‑93 (MQ=37)
  ccGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCcca                                                                                       >  2:4755601/1‑93 (MQ=37)
    gTCATCACCCACCACTTCGGGTCGTTAGCGCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCccacc                                                                                     >  2:6222988/1‑93 (MQ=255)
     tcatcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCccacca                                                                                    >  1:9429406/1‑93 (MQ=37)
     tcatcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCccacca                                                                                    >  2:4737691/1‑93 (MQ=37)
     tcatcaCCCACCACTTCGGGTCGCCAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCccacca                                                                                    >  2:8671451/1‑93 (MQ=255)
      catcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCaa                                                                                   <  2:874325/93‑1 (MQ=255)
      catcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCaa                                                                                   >  1:4937737/1‑93 (MQ=255)
       atcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAAt                                                                                  <  2:3082028/93‑1 (MQ=255)
       atcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAAt                                                                                  >  1:5603518/1‑93 (MQ=255)
       atcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAAt                                                                                  <  1:5272855/93‑1 (MQ=255)
       atcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAAt                                                                                  <  2:9126510/93‑1 (MQ=255)
        tcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATg                                                                                 >  1:6830628/1‑93 (MQ=255)
        tcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATg                                                                                 >  1:6830630/1‑93 (MQ=255)
           cccaccACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTaa                                                                              <  2:1124412/93‑1 (MQ=255)
           cccaccACTTCGGGTCGTCAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTaa                                                                              <  1:2387765/93‑1 (MQ=255)
            ccaccaCTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTaaa                                                                             >  1:3185352/1‑93 (MQ=255)
            ccaccaCTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTaaa                                                                             <  1:8537215/93‑1 (MQ=255)
            ccaccaCTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTaaa                                                                             <  2:3363433/93‑1 (MQ=255)
              accacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTacaaa                                                                           >  2:5091934/1‑93 (MQ=255)
                cacTTCGGTTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAg                                                                         <  2:5925005/93‑1 (MQ=255)
                cacTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAg                                                                         <  2:1748915/93‑1 (MQ=255)
                 acTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAgc                                                                        >  1:938611/1‑93 (MQ=25)
                 acTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAgc                                                                        <  1:8613844/93‑1 (MQ=25)
                 acTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAga                                                                        >  2:2174295/1‑92 (MQ=25)
                  cTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAgcg                                                                       <  1:5890934/93‑1 (MQ=21)
                  cTTCGGGTCGTCAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAgcg                                                                       >  2:4038266/1‑93 (MQ=255)
                    tCGGGTTGTTAGCTCAGTTGGTAAGGCCGGTTGATTTTAAACAAATGGTCCGAGGTTCCAAACCTGCACGGCCCACCAATGTAAAAAAGCGcc                                                                     <  2:1317460/93‑1 (MQ=2)
                     cGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGccc                                                                    <  2:5219238/93‑1 (MQ=14)
                     cGGGTCGTCAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGccc                                                                    <  1:1143760/93‑1 (MQ=9)
                                     ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGTTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTa                                                    <  1:5038055/93‑1 (MQ=9)
                                      tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCTAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTAc                                                   >  1:4179176/1‑93 (MQ=12)
                                        gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTa                                                 <  1:3917725/93‑1 (MQ=18)
                                         tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTAt                                                <  2:7489428/93‑1 (MQ=18)
                                         tAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTAt                                                <  2:2908198/93‑1 (MQ=18)
                                         tAGAGCAGTGGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTAt                                                <  1:4793295/93‑1 (MQ=16)
                                           agagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATc                                               <  2:328405/93‑1 (MQ=18)
                                           agagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATc                                               <  1:8374063/93‑1 (MQ=18)
                                             agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAACAGCGCCCTAAAGGGGCTTTTTTACTATCTg                                             >  2:2543955/1‑93 (MQ=11)
                                             agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTg                                             <  2:6909205/93‑1 (MQ=21)
                                              gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGc                                            >  1:4724061/1‑93 (MQ=18)
                                                aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGa                                          >  1:6487152/1‑93 (MQ=18)
                                                aGTTGACTTTGAATCAATTGGACGCAGGTACGAATCCTGCACGTCCCACCAAGGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGa                                          <  1:7604227/93‑1 (MQ=11)
                                                 gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGGCCTAAAGGCGCTTGTTTACTATCTGCGAt                                         >  1:9002767/1‑93 (MQ=255)
                                                 gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGAt                                         >  2:3798569/1‑93 (MQ=18)
                                                 gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGAt                                         >  1:1324188/1‑93 (MQ=18)
                                                 gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGAt                                         >  1:6658078/1‑93 (MQ=18)
                                                  ttGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATa                                        >  1:797068/1‑93 (MQ=18)
                                                   tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATAc                                       >  1:1058360/1‑93 (MQ=18)
                                                   tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATAc                                       <  1:986089/93‑1 (MQ=18)
                                                    gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACt                                      <  2:8494039/93‑1 (MQ=21)
                                                     aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAACAAAGCGCCCTAAAGGAGCTTTTTTACTATCTGCGATCCTc                                     >  1:1520431/1‑93 (MQ=11)
                                                      cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCa                                    <  1:7081399/93‑1 (MQ=18)
                                                       ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGGCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaa                                   <  1:9130124/93‑1 (MQ=255)
                                                        tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaaa                                  >  1:1898895/1‑93 (MQ=18)
                                                        tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaaa                                  >  2:9534562/1‑93 (MQ=18)
                                                        tttAATCAATTGGTCGCAGGTGCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaaa                                  <  2:7939763/93‑1 (MQ=16)
                                                        tttAATCAATGGGTAGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaaa                                  <  2:1464613/93‑1 (MQ=255)
                                                          tAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGa                                <  1:8095748/93‑1 (MQ=21)
                                                           aaTCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAt                               <  2:6255305/93‑1 (MQ=21)
                                                           aaTCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAt                               <  1:5950054/93‑1 (MQ=18)
                                                           aaTCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAt                               >  2:6355752/1‑93 (MQ=18)
                                                            aTCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGAtt                              >  1:1765634/1‑93 (MQ=18)
                                                              cAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCg                            >  1:2334720/1‑93 (MQ=18)
                                                              cAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCg                            >  1:6202531/1‑93 (MQ=18)
                                                               aaTTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGa                           <  2:10014472/93‑1 (MQ=21)
                                                               aaTTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGa                           >  2:5858213/1‑93 (MQ=18)
                                                                aTTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGaa                          >  2:2611144/1‑93 (MQ=18)
                                                                 ttGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAAc                         >  1:6256281/1‑93 (MQ=18)
                                                                 ttGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAAc                         <  2:1673249/93‑1 (MQ=18)
                                                                  tGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAAcc                        <  2:7987094/93‑1 (MQ=21)
                                                                   ggTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCt                       <  2:3623092/93‑1 (MQ=18)
                                                                      cGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgca                    >  2:6014044/1‑93 (MQ=18)
                                                                        cAGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgcacc                  >  1:4430159/1‑91 (MQ=16)
                                                                         aGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgcagca                 <  1:2469015/93‑1 (MQ=18)
                                                                         aGGTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgcagca                 <  2:9091637/93‑1 (MQ=18)
                                                                         aGGTTCGAATCCTGCACGCCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgcagca                 <  2:8677863/93‑1 (MQ=18)
                                                                          ggTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGcagcag                <  2:677202/93‑1 (MQ=18)
                                                                          ggTTCGAATCCGGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGcagcag                <  1:3717047/93‑1 (MQ=16)
                                                                           gTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAgg               >  1:7665542/1‑93 (MQ=18)
                                                                           gTTCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAgg               >  2:3472466/1‑93 (MQ=18)
                                                                            ttCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGt              >  2:6195287/1‑93 (MQ=18)
                                                                             tCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGtt             <  1:7240979/93‑1 (MQ=18)
                                                                             tCGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGtt             <  1:4860128/93‑1 (MQ=18)
                                                                              cGAATCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGttt            >  2:8662976/1‑93 (MQ=18)
                                                                                aaTCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgg          >  2:7265788/1‑92 (MQ=16)
                                                                                aaTCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga          <  1:3475691/93‑1 (MQ=21)
                                                                                 aTCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag         <  1:460073/93‑1 (MQ=18)
                                                                                 aTCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag         >  1:2045554/1‑93 (MQ=18)
                                                                                  tCCTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagt        >  1:4313522/1‑93 (MQ=18)
                                                                                   ccTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagtt       <  1:2697119/93‑1 (MQ=18)
                                                                                   ccTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagtt       >  2:7954666/1‑93 (MQ=18)
                                                                                    cTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTa                                                 >  2:6683563/1‑50 (MQ=255)
                                                                                    cTGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttg      >  1:7367548/1‑93 (MQ=18)
                                                                                     tGCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga     <  2:2003382/93‑1 (MQ=18)
                                                                                     tGCACGGCACACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga     <  2:9646331/93‑1 (MQ=16)
                                                                                      gCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag    <  2:7648966/93‑1 (MQ=18)
                                                                                      gCACGGCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag    <  1:6951251/93‑1 (MQ=18)
                                                                                        aCGGCCAACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcg  >  2:7204169/1‑93 (MQ=16)
                                                                                           |                                                                                         
TCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCG  >  NC_000913/2519251‑2519430

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