Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,352,673 T→C E118G (GAA→GGA)  arcB ← aerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,352,6730TC96.8% 326.2 / 0.3 94E118G (GAA→GGA) arcBaerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (44/47);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

GAGCTGAATTTGTGTCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTTCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGA  >  NC_000913/3352576‑3352768
                                                                                                 |                                                                                               
gAGCTGAATTTGTGTCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTTcc                                                                                               >  2:747651/1‑100 (MQ=255)
 aGCTGAATTTGTGTCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCt                                                                                              >  2:2032860/1‑100 (MQ=255)
  gCTGAATTTGTGTCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTg                                                                                             <  2:225022/100‑1 (MQ=255)
  gCTGAATTTGTGTCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTg                                                                                             <  2:29779/100‑1 (MQ=255)
    tGAATTTGTGTCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGAt                                                                                           <  1:2001006/100‑1 (MQ=255)
    tGAATTTGTGTCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCACTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGAt                                                                                           <  1:1744344/100‑1 (MQ=255)
     gAATTTGTGTCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGAtt                                                                                          >  2:657604/1‑100 (MQ=255)
      aaTTTGTGTCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGAttn                                                                                         <  1:1197113/100‑2 (MQ=255)
       aTTTGTGTCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTa                                                                                        <  1:2719094/100‑1 (MQ=255)
         ttGTGTCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGc                                                                                      <  2:1291001/100‑1 (MQ=255)
         ttGTGTCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGc                                                                                      >  2:1018526/1‑100 (MQ=255)
         ttGTGTCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGc                                                                                      <  2:875335/100‑1 (MQ=255)
          tgtgTCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCt                                                                                     <  1:731535/100‑1 (MQ=255)
           gtgtCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTg                                                                                    <  2:2232988/100‑1 (MQ=255)
           gtgtCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTg                                                                                    >  2:1320985/1‑100 (MQ=255)
             gtCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGATCGGCAATTCCCTGATTTAGCTggg                                                                                  >  1:2160413/1‑100 (MQ=255)
             gtCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTggg                                                                                  >  2:1628597/1‑100 (MQ=255)
              tctcTTCGCGCTCTTTGATTTCAATTTTCTGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTggg                                                                                  <  1:1149679/99‑1 (MQ=255)
              tctcTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGc                                                                                 <  2:1798851/100‑1 (MQ=255)
              tctcTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGc                                                                                 <  1:826906/100‑1 (MQ=255)
              tctcTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGc                                                                                 <  2:2599755/100‑1 (MQ=255)
               ctctTCGCGCTCTTTGATTTCAATTTTCCGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGc                                                                                 <  2:1187103/99‑1 (MQ=255)
               ctctTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCa                                                                                <  2:1855281/100‑1 (MQ=255)
                 ctTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAat                                                                              >  1:1617616/1‑100 (MQ=255)
                  ttCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAata                                                                             >  1:1379500/1‑100 (MQ=255)
                   tCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAatat                                                                            <  2:519278/100‑1 (MQ=255)
                   tCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAatat                                                                            <  2:34252/100‑1 (MQ=255)
                      cgcTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAAtattat                                                                         >  2:1667096/1‑100 (MQ=255)
                        ctctTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCt                                                                       >  1:2781256/1‑100 (MQ=255)
                        ctctTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCt                                                                       >  1:796070/1‑100 (MQ=255)
                           tttGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCtttt                                                                    >  1:19539/1‑100 (MQ=255)
                           tttGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCtttt                                                                    <  1:2453877/100‑1 (MQ=255)
                             tGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTaa                                                                  >  2:300314/1‑100 (MQ=255)
                              gATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCACTTCCCTGATTTAGCTGGGCAATATTATCTTTTAAc                                                                 >  1:1159897/1‑100 (MQ=255)
                              gATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAAc                                                                 <  2:1598557/100‑1 (MQ=255)
                              gATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAAc                                                                 <  1:1926172/100‑1 (MQ=255)
                               aTTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACt                                                                <  1:300310/100‑1 (MQ=255)
                                tttCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTg                                                               >  2:914802/1‑100 (MQ=255)
                                tttCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTg                                                               <  1:2318549/100‑1 (MQ=255)
                                 ttCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGa                                                              >  2:1239707/1‑100 (MQ=255)
                                     aTTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGt                                                          >  2:1812015/1‑100 (MQ=255)
                                        ttCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGa                                                       <  2:284646/100‑1 (MQ=255)
                                           aGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGc                                                    >  2:931508/1‑100 (MQ=255)
                                           aGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGc                                                    >  2:2528851/1‑100 (MQ=255)
                                            gTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCt                                                   >  2:2609535/1‑100 (MQ=255)
                                            gTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCt                                                   >  1:1311564/1‑100 (MQ=255)
                                            gTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCt                                                   >  2:452792/1‑100 (MQ=255)
                                              tGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCa                                                 <  2:865127/100‑1 (MQ=255)
                                               ggCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCATGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCaa                                                >  1:794691/1‑100 (MQ=255)
                                               ggCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCaa                                                >  1:91968/1‑100 (MQ=255)
                                                gCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCaaa                                               >  2:2533399/1‑100 (MQ=255)
                                                 ccGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAAt                                              >  2:2028226/1‑100 (MQ=255)
                                                  cGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAAtc                                             <  2:2789452/100‑1 (MQ=255)
                                                  cGAAGGTTTCCTGCAGTTCTGCTTCCGCCTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAAtc                                             <  1:1787120/100‑1 (MQ=255)
                                                   gAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAAtcg                                            <  2:2182462/100‑1 (MQ=255)
                                                   gAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAAtcg                                            >  2:1232791/1‑100 (MQ=255)
                                                     aGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAAtcgcg                                          <  2:2623440/100‑1 (MQ=255)
                                                     aGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAAtcgcg                                          <  1:980235/100‑1 (MQ=255)
                                                      ggTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAAtcgcgc                                         >  2:114124/1‑100 (MQ=255)
                                                       gTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAAtcgcgct                                        >  1:2621374/1‑100 (MQ=255)
                                                        tttCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAAtcgcgctc                                       <  1:1065462/100‑1 (MQ=255)
                                                           ccTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAAtcgcgctcgcg                                    <  2:1617628/100‑1 (MQ=255)
                                                            cTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAAtcgcgctcgcgc                                   <  2:984458/100‑1 (MQ=255)
                                                             tGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCa                                  <  2:2393604/100‑1 (MQ=255)
                                                             tGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCa                                  <  2:1972666/100‑1 (MQ=255)
                                                              gCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCAt                                 >  2:828101/1‑100 (MQ=255)
                                                               cAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCAtc                                <  2:2335193/100‑1 (MQ=255)
                                                                     tGCTTCCGCTTTTTCACGAACCGCAATTCCCTGATTTAGCTGTGCAATATTATCTTTTAACTGCACGCTGAGTCTCAAATCGCGCTCGCGCATctcctc                           <  1:2430178/99‑1 (MQ=255)
                                                                         tCCGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGtt                      <  2:1550781/100‑1 (MQ=255)
                                                                          ccGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGttt                     <  1:2362253/100‑1 (MQ=255)
                                                                           cGCTTTTTCACGAACGGCAATTTCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGtttt                    <  1:747645/100‑1 (MQ=255)
                                                                           cGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGtttt                    >  2:992060/1‑100 (MQ=255)
                                                                           cGCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGtttt                    <  1:1667084/100‑1 (MQ=255)
                                                                            gCTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGttttt                   <  2:2271110/100‑1 (MQ=255)
                                                                             cTTTTTCTCGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGtttttn                  <  2:760390/100‑2 (MQ=255)
                                                                             cTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTg                  <  2:2160428/100‑1 (MQ=255)
                                                                             cTTTTTCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTg                  <  1:250002/100‑1 (MQ=255)
                                                                               ttttCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCa                >  1:521209/1‑100 (MQ=255)
                                                                               ttttCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCa                <  1:1320977/100‑1 (MQ=255)
                                                                               ttttCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCa                >  1:2128576/1‑100 (MQ=255)
                                                                                 ttCACGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCAcc              >  1:1010851/1‑100 (MQ=255)
                                                                                     cGAACGGCAATTCCCTGATTTAGCTGGGCCATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGGGCATCTCCTCCAGTTTTTGCACCAGcc          >  2:371241/1‑100 (MQ=255)
                                                                                     cGAACGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGcc          >  2:149018/1‑100 (MQ=255)
                                                                                       aaCGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGt        >  1:1734456/1‑100 (MQ=255)
                                                                                        aCGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTg       >  2:435431/1‑100 (MQ=255)
                                                                                        aCGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTg       >  1:146289/1‑100 (MQ=255)
                                                                                        aCGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTg       >  1:949789/1‑100 (MQ=255)
                                                                                        aCGGCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTg       >  2:435960/1‑100 (MQ=255)
                                                                                          ggCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGAc     >  2:2387463/1‑100 (MQ=255)
                                                                                           gCAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACa    >  2:214189/1‑100 (MQ=255)
                                                                                            cAATTTCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAg   >  2:587001/1‑100 (MQ=255)
                                                                                            cAATTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAg   <  2:91969/100‑1 (MQ=255)
                                                                                             aaTTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAg   >  1:599290/1‑99 (MQ=255)
                                                                                             aaTTCCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGa  <  2:823316/100‑1 (MQ=255)
                                                                                                 |                                                                                               
GAGCTGAATTTGTGTCTCTTCGCGCTCTTTGATTTCAATTTTCAGTTGGCCGAAGGTTTCCTGCAGTTCTGCTTCCGCTTTTTCACGAACGGCAATTTCCTGATTTAGCTGGGCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGA  >  NC_000913/3352576‑3352768

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