Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,295,134 (C)7→8 coding (525/618 nt) ccmA ← ATP binding protein of heme exporter A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,295,1271.C81.1% 605.4 / 140.4 217coding (532/618 nt)ccmAATP binding protein of heme exporter A
Reads supporting (aligned to +/- strand):  ref base . (15/26);  new base C (95/81);  total (110/107)
Fisher's exact test for biased strand distribution p-value = 5.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 242 total aligned reads displayed.

AGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑T‑‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCC  >  NC_000913/2295035‑2295216
                                                                                              |                                                                                          
aGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑t                                                                                             >  1:7017571/1‑93 (MQ=255)
aGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑t                                                                                             <  2:6202006/93‑1 (MQ=255)
   aCATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑T‑‑ccc                                                                                        <  1:3215678/93‑1 (MQ=255)
      tcatGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑T‑‑cccccc                                                                                     <  1:1331476/93‑1 (MQ=255)
      tcatGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑T‑‑cccccc                                                                                     >  2:6977008/1‑93 (MQ=255)
        atGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑T‑‑cccccccc                                                                                   <  2:2745307/93‑2 (MQ=255)
           cGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑T‑‑cccccccctgc                                                                                >  2:4640033/1‑89 (MQ=255)
             gCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTc                                                                               <  2:1601341/93‑1 (MQ=255)
             gCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTc                                                                               >  2:661441/1‑93 (MQ=255)
             gCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑C‑‑CCCCCCCTGCTcc                                                                              >  2:5721842/1‑93 (MQ=255)
                cTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGt                                                                            >  1:6417096/1‑93 (MQ=255)
                cTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑C‑‑CCCCCCCTGCTCCGTa                                                                           <  2:1953407/93‑1 (MQ=255)
                 tCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTa                                                                           >  1:2121444/1‑93 (MQ=255)
                   gTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTAtg                                                                         <  2:3399685/93‑1 (MQ=255)
                   gTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTAtg                                                                         >  1:5365088/1‑93 (MQ=255)
                    tttGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTAtgc                                                                        >  1:5649901/1‑93 (MQ=255)
                    tttGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTAtgc                                                                        <  1:2818014/93‑1 (MQ=255)
                    tttGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTAtgc                                                                        >  1:3061926/1‑93 (MQ=255)
                    tttGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTAtgc                                                                        >  1:2467056/1‑93 (MQ=255)
                    tttGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTAtgc                                                                        >  1:272552/1‑93 (MQ=255)
                    tttGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑C‑‑CCCCCCCTGCTCCGTAtgct                                                                       >  2:5899335/1‑93 (MQ=255)
                       gCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTAtgctgc                                                                     <  1:1146451/93‑1 (MQ=255)
                          tCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGcc                                                                  >  1:1887293/1‑93 (MQ=255)
                          tCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGcc                                                                  >  2:2604693/1‑93 (MQ=255)
                          tCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑C‑‑CCCCCCCTGCTCCGTATGCTGCGCCa                                                                 <  2:3570770/93‑1 (MQ=255)
                           cAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑C‑‑CCCCCCCTGCTCCGTATGCTGCGCCAt                                                                <  2:5117355/93‑1 (MQ=255)
                             gTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑‑‑‑CCCCCCCTGCTCCGTATGCTGCGCCATACg                                                             >  1:2623815/1‑93 (MQ=255)
                               gAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACg                                                             >  1:3801282/1‑93 (MQ=255)
                                aaaTGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGc                                                            >  1:4346372/1‑93 (MQ=255)
                                aaaTGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGc                                                            >  1:7403697/1‑93 (MQ=255)
                                aaaGGCGGCGAATTTTCCTTTCAGCACCGTGGGGCGGCTGGTGGGTAGTCGGAACCACAACC‑‑CCCCCCCTGCTCCGTATGCTGCGCCATACGc                                                            <  2:2913545/93‑1 (MQ=255)
                                  aTGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTg                                                          <  1:6133974/93‑1 (MQ=255)
                                   tGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTgg                                                         <  1:585593/93‑1 (MQ=255)
                                   tGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTgg                                                         <  2:4220813/93‑1 (MQ=255)
                                   tGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTgg                                                         >  1:865962/1‑93 (MQ=255)
                                      ggcgAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑C‑‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCa                                                     >  1:1367133/1‑93 (MQ=255)
                                       gcgAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCa                                                     >  2:4131943/1‑93 (MQ=255)
                                       gcgAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCa                                                     >  2:3002866/1‑93 (MQ=255)
                                         gAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGa                                                   <  1:6498157/93‑1 (MQ=255)
                                          aaTTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGAc                                                  >  1:6797312/1‑93 (MQ=255)
                                          aaTTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGAc                                                  <  1:2414999/93‑1 (MQ=255)
                                          aaTTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGAc                                                  >  2:6673040/1‑93 (MQ=255)
                                           aTTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACg                                                 >  1:6395650/1‑93 (MQ=255)
                                           aTTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑C‑‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGa                                                <  1:5314580/93‑1 (MQ=255)
                                            ttttACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TCCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACg                                                 >  1:3115684/1‑93 (MQ=255)
                                            ttttACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGa                                                >  2:4805228/1‑93 (MQ=255)
                                            ttttACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGa                                                <  1:2295755/93‑1 (MQ=255)
                                            ttttACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGa                                                >  2:6014478/1‑93 (MQ=255)
                                             tttACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGAt                                               <  2:4529859/93‑1 (MQ=255)
                                             tttACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGAt                                               >  2:2436962/1‑93 (MQ=255)
                                             tttACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGAt                                               <  1:4769537/93‑1 (MQ=255)
                                             tttACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGAt                                               >  2:436989/1‑93 (MQ=255)
                                                 cTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGac                                           >  2:63712/1‑93 (MQ=255)
                                                 cTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGac                                           >  1:4536669/1‑93 (MQ=255)
                                                 cTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑C‑‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGaca                                          >  2:3748810/1‑93 (MQ=255)
                                                  tttCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑C‑‑CCCCCCCTGCTCCGTATGCTGCCCCATACGCTGGGTCAGACGATCGacac                                         <  2:6420062/93‑1 (MQ=255)
                                                   ttCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGacac                                         >  1:7129796/1‑93 (MQ=255)
                                                     cAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCTGACGATCGACACCg                                       <  2:5181640/93‑1 (MQ=255)
                                                     cAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑C‑‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGt                                      <  1:517483/93‑1 (MQ=255)
                                                       gCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGtt                                     >  2:5622110/1‑93 (MQ=255)
                                                          aCGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAAc                                  >  2:2650446/1‑93 (MQ=255)
                                                          aCGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAAc                                  <  1:3431056/93‑1 (MQ=255)
                                                            gTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGt                                >  1:5971268/1‑93 (MQ=255)
                                                            gTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGt                                >  1:4967023/1‑93 (MQ=255)
                                                            gTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑C‑‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTc                               >  1:1074229/1‑93 (MQ=255)
                                                              tGAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCa                              <  2:2876256/93‑1 (MQ=255)
                                                               gAGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCaa                             <  2:1155396/93‑1 (MQ=255)
                                                                aGCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAAt                            >  1:3589612/1‑93 (MQ=255)
                                                                 gCGGCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATc                           >  2:7186958/1‑93 (MQ=255)
                                                                 gCGGCTGGTGGGTAGTCAGAATAACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATc                           >  1:6574413/1‑93 (MQ=255)
                                                                  cGGCTGGTGGGTAGTCAGAATCACAACC‑‑CCCCCCCTGCTCCGTATGCTGCGCCATCCGCTGGGTCAGACGATCGACACCGTTAACGTCAATcg                          <  1:6259701/93‑1 (MQ=255)
                                                                   ggCTGGTGGGTAGTCAGAATCACAA‑C‑‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATcgcg                        <  1:1776091/93‑1 (MQ=255)
                                                                    gCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATcgcg                        <  2:826216/93‑1 (MQ=255)
                                                                    gCTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATcgcg                        >  2:2099686/1‑93 (MQ=255)
                                                                    gCTGGGGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAAAGTCAATcgcg                        >  1:1485636/1‑93 (MQ=255)
                                                                     cTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCgg                       <  2:3048012/93‑1 (MQ=255)
                                                                     cTGGTGGGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCgg                       <  1:1661455/93‑1 (MQ=255)
                                                                     cTGGTGGGTAGTCAGAATCACAA‑C‑‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGt                      <  2:3342298/93‑1 (MQ=255)
                                                                      tggtggGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGt                      >  1:5225675/1‑93 (MQ=255)
                                                                       ggtggGTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCCGCGCCATACGCTGGGTCAGACGATCCACACCGTTAACGTCAATCGCGGTa                     >  2:6841726/1‑93 (MQ=255)
                                                                          gggTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTaaaa                  <  2:2202880/93‑1 (MQ=255)
                                                                          gggTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTaaaa                  >  1:4170258/1‑93 (MQ=255)
                                                                          gggTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGAGCGACACCGTTAACGTCAATCGCGGTaaaa                  <  1:6195754/93‑1 (MQ=255)
                                                                            gTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAgg                <  2:2322234/93‑1 (MQ=255)
                                                                            gTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAgg                >  1:970597/1‑93 (MQ=255)
                                                                            gTAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAgg                <  2:1404325/93‑1 (MQ=255)
                                                                             tAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGc               <  1:5910118/93‑1 (MQ=255)
                                                                             tAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGc               <  2:4007189/93‑1 (MQ=255)
                                                                             tAGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGc               <  1:3891082/93‑1 (MQ=255)
                                                                              aGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCt              <  1:7265102/93‑1 (MQ=255)
                                                                              aGTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCt              <  1:7197450/93‑1 (MQ=255)
                                                                               gTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCtc             <  1:4363966/93‑1 (MQ=255)
                                                                               gTCAGAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCtc             >  1:3578682/1‑93 (MQ=255)
                                                                                tCAGAATCACAA‑C‑‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCtcgt           >  1:7506869/1‑93 (MQ=255)
                                                                                   gAATCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCtcgtcg         >  2:5037299/1‑93 (MQ=255)
                                                                                    aaTCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCTCCATACGCTGGGTCAGACGATCTACACCGTTAACGTCAATCGCGGTAAAAGCCTCGTCGa        <  1:1024694/93‑1 (MQ=255)
                                                                                     aTCACAA‑TCCCCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGa        >  2:7233371/1‑93 (MQ=255)
                                                                                      tCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAgg      <  2:109533/93‑1 (MQ=255)
                                                                                      tCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAgg      <  1:505436/93‑1 (MQ=255)
                                                                                      tCACAA‑TC‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAgg      >  1:5446931/1‑93 (MQ=255)
                                                                                              |                                                                                          
AGAACATCATGCGGCCCTCGTTTGCGTCAGTGAAATGCGGCGAATTTTACTTTCAGCAACGTTGAGCGGCTGGTGGGTAGTCAGAATCACAA‑T‑‑CCCCCCCTGCTCCGTATGCTGCGCCATACGCTGGGTCAGACGATCGACACCGTTAACGTCAATCGCGGTAAAAGGCTCGTCGAGGATCC  >  NC_000913/2295035‑2295216

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