Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 780,918 T→C intergenic (+43/‑390) lysQ → / → nadA tRNA‑Lys/quinolinate synthetase A protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913780,9180TC99.3% 355.1 / ‑4.8 148intergenic (+43/‑390)lysQ/nadAtRNA‑Lys/quinolinate synthetase A protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (82/65);  total (82/66)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

Only 100 of 152 total aligned reads displayed.

TGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGC  >  NC_000913/780828‑780999
                                                                                          |                                                                                 
tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGacac                                                                                 <  1:4689229/93‑1 (MQ=21)
 gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACt                                                                                >  1:200529/1‑93 (MQ=18)
 gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACt                                                                                >  2:4453487/1‑93 (MQ=18)
   cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAACAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCa                                                                              >  1:3133998/1‑93 (MQ=16)
   cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCa                                                                              >  1:3133997/1‑93 (MQ=18)
   cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCa                                                                              <  1:1919974/93‑1 (MQ=21)
    ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCaa                                                                             >  1:3333791/1‑93 (MQ=18)
    ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCaa                                                                             >  1:4712266/1‑93 (MQ=18)
     tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCaaa                                                                            <  1:3955483/93‑1 (MQ=18)
     tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCaaa                                                                            >  1:2690213/1‑93 (MQ=18)
       tAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGa                                                                          <  1:2243934/93‑1 (MQ=21)
        aaTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGAt                                                                         <  2:456719/93‑1 (MQ=21)
        aaTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGAt                                                                         <  1:610865/93‑1 (MQ=18)
         aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGAtt                                                                        <  2:4420547/93‑1 (MQ=21)
         aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGAtt                                                                        <  2:3011705/93‑1 (MQ=21)
          tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTc                                                                       >  1:4327904/1‑93 (MQ=18)
           cAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCg                                                                      <  1:1268764/93‑1 (MQ=21)
           cAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCg                                                                      >  2:1134499/1‑93 (MQ=18)
            aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGa                                                                     >  2:2194012/1‑93 (MQ=18)
             aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGaa                                                                    >  2:3411125/1‑93 (MQ=18)
             aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGaa                                                                    <  2:489155/93‑1 (MQ=21)
              ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAAc                                                                   <  2:2110585/93‑1 (MQ=21)
              ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAAc                                                                   >  1:726433/1‑93 (MQ=18)
                ggTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCt                                                                 <  1:4138384/93‑1 (MQ=18)
                ggTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCt                                                                 >  2:4388275/1‑93 (MQ=18)
                   cGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTgca                                                              <  2:1544727/93‑1 (MQ=18)
                   cGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTgca                                                              <  1:3970416/93‑1 (MQ=18)
                    gCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTgcag                                                             >  2:4557381/1‑93 (MQ=18)
                     cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTgcagc                                                            >  1:699818/1‑93 (MQ=18)
                     cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTgcagc                                                            >  2:796877/1‑93 (MQ=18)
                      aGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTgcagca                                                           <  1:2299454/93‑1 (MQ=18)
                      aGGTGCGAATCCTGCACGACCCACCAATGTAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGcagcag                                                          <  2:3802039/93‑1 (MQ=11)
                       ggTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGcagcag                                                          <  1:4580890/93‑1 (MQ=18)
                         ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGt                                                        <  1:2183274/93‑1 (MQ=21)
                          tCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAggct                                                       >  1:1050161/1‑91 (MQ=16)
                          tCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGtt                                                       <  2:1446497/93‑1 (MQ=18)
                          tCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGtt                                                       <  1:1265187/93‑1 (MQ=21)
                          tCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGtt                                                       <  1:2525894/93‑1 (MQ=18)
                           cGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGttt                                                      <  2:4653642/93‑1 (MQ=21)
                             aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTttga                                                    <  1:1472168/93‑1 (MQ=21)
                              aTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTttgag                                                   <  1:3831232/93‑1 (MQ=18)
                              aTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTttgag                                                   >  1:1543785/1‑93 (MQ=18)
                                 cTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTttgagttg                                                >  2:1016730/1‑93 (MQ=18)
                                     aCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcg                                            >  1:4302508/1‑93 (MQ=18)
                                     aCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcg                                            <  1:4591477/93‑1 (MQ=18)
                                      cGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcgc                                           >  1:2493140/1‑93 (MQ=18)
                                       gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                          <  1:107368/93‑1 (MQ=18)
                                        aCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAg                                         >  2:4383211/1‑93 (MQ=18)
                                          ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCg                                       >  1:2796715/1‑93 (MQ=18)
                                            accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaa                                     >  1:867550/1‑93 (MQ=18)
                                            accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaa                                     >  1:2366547/1‑93 (MQ=18)
                                              caATTTAAAAAAGCGCCCTAAAGGCGTTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                   <  2:3841152/93‑1 (MQ=255)
                                              caATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                   <  2:2625055/93‑1 (MQ=18)
                                               aaTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaaca                                  >  1:4450858/1‑93 (MQ=21)
                                               aaTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaaca                                  >  2:184159/1‑93 (MQ=18)
                                                 tGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaac                                <  1:2643275/93‑1 (MQ=18)
                                                 tGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaac                                >  1:3614969/1‑93 (MQ=21)
                                                  gTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACg                               >  2:3913842/1‑93 (MQ=21)
                                                  gTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACg                               >  1:3323736/1‑93 (MQ=21)
                                                   tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACgg                              >  1:1488368/1‑93 (MQ=21)
                                                      aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGc                           <  1:1888214/93‑1 (MQ=18)
                                                      aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGc                           <  2:2878685/93‑1 (MQ=18)
                                                      aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGc                           <  2:2354269/93‑1 (MQ=18)
                                                      aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGc                           <  2:1692391/93‑1 (MQ=18)
                                                          gcgcCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCt                       <  1:3863379/93‑1 (MQ=18)
                                                           cgcCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTc                      <  1:3439427/93‑1 (MQ=18)
                                                             cccTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                    <  1:3142631/93‑1 (MQ=18)
                                                             cccTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                    >  1:204465/1‑93 (MQ=18)
                                                             cccTAAAGGCGCTTGTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                    <  1:3142633/93‑1 (MQ=16)
                                                              ccTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgcg                   <  2:3096021/93‑1 (MQ=18)
                                                              ccTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgcg                   >  2:3585032/1‑93 (MQ=18)
                                                               cTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcgg                  >  1:3079676/1‑93 (MQ=21)
                                                               cTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcgg                  <  2:416136/93‑1 (MQ=18)
                                                               cTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcgg                  <  2:3494703/93‑1 (MQ=18)
                                                                tAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgcgtc                 >  2:1490074/1‑91 (MQ=16)
                                                                  aaGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGa               >  1:291462/1‑93 (MQ=21)
                                                                   aGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc              >  1:3660173/1‑93 (MQ=21)
                                                                   aGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc              >  2:1161896/1‑93 (MQ=18)
                                                                   aGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc              >  2:1983032/1‑93 (MQ=21)
                                                                     gcgcTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg            >  1:4547269/1‑93 (MQ=18)
                                                                      cgcTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGc           >  2:3649690/1‑93 (MQ=21)
                                                                       gcTTTTTTGCTATCTGCGACGCTCAAAGATTCGAGCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGcc          <  2:3149419/93‑1 (MQ=255)
                                                                       gcTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGcc          <  2:3149420/93‑1 (MQ=18)
                                                                       gcTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGcc          >  2:192826/1‑93 (MQ=18)
                                                                       gcTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGcc          >  2:575950/1‑93 (MQ=21)
                                                                       gcTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGcc          >  2:1656281/1‑93 (MQ=18)
                                                                       gcTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGcc          >  2:2724481/1‑93 (MQ=21)
                                                                       gcTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGcc          >  1:15792/1‑93 (MQ=18)
                                                                        cTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGccc         <  2:1467959/93‑1 (MQ=18)
                                                                        cTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGccc         <  2:1711431/93‑1 (MQ=18)
                                                                        cTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGccc         >  2:3437766/1‑93 (MQ=21)
                                                                         ttttttGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCg        >  2:3829087/1‑93 (MQ=21)
                                                                         ttttttGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCg        >  2:4700036/1‑93 (MQ=21)
                                                                         ttttttGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCg        >  1:3692469/1‑93 (MQ=21)
                                                                            tttGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAg     <  2:1430833/93‑1 (MQ=18)
                                                                              tGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg   <  2:2779738/93‑1 (MQ=17)
                                                                              tGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg   >  1:4348052/1‑93 (MQ=17)
                                                                              tGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg   >  1:2669887/1‑93 (MQ=17)
                                                                               gCTATCTGCGACACTCAAAGATTCGAGCCTGCAGCAGGTTTGATTGGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgc  <  1:2865735/93‑1 (MQ=11)
                                                                               gCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgc  >  2:2373585/1‑93 (MQ=14)
                                                                                          |                                                                                 
TGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGC  >  NC_000913/780828‑780999

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