Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 298,837 C→T A700T (GCG→ACG)  paoC ← aldehyde dehydrogenase: molybdenum cofactor‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913298,8370CT99.7% 1290.5 / ‑4.4 306A700T (GCG→ACG) paoCaldehyde dehydrogenase: molybdenum cofactor‑binding subunit
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (138/167);  total (138/168)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.05e-01

Only 100 of 308 total aligned reads displayed.

AACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG  >  NC_000913/298739‑298934
                                                                                                  |                                                                                                 
aaCCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCgtgc                                                                                                 >  2:509574/1‑101 (MQ=255)
 aCCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCt                                                                                                <  1:509021/101‑1 (MQ=255)
  ccACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTc                                                                                               <  2:509073/101‑1 (MQ=255)
     caTCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCAcg                                                                                            <  2:509267/101‑1 (MQ=255)
      aTCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCa                                                                                              >  2:509583/1‑98 (MQ=255)
      aTCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCAcgc                                                                                           >  1:509582/1‑101 (MQ=255)
       tCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCAcgcc                                                                                          <  1:508914/101‑1 (MQ=255)
       tCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCAcgcc                                                                                          <  1:509001/101‑1 (MQ=255)
          ggCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCAcgcc                                                                                          >  1:509589/1‑98 (MQ=255)
          ggCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGca                                                                                       >  2:509588/1‑101 (MQ=255)
             aGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAg                                                                                    >  1:509593/1‑101 (MQ=255)
                ttATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGcc                                                                                 >  1:509596/1‑101 (MQ=255)
                ttATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGcc                                                                                 <  1:509214/101‑1 (MQ=255)
                  aTCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCa                                                                               <  2:509198/101‑1 (MQ=255)
                  aTCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCa                                                                               >  1:509600/1‑101 (MQ=255)
                  aTCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCa                                                                               <  1:509055/101‑1 (MQ=255)
                    cGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGc                                                                             <  2:509146/101‑1 (MQ=255)
                    cGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGc                                                                             >  2:509601/1‑101 (MQ=255)
                      agcagcTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTc                                                                           >  1:509605/1‑101 (MQ=255)
                      agcagcTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTc                                                                           <  2:508906/101‑1 (MQ=255)
                         agcTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCAcc                                                                        >  2:509610/1‑101 (MQ=255)
                           cTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGa                                                                      >  1:509612/1‑101 (MQ=255)
                              aTCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGacac                                                                   >  2:509615/1‑101 (MQ=255)
                                ccAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCt                                                                 >  2:509621/1‑101 (MQ=255)
                                ccAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCt                                                                 >  1:509620/1‑101 (MQ=255)
                                  agagTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCttt                                                               <  1:509118/101‑1 (MQ=255)
                                  agagTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCttt                                                               <  2:509097/101‑1 (MQ=255)
                                    agTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTgg                                                             <  2:509229/101‑1 (MQ=255)
                                    agTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTgg                                                             >  1:509625/1‑101 (MQ=255)
                                     gTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGc                                                            <  1:509129/101‑1 (MQ=255)
                                       gATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCt                                                          >  2:509629/1‑101 (MQ=255)
                                         tGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTc                                                        <  1:509034/101‑1 (MQ=255)
                                          gggATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCa                                                       >  1:509631/1‑101 (MQ=255)
                                           ggATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCAt                                                      <  2:509248/101‑1 (MQ=255)
                                           ggATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCAt                                                      <  1:509238/101‑1 (MQ=255)
                                            gATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATc                                                     >  1:509634/1‑101 (MQ=255)
                                              tAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCgg                                                   <  1:509079/101‑1 (MQ=255)
                                                aTCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCgggg                                                 >  1:509638/1‑101 (MQ=255)
                                                  cccGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGggag                                               >  2:509642/1‑101 (MQ=255)
                                                  cccGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGggag                                               >  1:509643/1‑101 (MQ=255)
                                                   ccGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGggagg                                              <  1:509189/101‑1 (MQ=255)
                                                   ccGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGggagg                                              <  2:509105/101‑1 (MQ=255)
                                                    cGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGggagga                                             <  1:509179/101‑1 (MQ=255)
                                                     gTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGat                                            >  1:509648/1‑101 (MQ=255)
                                                      tACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGata                                           <  2:509319/101‑1 (MQ=255)
                                                      tACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGata                                           >  1:509651/1‑101 (MQ=255)
                                                      tACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGata                                           >  2:509652/1‑101 (MQ=255)
                                                       aCCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGat                                            >  2:509656/1‑99 (MQ=255)
                                                        cccGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGgga                                                <  2:509114/94‑1 (MQ=255)
                                                         ccGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATgg                                        >  1:509660/1‑101 (MQ=255)
                                                         ccGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATgg                                        <  1:508826/101‑1 (MQ=255)
                                                          cGAATACCGGTGGCGTTATACACCGCGTTGGCGATAACCGTGCTCACGCCGCACAGGCCCAGCTCa                                                                          >  1:509663/1‑66 (MQ=255)
                                                           gAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATgggg                                      >  2:509664/1‑101 (MQ=255)
                                                            aaTACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGt                                     <  2:509127/101‑1 (MQ=255)
                                                            aaTACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGt                                     <  1:509308/101‑1 (MQ=255)
                                                            aaTACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGt                                     <  1:509520/101‑1 (MQ=255)
                                                             aTACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTc                                    <  1:509259/101‑1 (MQ=255)
                                                              tACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTg                                                              <  1:509087/74‑1 (MQ=255)
                                                              tACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCg                                   >  1:509668/1‑101 (MQ=255)
                                                               aCCGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCgggg                                                 >  2:509670/1‑86 (MQ=255)
                                                                ccGGTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGTGTCGGt                                 <  1:509157/101‑1 (MQ=255)
                                                                   gTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCgg                                                   <  2:509049/80‑1 (MQ=255)
                                                                   gTGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTatc                              <  1:509269/101‑1 (MQ=255)
                                                                    tGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTatca                             <  1:509351/101‑1 (MQ=255)
                                                                    tGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTatca                             <  1:509228/101‑1 (MQ=255)
                                                                    tGGCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTatca                             >  2:509676/1‑101 (MQ=255)
                                                                     ggCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTatcat                            <  2:509092/101‑1 (MQ=255)
                                                                     ggCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTatcat                            <  1:509140/101‑1 (MQ=255)
                                                                      gCGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTatcatc                           <  2:509271/101‑1 (MQ=255)
                                                                       cGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATcc                          >  2:509679/1‑101 (MQ=255)
                                                                       cGTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATcc                          >  1:509680/1‑101 (MQ=255)
                                                                        gTTATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCa                         <  2:509011/101‑1 (MQ=255)
                                                                         ttATACACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCgg                                  >  1:509685/1‑91 (MQ=255)
                                                                          tataCACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAgg                       <  1:509075/101‑1 (MQ=255)
                                                                          tataCACCGCGTTGGAGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGgga                                                >  2:509687/1‑76 (MQ=255)
                                                                            taCACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGaa                     <  2:509183/101‑1 (MQ=255)
                                                                            taCACCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGaa                     <  1:509283/101‑1 (MQ=255)
                                                                             acacCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGaaa                    >  1:509690/1‑101 (MQ=255)
                                                                              cacCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGaaaa                   <  1:509442/101‑1 (MQ=255)
                                                                              cacCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGaaaa                   <  1:509294/101‑1 (MQ=255)
                                                                               acCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAt                  >  1:509693/1‑101 (MQ=255)
                                                                               acCGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAAt                  <  1:509357/101‑1 (MQ=255)
                                                                                ccGCGTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATc                 <  2:508994/101‑1 (MQ=255)
                                                                                 cgcgTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCa                <  1:509262/101‑1 (MQ=255)
                                                                                   cgTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCAACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCAcc              >  2:509695/1‑101 (MQ=255)
                                                                                    gTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCAcct             >  2:509698/1‑101 (MQ=255)
                                                                                    gTTGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCAcct             <  2:509207/101‑1 (MQ=255)
                                                                                      tGGCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCAcctcc           >  2:509702/1‑101 (MQ=255)
                                                                                       ggCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCAcctcct          <  1:509197/101‑1 (MQ=255)
                                                                                       ggCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCAc               <  1:509332/96‑1 (MQ=255)
                                                                                        gCGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTg         <  2:509160/101‑1 (MQ=255)
                                                                                         cGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGt        <  1:509577/101‑1 (MQ=255)
                                                                                         cGATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGt        <  1:509323/101‑1 (MQ=255)
                                                                                          gATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGtt       >  1:509707/1‑101 (MQ=255)
                                                                                          gATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGtt       <  2:509360/101‑1 (MQ=255)
                                                                                          gATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGtt       <  2:509290/101‑1 (MQ=255)
                                                                                          gATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGtt       <  2:509136/101‑1 (MQ=255)
                                                                                          gATAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGtt       >  2:509706/1‑101 (MQ=255)
                                                                                           aTAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGttt      >  1:509711/1‑101 (MQ=255)
                                                                                            tAGCCGTGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGtttt     <  2:509239/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
AACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG  >  NC_000913/298739‑298934

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑