Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,337,699 T→C Q3R (CAG→CGG)  yrbF ← putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,337,6990TC98.4% 799.1 / 0.9 254Q3R (CAG→CGG) yrbFputative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base T (0/4);  new base C (144/106);  total (144/110)
Fisher's exact test for biased strand distribution p-value = 3.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

Only 100 of 259 total aligned reads displayed.

CGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACC  >  NC_000913/3337609‑3337788
                                                                                          |                                                                                         
cgcgGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGc                                                                                         >  1:1143745/1‑93 (MQ=255)
 gcgGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCt                                                                                        >  1:1082895/1‑93 (MQ=255)
 gcgGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCt                                                                                        >  1:246644/1‑93 (MQ=255)
  cgGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTc                                                                                       >  2:1696395/1‑93 (MQ=255)
   ggCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTcc                                                                                      >  1:478191/1‑93 (MQ=255)
   ggCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTcc                                                                                      >  1:1374380/1‑93 (MQ=255)
   ggCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTcc                                                                                      >  1:749325/1‑93 (MQ=255)
   ggCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTcc                                                                                      >  2:3208891/1‑93 (MQ=255)
       cGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATaa                                                                                  >  1:1193419/1‑93 (MQ=255)
         gTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAAtt                                                                                >  1:1028268/1‑93 (MQ=255)
          tCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAAtta                                                                               >  1:2722260/1‑92 (MQ=255)
          tCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTc                                                                               >  1:196252/1‑93 (MQ=255)
            aGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCAc                                                                             >  2:2184984/1‑93 (MQ=255)
            aGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCAc                                                                             >  2:2653260/1‑93 (MQ=255)
            aGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCAc                                                                             >  2:4492327/1‑93 (MQ=255)
             gggAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCAcc                                                                            >  1:1588569/1‑93 (MQ=255)
             gggAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCAcc                                                                            >  2:47693/1‑93 (MQ=255)
             gggAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCAcc                                                                            >  1:4479211/1‑93 (MQ=255)
              ggAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCAccc                                                                           >  2:35116/1‑93 (MQ=255)
              ggAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCAccc                                                                           >  1:3888927/1‑93 (MQ=255)
               gAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCt                                                                          >  1:250953/1‑93 (MQ=255)
                aaaTATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCtt                                                                         >  2:4686244/1‑93 (MQ=255)
                 aaTATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTc                                                                        >  2:4172065/1‑93 (MQ=255)
                 aaTATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTc                                                                        >  2:2429093/1‑93 (MQ=255)
                  atatTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTAACCCTTCg                                                                       >  1:1974041/1‑93 (MQ=255)
                    attacCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTc                                                                     <  1:3630544/93‑1 (MQ=255)
                     ttatCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCt                                                                    >  1:716221/1‑93 (MQ=255)
                      tatCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCtt                                                                   <  1:1535749/93‑1 (MQ=255)
                        tCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGc                                                                 >  1:4388383/1‑93 (MQ=255)
                        tCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGc                                                                 >  2:150850/1‑93 (MQ=255)
                            aGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTg                                                             >  2:4231010/1‑93 (MQ=255)
                             gATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGa                                                            >  1:3402221/1‑93 (MQ=255)
                             gATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGa                                                            >  1:2226138/1‑93 (MQ=255)
                             aaTGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGa                                                            <  2:3889221/92‑1 (MQ=255)
                              aTGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGAt                                                           <  2:1185573/93‑1 (MQ=255)
                               tGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGAtt                                                          >  2:2103634/1‑93 (MQ=255)
                                  aGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTc                                                       >  1:2819759/1‑93 (MQ=255)
                                  aGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTc                                                       >  1:3992313/1‑93 (MQ=255)
                                    cGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTa                                                     <  2:1818619/93‑1 (MQ=255)
                                     gATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTaa                                                    >  1:3678195/1‑93 (MQ=255)
                                     gATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTaa                                                    >  2:2950075/1‑93 (MQ=255)
                                     gATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTaa                                                    >  2:4378758/1‑93 (MQ=255)
                                      aTTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAg                                                   >  2:1084996/1‑93 (MQ=255)
                                      aTTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAg                                                   >  2:912010/1‑93 (MQ=255)
                                       ttGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGc                                                  >  1:3415905/1‑93 (MQ=255)
                                       ttGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGc                                                  >  1:151527/1‑93 (MQ=255)
                                       ttGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGc                                                  >  1:2154115/1‑93 (MQ=255)
                                       ttGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATGCGCCACAGACCGCGCCATAATTCCCCCTTCGTCTTGCGTTGATTTTCTAAGc                                                  <  2:4437434/93‑1 (MQ=255)
                                        tGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCCCCCTTCGTCTTGCGTTGATTTTCTAAGCa                                                 <  1:2321972/93‑1 (MQ=255)
                                         gCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCAt                                                >  1:2750712/1‑93 (MQ=255)
                                          ccACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATg                                               <  2:3946845/93‑1 (MQ=255)
                                           cACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATgg                                              <  2:1724062/93‑1 (MQ=255)
                                            aCGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGgc                                             >  2:712857/1‑93 (MQ=255)
                                             cgcgTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGgcg                                            <  2:2918573/93‑1 (MQ=255)
                                                gTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTc                                         <  2:3786407/93‑1 (MQ=255)
                                                 tAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCa                                        >  1:392449/1‑93 (MQ=255)
                                                 tAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCa                                        >  2:2514475/1‑93 (MQ=255)
                                                  aaaaCTGTCAGCGGACATATCGACTAAATTCGCCACAGACTGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCaa                                       <  1:3862785/93‑1 (MQ=255)
                                                  aaaaCTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCaa                                       <  2:757854/93‑1 (MQ=255)
                                                   aaaCTGACATCGCGCATATCGACTGAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAAt                                      >  1:3724277/1‑93 (MQ=255)
                                                    aaCTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAAtt                                     <  1:3085621/93‑1 (MQ=255)
                                                      cTGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTa                                   <  2:1802083/93‑1 (MQ=255)
                                                       tGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTac                                  >  1:1237735/1‑92 (MQ=255)
                                                       tGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTaa                                  >  2:878965/1‑93 (MQ=255)
                                                       tGACATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTaa                                  <  2:100313/93‑1 (MQ=255)
                                                          cATCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCt                               <  1:1744328/93‑1 (MQ=255)
                                                            tCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTg                             >  1:3186237/1‑93 (MQ=255)
                                                            tCGCGCATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTg                             >  2:31011/1‑93 (MQ=255)
                                                              gcgcATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGaa                           >  1:323269/1‑93 (MQ=255)
                                                              gcgcATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGaa                           <  1:445278/93‑1 (MQ=255)
                                                               cgcATATCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAAc                          >  1:891726/1‑93 (MQ=255)
                                                                  atatCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCa                       >  2:196738/1‑93 (MQ=255)
                                                                   tatCGACTAAATTCGCCACAGACCGCTCCCTAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCaa                      >  1:1300956/1‑93 (MQ=255)
                                                                    atCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAAc                     <  2:559612/93‑1 (MQ=255)
                                                                    atCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAAc                     >  1:3774143/1‑93 (MQ=255)
                                                                    atCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAAc                     >  2:636881/1‑93 (MQ=255)
                                                                    atCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAAc                     <  1:3579395/93‑1 (MQ=255)
                                                                     tCGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACa                    >  1:3812036/1‑93 (MQ=255)
                                                                      cGACTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACat                   <  2:2382210/93‑1 (MQ=255)
                                                                         cTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATAtt                >  1:1799346/1‑93 (MQ=255)
                                                                         cTAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATAtt                >  2:3854964/1‑93 (MQ=255)
                                                                          tAAATTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATAttt               >  1:1912094/1‑93 (MQ=255)
                                                                            aaTTCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTAc             <  2:1341240/93‑1 (MQ=255)
                                                                              ttCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAg           <  2:278787/93‑1 (MQ=255)
                                                                              ttCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAg           <  1:4748952/93‑1 (MQ=255)
                                                                               tCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGa          <  1:583743/93‑1 (MQ=255)
                                                                               tCGCCACAGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGa          <  2:3412577/93‑1 (MQ=255)
                                                                                   cacaGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAatat      >  1:2017509/1‑93 (MQ=255)
                                                                                   cacaGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAatat      >  1:2854246/1‑93 (MQ=255)
                                                                                   cacaGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAatat      >  2:1952536/1‑93 (MQ=255)
                                                                                    acaGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATAtt     >  2:3314028/1‑93 (MQ=255)
                                                                                    acaGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATAtt     <  1:3482428/93‑1 (MQ=255)
                                                                                    acaGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATAtt     >  1:1544441/1‑93 (MQ=255)
                                                                                    acaGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATAtt     >  1:2972217/1‑93 (MQ=255)
                                                                                     caGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTa    <  2:354840/93‑1 (MQ=255)
                                                                                      aGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTAc   >  2:3621943/1‑93 (MQ=255)
                                                                                      aGACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTAc   <  1:4146486/93‑1 (MQ=255)
                                                                                       gACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACGCAACATATTTACAGAATATTAcc  >  2:2630019/1‑93 (MQ=255)
                                                                                       gACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTAcc  >  2:1569504/1‑93 (MQ=255)
                                                                                       gACCGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTAcc  >  1:3037569/1‑93 (MQ=255)
                                                                                          |                                                                                         
CGCGGCACGGTCAGGGAAATATTATCGAAGATGCAGCGATTGCCACGCGTAAAACTGACATCGCGCATATCGACTAAATTCGCCACAGACTGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACC  >  NC_000913/3337609‑3337788

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