Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,228 C→A 17.1% intergenic (+23/+29) 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,2280CA17.1% 232.3 / 36.0 88intergenic (+23/+29)yfeD/gltXDUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (41/32);  new base A (12/3);  total (53/35)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

AAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTA  >  NC_000913/2519129‑2519326
                                                                                                   |                                                                                                  
aaaaaGTCGGGCACTTTGTGTCATGTGAAGGCAATCGCAATTACCTCGCCCGTATCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtc                                                                                                    <  1:1017608/100‑1 (MQ=255)
    aGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCt                                                                                                >  1:1621147/1‑100 (MQ=255)
      tCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGc                                                                                              <  1:483591/100‑1 (MQ=255)
        gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                            <  1:747473/100‑1 (MQ=255)
        gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                            >  2:233689/1‑100 (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:2081775/100‑1 (MQ=255)
                                            cTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcg                                                        >  2:2290331/1‑100 (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                                             >  2:248390/1‑100 (MQ=255)
                                                       cGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAAt                                             >  1:2377965/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:1490186/100‑1 (MQ=255)
                                                                       gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttg                             >  2:425166/1‑100 (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:1928473/1‑100 (MQ=39)
                                                                                   aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg                 >  2:326951/1‑100 (MQ=39)
                                                                                   aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg                 >  2:2492366/1‑100 (MQ=39)
                                                                                     cGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt               >  2:318393/1‑100 (MQ=39)
                                                                                     cGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt               >  1:1762274/1‑100 (MQ=39)
                                                                                      ggCAGGAGATAATATCCTGCCGTTTATCTTTTTATTTTTTATTTACTGCTGATTTTAACGTTCTGCAATAAAATCCTGCGCTTTGTTTATATGCTCGAt               >  1:2341854/1‑99 (MQ=16)
                                                                                       gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATAc             >  2:2044183/1‑100 (MQ=255)
                                                                                       gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATAc             >  1:1057410/1‑100 (MQ=255)
                                                                                       gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTAGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATAc             >  1:1290383/1‑100 (MQ=255)
                                                                                           gATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCg         >  2:992450/1‑100 (MQ=255)
                                                                                               ttatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTc     <  1:1769317/100‑1 (MQ=255)
                                                                                                tatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCt    <  2:2188415/100‑1 (MQ=255)
                                                                                                 atCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCtt   >  2:2001952/1‑100 (MQ=255)
                                                                                                  tctcCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTa  <  2:1065555/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
AAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTA  >  NC_000913/2519129‑2519326

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: