Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,537,405 +GC intergenic (+61/‑394) b3407 → / → feoA hypothetical protein/ferrous iron transport protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,537,4041.C99.7% 975.4 / ‑5.3 316intergenic (+60/‑395)b3407/feoAhypothetical protein/ferrous iron transport protein A
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base C (159/156);  minor base T (1/0);  total (160/156)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.30e-01
*NC_0009133,537,4042.G100.0% 991.2 / NA 318intergenic (+60/‑395)b3407/feoAhypothetical protein/ferrous iron transport protein A
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (160/158);  total (160/158)

Only 100 of 344 total aligned reads displayed.

GCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑CCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATA  >  NC_000913/3537312‑3537487
                                                                                             ||                                                                                   
gcgcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCa                                                                                       <  1:17203/93‑1 (MQ=255)
 cgcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑c                                                                                    >  1:6555559/1‑93 (MQ=255)
 cgcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑c                                                                                    <  1:619757/93‑1 (MQ=255)
  gcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑cg                                                                                   >  1:5864902/1‑92 (MQ=255)
   cTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCGACGGTGCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑cgc                                                                                  <  2:853829/93‑3 (MQ=255)
    tGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑cgcc                                                                                 <  2:1410499/93‑4 (MQ=255)
    tGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑cgcc                                                                                 >  2:1969184/1‑90 (MQ=255)
            gcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTgg                                                                           <  1:4394881/93‑1 (MQ=38)
            gcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTgg                                                                           >  2:4846586/1‑93 (MQ=38)
              aaTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCa                                                                         >  1:986726/1‑93 (MQ=37)
              aaTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCa                                                                         >  1:3621816/1‑93 (MQ=37)
              aaTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCa                                                                         >  2:2709113/1‑93 (MQ=37)
               aTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaa                                                                        <  1:1056887/93‑1 (MQ=38)
                 gggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacg                                                                      <  2:3092415/93‑1 (MQ=37)
                  ggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacga                                                                     >  1:5344724/1‑93 (MQ=37)
                  ggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacga                                                                     <  1:1989284/93‑1 (MQ=37)
                   gCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacgaa                                                                    <  1:1646019/93‑1 (MQ=37)
                   gCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacgaa                                                                    <  1:2213543/93‑1 (MQ=37)
                   gCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacgaa                                                                    <  1:3548396/93‑1 (MQ=38)
                     aaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacgaacg                                                                  >  1:6980270/1‑93 (MQ=37)
                     aaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCaacgaacg                                                                  >  1:6042427/1‑93 (MQ=37)
                      aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGc                                                                 >  2:6446618/1‑93 (MQ=37)
                      aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGc                                                                 >  1:4325511/1‑93 (MQ=37)
                      aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGc                                                                 >  2:7910033/1‑93 (MQ=37)
                      aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGc                                                                 <  1:6625772/93‑1 (MQ=37)
                       aaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGcc                                                                >  1:2354512/1‑93 (MQ=37)
                       aaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGcc                                                                <  2:4592586/93‑1 (MQ=37)
                        aaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCt                                                               <  2:1701476/93‑1 (MQ=37)
                        aaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCt                                                               <  1:7437486/93‑1 (MQ=37)
                        aaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCt                                                               <  1:4884965/93‑1 (MQ=37)
                         aaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCtt                                                              >  2:5417553/1‑93 (MQ=37)
                          aCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTg                                                             >  1:4653355/1‑93 (MQ=37)
                               aaaCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGa                                                        >  2:5360466/1‑93 (MQ=37)
                                aaCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGAc                                                       >  1:7562454/1‑93 (MQ=38)
                                     ccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaaa                                                  >  2:3003830/1‑93 (MQ=37)
                                      cGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaaaa                                                 <  1:5511305/93‑1 (MQ=37)
                                      cGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaaaa                                                 <  2:640696/93‑1 (MQ=37)
                                       gtgcTGGCCTAGGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAg                                                <  2:5852429/90‑1 (MQ=25)
                                       gTACTGGCCTACGGTTCGAATTTGCACGAAATCGTTGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAg                                                >  1:3308941/1‑93 (MQ=38)
                                       gTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAg                                                <  1:3489154/93‑1 (MQ=37)
                                        tACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGt                                               >  2:1291100/1‑93 (MQ=37)
                                        tACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGt                                               >  1:7117809/1‑93 (MQ=37)
                                        tACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGt                                               >  1:1219533/1‑93 (MQ=38)
                                         aCTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTg                                              >  1:5682833/1‑93 (MQ=37)
                                          cTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCATGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGc                                             >  1:1596339/1‑93 (MQ=38)
                                          cTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGc                                             <  2:3763458/93‑1 (MQ=37)
                                           tGGCGTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGcc                                            <  2:7606551/93‑1 (MQ=37)
                                           tGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGcc                                            <  1:7207358/93‑1 (MQ=37)
                                           tGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGcc                                            >  2:2599316/1‑93 (MQ=37)
                                              ccTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGga                                         >  1:195972/1‑93 (MQ=38)
                                              ccTACGATTCTATTTAGATCGAAATCGTAGGTCAGATAAGGGGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCTGga                                         <  2:2288104/93‑1 (MQ=16)
                                               gtACGGTTCTAATTTGCACGAAATCGTAGGTCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaa                                        <  1:3900550/92‑2 (MQ=25)
                                               cTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgag                                        >  2:6598227/1‑93 (MQ=38)
                                                tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                                       <  1:6401575/93‑1 (MQ=37)
                                                tACGATTCAAATTTCCAGAAAACCGTGGGCAAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                                       <  1:2851403/93‑1 (MQ=16)
                                                 aCGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGaa                                      <  2:2755562/93‑1 (MQ=37)
                                                   ggTTCGATTTTGCACAAAATCGTAGGGCAAATAAGGCTTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAata                                    <  1:1452375/93‑1 (MQ=37)
                                                   ggTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAata                                    <  2:1203689/93‑1 (MQ=37)
                                                   ggTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAata                                    <  2:572144/93‑1 (MQ=37)
                                                    gTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAatat                                   <  1:2649201/93‑1 (MQ=37)
                                                    gTTCGAATTTGAGCGAAATCGTAGGTCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAatat                                   <  2:3544635/93‑1 (MQ=25)
                                                      tCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtct                                 <  2:5864903/93‑1 (MQ=37)
                                                       cGAATTTGCACGAAATCGTAGGACAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtctc                                <  2:666326/93‑1 (MQ=37)
                                                        gAATTTGCACTAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                               <  1:6809232/93‑1 (MQ=37)
                                                        gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                               >  1:5001150/1‑93 (MQ=38)
                                                        gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                               >  2:736286/1‑93 (MQ=38)
                                                        gAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                               <  1:3195411/93‑1 (MQ=37)
                                                        gAATTTGCACGAAATCGTAGGGAAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                               <  2:7248129/93‑1 (MQ=37)
                                                         aaTTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCg                              >  1:5426427/1‑93 (MQ=38)
                                                         aaTTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCg                              <  1:1905796/93‑1 (MQ=37)
                                                            ttGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCa                           <  1:865300/93‑1 (MQ=38)
                                                            ttGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCa                           <  2:2127149/93‑1 (MQ=38)
                                                              ccACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCTTCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtt                         <  2:3143337/92‑1 (MQ=37)
                                                                acaaaaaCGAAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGTCATAAAAGTGCCGGAGAATATCTCCGGCAtttt                       <  2:1068843/90‑1 (MQ=25)
                                                                  gAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTa                     <  1:7273119/93‑1 (MQ=38)
                                                                  gAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTa                     <  1:2583088/93‑1 (MQ=38)
                                                                  gAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTa                     <  1:3453013/93‑1 (MQ=38)
                                                                   aaaTCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAt                    <  2:4726284/93‑1 (MQ=38)
                                                                    aaTCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAtt                   <  1:5971393/93‑1 (MQ=38)
                                                                      tCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGGGAAGATTTTCGGCGtttttttttc                 >  2:4215862/1‑91 (MQ=21)
                                                                      tCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTcc                 >  2:7055932/1‑93 (MQ=38)
                                                                       cGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCca                <  1:79084/93‑1 (MQ=38)
                                                                        gTAGGTCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcac               <  2:515913/93‑1 (MQ=37)
                                                                        gTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcac               <  1:5125758/93‑1 (MQ=38)
                                                                         tAGGGTAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcaca              <  1:774988/93‑1 (MQ=38)
                                                                         tAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcaca              >  1:2756653/1‑93 (MQ=38)
                                                                         tAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcaca              <  1:3417226/93‑1 (MQ=38)
                                                                           gggCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGc            >  1:2995583/1‑93 (MQ=38)
                                                                             gCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCa          <  2:6047286/93‑1 (MQ=38)
                                                                             gCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCa          <  1:5290203/93‑1 (MQ=38)
                                                                              cAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCaa         >  2:2174355/1‑93 (MQ=38)
                                                                               aGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCaaa        >  2:3379789/1‑93 (MQ=38)
                                                                               aGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCaaa        <  2:5515696/93‑1 (MQ=38)
                                                                               aGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGAAATAAAAGTACCGGCGAAAAACTCCGGCATTTTTATTTCACAACCaaa        >  2:6915384/1‑93 (MQ=21)
                                                                                gATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAAc       >  2:596160/1‑93 (MQ=39)
                                                                                  tAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTc     >  1:4086762/1‑93 (MQ=39)
                                                                                  tAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTc     >  1:632314/1‑93 (MQ=39)
                                                                                  tAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTc     >  1:647447/1‑93 (MQ=39)
                                                                                    aGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCat   >  2:4410868/1‑93 (MQ=255)
                                                                                     ggCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCata  <  1:1795599/93‑1 (MQ=39)
                                                                                             ||                                                                                   
GCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCA‑‑CCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATA  >  NC_000913/3537312‑3537487

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: