Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,222 T→G 17.7% intergenic (+17/+35) yfeD → / ← gltX DUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,519,2220TG17.7% 197.5 / 41.3 86intergenic (+17/+35)yfeD/gltXDUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (40/30);  new base G (12/3);  total (52/34)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.90e-01

AGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCG  >  NC_000913/2519123‑2519319
                                                                                                   |                                                                                                 
aGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGat                                                                                                   >  2:720064/1‑100 (MQ=255)
 gCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGata                                                                                                  >  2:747476/1‑100 (MQ=255)
 gCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGata                                                                                                  <  2:1244559/100‑1 (MQ=255)
  cAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGatat                                                                                                 <  2:1240939/100‑1 (MQ=255)
    ggAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatta                                                                                               <  2:1950484/100‑1 (MQ=255)
    ggAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatta                                                                                               <  1:1056634/100‑1 (MQ=255)
    ggAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatta                                                                                               >  2:1980534/1‑100 (MQ=255)
    ggAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatta                                                                                               <  2:2290475/100‑1 (MQ=255)
      aaaaaGTCGGGCACTTTGTGTCATGTGAAGGCAATCGCAATTACCTCGCCCGTATCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtc                                                                                             <  1:1017608/100‑1 (MQ=255)
      aaaaaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtattat                                                                                              <  1:8835/99‑1 (MQ=255)
          aGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCt                                                                                         >  1:1621147/1‑100 (MQ=255)
            tCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGc                                                                                       <  1:483591/100‑1 (MQ=255)
              gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                     >  2:233689/1‑100 (MQ=255)
              gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                     <  1:747473/100‑1 (MQ=255)
               cgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGt                                                                                    >  2:1769326/1‑100 (MQ=255)
                gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                   <  1:966290/100‑1 (MQ=255)
                  aaTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtttt                                                                                 <  2:531611/100‑2 (MQ=38)
                   aTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttt                                                                                <  2:1068736/100‑3 (MQ=38)
                       ttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttt                                                                            <  1:1754695/100‑1 (MQ=255)
                        tGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCtttt                                                                           <  2:2574385/100‑1 (MQ=255)
                         gTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttttt                                                                          >  1:2112270/1‑100 (MQ=255)
                           cccGTGACGGAATTCGCAATTTCCTCGCCCGTCTCGAGATTGACGAATCAGCATAAATAAAAAACGGCGGGATATTATCTCCTGCCGTTAATCTTTTTa                                                                         >  1:2414416/3‑99 (MQ=255)
                           cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTac                                                                        <  2:1268473/100‑1 (MQ=255)
                           cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttttatttatg                                                                        <  2:147630/100‑11 (MQ=25)
                             cGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTacac                                                                      >  2:102769/1‑100 (MQ=255)
                              gTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACg                                                                     >  2:1396756/1‑100 (MQ=255)
                              gTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGAGATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTacac                                                                      >  1:27279/1‑99 (MQ=255)
                                  aGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGACTCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTaa                                                                 >  1:863514/1‑100 (MQ=255)
                                  aGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTaa                                                                 <  1:1167298/100‑1 (MQ=255)
                                   ggAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTaa                                                                 <  1:1056984/99‑1 (MQ=255)
                                     aaTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTa                                                              <  1:2595089/100‑1 (MQ=255)
                                       ttCGCAATTTCCTCGCTCGTCTCGATATTGACGAATCAGCATAAATAAAAATCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAct                                                            <  1:2116842/100‑1 (MQ=255)
                                       ttCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAct                                                            >  2:760304/1‑100 (MQ=255)
                                        tCGCAATTTCTTCTCCCGTCTCGATATTGACTAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActg                                                           <  1:1371451/100‑1 (MQ=255)
                                         cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                          >  2:888937/1‑100 (MQ=255)
                                         cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                          <  1:233687/100‑1 (MQ=255)
                                          gCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgct                                                         <  2:1689618/100‑1 (MQ=255)
                                            aaTTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGa                                                       <  1:239574/100‑1 (MQ=255)
                                                  cTCGCCCGTCTCGATATTGACGAATTAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcg                                                 <  2:2064226/100‑1 (MQ=255)
                                                  cTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcg                                                 >  2:2290331/1‑100 (MQ=255)
                                                  cTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcg                                                 <  2:2081775/100‑1 (MQ=255)
                                                    tccgatcttcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgcg                                               <  2:1156357/92‑1 (MQ=255)
                                                    cGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgcg                                               <  1:809346/100‑1 (MQ=255)
                                                     gCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGt                                              <  2:1145383/100‑1 (MQ=255)
                                                      cccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGtt                                             >  2:1461433/1‑100 (MQ=255)
                                                      cccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGtt                                             >  1:1032037/1‑100 (MQ=255)
                                                      cccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGtt                                             >  1:1032407/1‑100 (MQ=255)
                                                       ccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTc                                            <  1:966394/100‑1 (MQ=255)
                                                        cGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGt                                              <  1:157621/97‑1 (MQ=255)
                                                        cGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCa                                           >  2:1195342/1‑100 (MQ=255)
                                                           ctcGATATTGACCAATCAGCATATCTAAAAAACGGTCTGATATTATCTCCTGCCGTTCATCTTTTTACTCGCTAATTACTGCTGATTTTCGCGCTCAg                                          <  1:839967/98‑1 (MQ=255)
                                                            tcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgcgagatcgga                                       >  1:1156353/1‑92 (MQ=255)
                                                             cGATCTTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGCTTTTCGCTTTCAGCAAt                                      >  1:2420394/1‑100 (MQ=255)
                                                             cGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAAt                                      >  1:2377965/1‑100 (MQ=255)
                                                             cGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAAt                                      >  2:248390/1‑100 (MQ=255)
                                                             cGATATTGAAGAATCAGCATAACTACAAAACGGCCGGACATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCa                                        <  2:116982/98‑1 (MQ=255)
                                                              gATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATa                                     >  1:577764/1‑100 (MQ=255)
                                                                tatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCACTaaa                                   >  1:755273/1‑100 (MQ=255)
                                                                tatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaa                                   >  2:897243/1‑100 (MQ=255)
                                                                 cttGACGAATCAGGATAAATAAAAAACGGAAGGATATTAGCTCGTGACGTTTAGCTTTTTACACGCTCATTACTGCTGATTTTCTCTTTCAgaacta                                     >  1:1726271/2‑92 (MQ=255)
                                                                   tGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATc                                >  1:2081665/1‑100 (MQ=255)
                                                                   tGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATc                                >  1:1990948/1‑100 (MQ=255)
                                                                      cGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                             <  2:500477/100‑1 (MQ=255)
                                                                       gAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgc                            >  2:768057/1‑100 (MQ=255)
                                                                        aaTCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgcg                           >  2:456894/1‑100 (MQ=255)
                                                                         aTCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgcgc                          <  2:1533734/100‑1 (MQ=255)
                                                                             gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttg                      >  2:425166/1‑100 (MQ=255)
                                                                             gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttg                      <  2:1490186/100‑1 (MQ=255)
                                                                              cATAAATAAAAAAGGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                     >  1:1602394/1‑100 (MQ=255)
                                                                              cATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                     <  1:1869297/100‑1 (MQ=255)
                                                                              cATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                     >  2:857861/1‑100 (MQ=255)
                                                                               taaaaagataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                    >  2:189293/10‑100 (MQ=25)
                                                                                taaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACTCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgttg                   <  1:532217/100‑1 (MQ=255)
                                                                                taaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgttg                   >  2:542511/1‑100 (MQ=255)
                                                                                aaaaagataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgttg                   >  2:1080437/9‑100 (MQ=25)
                                                                                 aaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGa                  >  2:1258371/1‑100 (MQ=255)
                                                                                 aaaagataaaCGGCAGGAGATAATATCCTGCCGATTATCTTTTTACACGCTAATTACTGCTGATTTTGACGTTTAGCAATAAAATTCAGCGCTTTGTTGa                  >  1:2405137/8‑100 (MQ=21)
                                                                                  aataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAt                 <  1:1461425/100‑1 (MQ=255)
                                                                                  aaagataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAt                 >  2:2542581/7‑100 (MQ=37)
                                                                                   aagataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                >  2:2387758/6‑100 (MQ=37)
                                                                                    taaaAAACGGCAGGATATTATCTCGTGCCGTTTATCTTTTAAACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTGGTTGAt                 >  1:2559016/1‑97 (MQ=255)
                                                                                       taaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCCATAAAATCCAGCGCTTTGTTGATACGCt            >  1:1258886/2‑100 (MQ=38)
                                                                                       aaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCt            >  2:2523246/1‑100 (MQ=255)
                                                                                        aaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGGTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTc           <  1:1816361/100‑1 (MQ=255)
                                                                                         aaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg          >  1:2335850/1‑100 (MQ=255)
                                                                                         aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg          >  2:2492366/1‑100 (MQ=39)
                                                                                         aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg          >  2:1928473/1‑100 (MQ=39)
                                                                                         aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg          >  2:326951/1‑100 (MQ=39)
                                                                                           cGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt        >  1:1762274/1‑100 (MQ=39)
                                                                                           cGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt        >  2:318393/1‑100 (MQ=39)
                                                                                            ggCAGGAGATAATATCCTGCCGTTTATCTTTTTATTTTTTATTTACTGCTGATTTTAACGTTCTGCAATAAAATCCTGCGCTTTGTTTATATGCTCGAt        >  1:2341854/1‑99 (MQ=16)
                                                                                             gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATAc      >  1:1057410/1‑100 (MQ=255)
                                                                                             gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATAc      >  2:2044183/1‑100 (MQ=255)
                                                                                             gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTAGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATAc      >  1:1290383/1‑100 (MQ=255)
                                                                                                 gATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCg  >  2:992450/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
AGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCG  >  NC_000913/2519123‑2519319

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