Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,764,687 T→C 16.3% intergenic (‑301/+247) sufA ← / ← ydiH Fe‑S cluster assembly protein/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,764,6870TC16.3% 294.1 / 48.0 135intergenic (‑301/+247)sufA/ydiHFe‑S cluster assembly protein/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (51/62);  new base C (10/12);  total (61/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

Only 100 of 137 total aligned reads displayed.

GCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCAT  >  NC_000913/1764588‑1764785
                                                                                                   |                                                                                                   
 cATTAAGGCTGTGATGAAGAAGAACAAGTTGTATTTAAAATGACCATTCAAGAGACAGGAAAGAAATTGCAGCTTTCAGAAATAAGAAAAACCTAAAGTc                                                                                                    >  2:1062175/1‑100 (MQ=255)
  aTTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGGAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTc                                                                                                    >  2:2913239/1‑99 (MQ=255)
      aGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑c                                                                                              <  2:968495/100‑1 (MQ=255)
       ggCTGTGATGAAAAAGAAAAAGTTTTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGAGTTCAGATATAAGAAAACTCTTAAGTCTGAG‑Cg                                                                                             >  2:765577/1‑100 (MQ=255)
       ggCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑Cg                                                                                             >  2:3689593/1‑100 (MQ=255)
       ggCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑Cg                                                                                             <  1:3689621/100‑1 (MQ=255)
         cTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGac                                                                                           <  1:682571/100‑1 (MQ=255)
         cTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGac                                                                                           <  1:1905480/100‑1 (MQ=255)
            tgatgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGAcaca                                                                                        >  2:1727991/1‑100 (MQ=255)
             gatgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAg                                                                                       >  1:3266259/1‑100 (MQ=255)
             gatgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAg                                                                                       >  1:641237/1‑100 (MQ=255)
             gatgaAAAAGAAAAACTTGTATTGAGAATGACTATTTAAGAGTTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAg                                                                                       <  1:3756791/100‑1 (MQ=255)
                 aaaaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGc                                                                                     <  1:3572362/98‑1 (MQ=255)
                  aaaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTa                                                                                  >  1:2266347/1‑100 (MQ=255)
                   aaagaaaaagTTGTATTGATATTGACTAGTTAAGGGATAGGTAAAAAAGTGCACCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCt                                                                                    <  2:2575264/97‑1 (MQ=255)
                    aagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAg                                                                                >  1:160916/1‑100 (MQ=255)
                    aagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTG‑TCGCACAGACTTAAg                                                                                <  1:3002920/100‑1 (MQ=37)
                     agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAgg                                                                               >  1:3245556/1‑100 (MQ=255)
                     agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAgg                                                                               >  1:2115028/1‑100 (MQ=255)
                     agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAgg                                                                               >  1:1501054/1‑100 (MQ=255)
                     agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAgg                                                                               >  1:1501457/1‑100 (MQ=255)
                     agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAgg                                                                               <  1:1695803/100‑1 (MQ=255)
                     agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAgg                                                                               <  2:3075426/100‑1 (MQ=37)
                      gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAggg                                                                              >  1:3564510/1‑100 (MQ=255)
                      gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAggg                                                                              >  2:1104571/1‑100 (MQ=255)
                      gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAggg                                                                              <  1:2716003/100‑1 (MQ=37)
                      gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAggg                                                                              <  1:210569/100‑1 (MQ=37)
                       aaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGt                                                                             <  2:1956248/100‑1 (MQ=255)
                       aaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAGGGt                                                                             <  2:130453/100‑1 (MQ=37)
                       aaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTG‑TCGCACAGACTTAAgg                                                                               <  2:1320404/98‑1 (MQ=37)
                          aagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTc                                                                          <  1:2307808/100‑1 (MQ=255)
                          aagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTc                                                                          <  2:117774/100‑1 (MQ=255)
                          aagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAGGGtttt                                                                          <  2:3596908/100‑2 (MQ=25)
                            gTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTa                                                                        >  1:3115118/1‑100 (MQ=255)
                            gTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAGGGttttct                                                                        <  2:1537375/100‑4 (MQ=25)
                              tGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTAcc                                                                      <  1:1410462/100‑1 (MQ=255)
                              tGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTAcc                                                                      >  2:3731330/1‑100 (MQ=255)
                                 aTTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCa                                                                   >  2:997927/1‑100 (MQ=255)
                                 aTTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCa                                                                   >  1:2074398/1‑100 (MQ=255)
                                   tGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGttttcttatttct                                                                 >  1:3215266/1‑90 (MQ=255)
                                   tGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATc                                                                 <  1:664500/100‑1 (MQ=255)
                                    gAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATcc                                                                >  1:2849973/1‑100 (MQ=255)
                                    gAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATcc                                                                >  1:505084/1‑100 (MQ=255)
                                     aaaaTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCg                                                               <  2:2702753/100‑1 (MQ=255)
                                      aaaTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCgg                                                              <  1:1899724/100‑1 (MQ=255)
                                        aTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGgcg                                                            <  1:3024710/100‑1 (MQ=255)
                                         tGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGgcgc                                                           <  1:951574/100‑1 (MQ=255)
                                          gACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACATGCTTAAGGGTTTCTACCCCATCCGgcg                                                            <  2:2772026/98‑1 (MQ=255)
                                          gACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCt                                                          <  1:3557052/100‑1 (MQ=255)
                                           aCTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCtt                                                         >  1:190459/1‑100 (MQ=255)
                                            cTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTa                                                        <  2:1806856/100‑1 (MQ=255)
                                            cTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTa                                                        <  2:293916/100‑1 (MQ=255)
                                              aTTTAAGAGATAGTTAAAAAAGGGCAGCTGTCAGAAATAAGATAACCCTTAATTCTGTT‑CTACACATTCTTAAGGGTTTCTACCCCATCCGGCGCTTAtc                                                      <  1:2224504/100‑1 (MQ=255)
                                              aTTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTAtc                                                      >  1:1259844/1‑100 (MQ=255)
                                                 tAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTT‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTcc                                                   >  1:2914623/1‑100 (MQ=255)
                                                 tAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTcc                                                   >  1:3447696/1‑100 (MQ=255)
                                                 tAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTcc                                                   <  2:324091/100‑1 (MQ=255)
                                                   agagATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCgg                                                 <  2:2392400/100‑1 (MQ=255)
                                                    gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGc                                                <  1:1080802/100‑1 (MQ=255)
                                                    gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGc                                                >  1:3301544/1‑100 (MQ=255)
                                                    gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTGAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGc                                                <  2:1978801/100‑1 (MQ=255)
                                                            tAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCACATCCGGCGCTTATCTCCGGCACTCTCAg                                        >  2:866063/1‑100 (MQ=255)
                                                             aaaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGt                                       >  1:1679749/1‑100 (MQ=255)
                                                              aaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACTCCATCCGGCGCTTATCTCCGACACTCTCAGTg                                      >  2:2964052/1‑100 (MQ=255)
                                                              aaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTg                                      <  1:992525/100‑1 (MQ=255)
                                                              aaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTg                                      >  2:435098/1‑100 (MQ=255)
                                                              aaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTg                                      <  1:1314196/100‑1 (MQ=255)
                                                              aaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTg                                      >  1:1755066/1‑100 (MQ=255)
                                                                aaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGTCACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGc                                    <  2:2442860/100‑1 (MQ=255)
                                                                     gCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGc                               >  2:1786887/1‑100 (MQ=255)
                                                                      naGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGct                              >  2:3535402/2‑100 (MQ=255)
                                                                      cAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGct                              <  1:2135544/100‑1 (MQ=255)
                                                                      cAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGc                               >  2:2851090/1‑99 (MQ=255)
                                                                         cGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCtt                           <  1:3484619/100‑1 (MQ=255)
                                                                          gTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTg                          <  2:2737492/100‑1 (MQ=255)
                                                                           ttCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGa                         <  2:340342/100‑1 (MQ=255)
                                                                              aGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAgg                      <  1:2438006/100‑1 (MQ=255)
                                                                                 aaTAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAgggga                   <  1:2622074/100‑2 (MQ=255)
                                                                                 aaTAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGc                   >  1:890023/1‑100 (MQ=255)
                                                                                   tAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCg                  >  2:3507271/1‑99 (MQ=255)
                                                                                    gtagaaaCCCTTAAGCCTGTG‑TCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTCGCTTAGCTCTTGAAGGGGCGa                 >  2:2651532/5‑99 (MQ=21)
                                                                                    aaGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGAt                <  2:3793661/100‑1 (MQ=255)
                                                                                     aGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGAt                <  2:2914602/99‑1 (MQ=255)
                                                                                      gAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATaa              >  1:1304731/1‑100 (MQ=255)
                                                                                       aaaaCCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAg             <  2:1501190/100‑1 (MQ=255)
                                                                                         aaCCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGTCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGa            <  2:1969594/99‑1 (MQ=255)
                                                                                         aaCCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaa           <  2:3115094/100‑1 (MQ=255)
                                                                                          aCCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaat          <  1:3499576/100‑1 (MQ=255)
                                                                                          aCCCTTAAGCCTGTGTCG‑CACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaat          >  1:1668922/1‑100 (MQ=37)
                                                                                           cccTTAAGCCTGTGTCG‑CACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaata         >  1:884257/1‑100 (MQ=37)
                                                                                           cccTTAAGCCTGTGTCG‑CACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaata         >  1:2522086/1‑100 (MQ=37)
                                                                                           cccTTAAGCCTGTGTCG‑CACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaata         >  1:1093315/1‑100 (MQ=37)
                                                                                           cccTTAAGCCTGTGTCG‑CACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaata         >  1:1066257/1‑100 (MQ=37)
                                                                                           cccTTAAGCCTGTG‑TCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaata         >  2:873997/1‑100 (MQ=37)
                                                                                            ccTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataa        >  1:678403/1‑100 (MQ=255)
                                                                                            ccTTAAGCCTGTGTCG‑CACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataa        >  1:606890/1‑100 (MQ=37)
                                                                                            ccTTAAGCCTGTGTCG‑CACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataa        >  1:3320387/1‑100 (MQ=37)
                                                                                             cTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataat       >  2:2545043/1‑100 (MQ=255)
                                                                                              ttAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtc      <  2:134383/100‑1 (MQ=255)
                                                                                                aaGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtctc    <  1:998507/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCAT  >  NC_000913/1764588‑1764785

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