Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,791,400 T→C Q61Q (CAA→CAG btuD ← ATP‑binding component of vitamin B12 transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,791,4000TC99.0% 302.2 / ‑4.7 100Q61Q (CAA→CAGbtuDATP‑binding component of vitamin B12 transport system
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (55/44);  total (55/45)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 101 total aligned reads displayed.

CAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCAC  >  NC_000913/1791311‑1791489
                                                                                         |                                                                                         
cAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGc                                                                                        <  2:3647743/93‑1 (MQ=255)
cAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGc                                                                                        >  1:3265739/1‑93 (MQ=255)
 aGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGcc                                                                                       >  2:4230092/1‑93 (MQ=255)
  gACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGccc                                                                                      >  1:2905269/1‑93 (MQ=255)
     cGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCcgc                                                                                   <  2:670314/93‑1 (MQ=255)
        cGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGaa                                                                                <  2:2609751/93‑1 (MQ=255)
         gTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGGGAAc                                                                               >  1:2710041/1‑93 (MQ=255)
         gTTGCAAACGGCGGCGTCTGCTGTGGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGagc                                                                               <  2:2973787/93‑3 (MQ=255)
          ttGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTTCGGACCATGCTTCCAGTGGCTGCCCCGCGAACt                                                                              >  1:4580072/1‑93 (MQ=255)
          ttGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACt                                                                              >  2:2178965/1‑93 (MQ=255)
           tGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTg                                                                             >  1:3089440/1‑93 (MQ=255)
            gCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGa                                                                            <  1:4610166/93‑1 (MQ=255)
             cAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGaa                                                                           <  1:3638738/93‑1 (MQ=255)
             cAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGaa                                                                           >  1:2956464/1‑93 (MQ=255)
             cAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGaa                                                                           <  1:552829/93‑1 (MQ=255)
               aaCGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATg                                                                         >  1:3076424/1‑93 (MQ=255)
                 cggcggCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCt                                                                       <  1:1941648/93‑1 (MQ=255)
                  ggcggcGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCtt                                                                      >  2:3174470/1‑93 (MQ=255)
                  ggcggcGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCtt                                                                      <  2:4655362/93‑1 (MQ=255)
                     ggcgTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTccc                                                                   >  2:3108921/1‑93 (MQ=255)
                     ggcgTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTccc                                                                   <  2:2316700/93‑1 (MQ=255)
                        gTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTa                                                                >  1:307420/1‑93 (MQ=255)
                         tCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTAc                                                               <  1:330943/93‑1 (MQ=255)
                         tCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTAc                                                               >  1:3369542/1‑93 (MQ=255)
                         tCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTAc                                                               >  1:2394384/1‑93 (MQ=255)
                         tCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTAc                                                               >  2:345109/1‑93 (MQ=255)
                          ctgctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGGCCCGCGAACTGAATGCTTCCCTTAcc                                                              >  2:1144691/1‑93 (MQ=255)
                          ctgctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTAcc                                                              >  2:3389257/1‑93 (MQ=255)
                          ctgctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTAcc                                                              <  2:2811475/93‑1 (MQ=255)
                          ctgctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTAcc                                                              <  1:3752032/93‑1 (MQ=255)
                           tgctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCg                                                             >  1:2786853/1‑93 (MQ=255)
                            tctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGc                                                            <  1:3644412/92‑1 (MQ=255)
                             ctgTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCt                                                           <  2:1126827/93‑1 (MQ=255)
                              tgttgtGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTg                                                          >  2:1284068/1‑93 (MQ=255)
                               gttgtGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTgg                                                         <  2:108091/93‑1 (MQ=255)
                                  gtgAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCa                                                      <  1:1942721/93‑1 (MQ=255)
                                  gtgAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCa                                                      <  2:2860219/93‑1 (MQ=255)
                                   tgAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCAt                                                     >  2:1396876/1‑93 (MQ=255)
                                    gAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGCACTGAATGCTTCCCTTACCGCTGGTCAtt                                                    >  1:1041622/1‑93 (MQ=255)
                                    gAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCAtt                                                    >  2:163979/1‑93 (MQ=255)
                                    gAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTACCTTACCGCTGGTCAtt                                                    >  2:1415245/1‑93 (MQ=255)
                                     aaaGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTTGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTc                                                   <  1:3294160/93‑1 (MQ=255)
                                     aaaGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTc                                                   <  1:1931612/93‑1 (MQ=255)
                                      aaGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTcc                                                  <  1:3359931/93‑1 (MQ=255)
                                      aaGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTcc                                                  <  1:3218947/93‑1 (MQ=255)
                                      aaGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTcc                                                  >  2:406316/1‑93 (MQ=255)
                                      aaGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTcc                                                  <  1:2941712/93‑1 (MQ=255)
                                       aGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCg                                                 >  2:4423379/1‑93 (MQ=255)
                                       aGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCg                                                 >  1:2260199/1‑93 (MQ=255)
                                       aGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCCCGAACTGAATGCTTCCCTTACCGCTGGCCATTCCg                                                 >  1:2260200/1‑93 (MQ=255)
                                         aTAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGc                                               <  2:2762118/93‑1 (MQ=255)
                                           aGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGTCTGCCCCGCGAACTGAATGCTTCCCTTTCCGCTGGTCATTCCGGCCa                                             <  2:4128577/93‑1 (MQ=255)
                                             gcgcGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCAtt                                           <  1:2747488/93‑1 (MQ=255)
                                              cgcgATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTc                                          <  2:2657009/93‑1 (MQ=255)
                                              cgcgATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTc                                          <  1:3168381/93‑1 (MQ=255)
                                                 gATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTcgcg                                       >  2:3031148/1‑93 (MQ=255)
                                                    gCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCa                                    >  1:4476887/1‑93 (MQ=255)
                                                      aGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGt                                  >  2:3964701/1‑93 (MQ=255)
                                                         gcgAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTTCCTTACCGCTGGTCATTCCGGCCATCGGCGGCAGTAAg                               >  2:3341194/1‑93 (MQ=255)
                                                         gcgAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAg                               >  2:2021718/1‑93 (MQ=255)
                                                          cgAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAgg                              <  1:286723/93‑1 (MQ=255)
                                                           gAGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAgga                             >  2:2347038/1‑92 (MQ=255)
                                                            aGTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTa                            <  2:4017613/93‑1 (MQ=255)
                                                             gTTTTGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTAc                           <  1:3061950/93‑1 (MQ=255)
                                                              ttttGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAATGTAct                          >  1:46062/1‑93 (MQ=255)
                                                               tttGTTGCGGAGCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTActc                         <  2:3865834/93‑1 (MQ=255)
                                                               tttGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTActc                         <  1:1511532/93‑1 (MQ=255)
                                                               tttGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTActc                         <  2:3263439/93‑1 (MQ=255)
                                                               tttGTTGCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTActc                         >  2:1210185/1‑93 (MQ=255)
                                                                ttgttgCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTActct                        >  1:756350/1‑93 (MQ=255)
                                                                ttgttgCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTActct                        <  2:1581378/93‑1 (MQ=255)
                                                                ttgttgCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTActct                        <  1:1357608/93‑1 (MQ=255)
                                                                 tgttgCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCtt                       >  2:145/1‑93 (MQ=255)
                                                                  gttgCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTa                      <  1:2463807/93‑1 (MQ=255)
                                                                  gttgCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTa                      >  2:117088/1‑93 (MQ=255)
                                                                   ttgCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTTGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTAc                     >  1:2058770/1‑93 (MQ=255)
                                                                   ttgCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTAc                     <  2:3784693/93‑1 (MQ=255)
                                                                   ttgCGGACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTAc                     >  1:525400/1‑93 (MQ=255)
                                                                     gCGTACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTAccc                   <  2:2272037/93‑1 (MQ=255)
                                                                       ggACCATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCcgc                 >  1:26295/1‑93 (MQ=255)
                                                                          ccATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGcc              >  1:3739223/1‑93 (MQ=255)
                                                                          ccATGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGcc              >  2:1624935/1‑93 (MQ=255)
                                                                            aTGCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCAt            >  2:1227253/1‑93 (MQ=255)
                                                                              gCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTc          <  1:1253794/93‑1 (MQ=255)
                                                                              gCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTc          <  1:1788115/93‑1 (MQ=255)
                                                                              gCTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTc          >  1:1812011/1‑93 (MQ=255)
                                                                               cTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCg         <  2:3083386/93‑1 (MQ=255)
                                                                               cTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCg         <  2:1259257/93‑1 (MQ=255)
                                                                               cTTCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCg         >  1:2186581/1‑93 (MQ=255)
                                                                                ttCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCgg        >  1:2576035/1‑93 (MQ=255)
                                                                                 tCCAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGc       <  2:4238570/93‑1 (MQ=255)
                                                                                  ccAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGcc      >  1:1713125/1‑93 (MQ=255)
                                                                                  ccAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGcc      >  1:2824017/1‑93 (MQ=255)
                                                                                   cAGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGccc     >  1:649165/1‑93 (MQ=255)
                                                                                    aGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGcccc    >  1:1697665/1‑93 (MQ=255)
                                                                                    aGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGcccc    >  1:1939227/1‑93 (MQ=255)
                                                                                    aGTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGcccc    >  1:2237850/1‑93 (MQ=255)
                                                                                     gTGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCcca   >  1:3134503/1‑93 (MQ=255)
                                                                                      tGGTTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCccac  <  1:212006/93‑1 (MQ=255)
                                                                                      tGGCTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGGCAGTAAGGTACTCTTACCCGCGCCATTCGGCCccac  >  1:3109355/1‑93 (MQ=255)
                                                                                         |                                                                                         
CAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCAC  >  NC_000913/1791311‑1791489

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