Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,606,898 T→C intergenic (+6/‑85) yjjG → / → prfC putative phosphatase/peptide chain release factor RF‑3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,8980TC90.3% 453.6 / 40.4 144intergenic (+6/‑85)yjjG/prfCputative phosphatase/peptide chain release factor RF‑3
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (56/74);  minor base . (7/7);  total (63/81)
Fisher's exact test for biased strand distribution p-value = 7.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

Only 100 of 175 total aligned reads displayed.

CCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGA  >  NC_000913/4606808‑4606986
                                                                                            |                                                                                        
ccATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCcccc                                                                                        >  2:526252/1‑93 (MQ=255)
  aTCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCcccccc                                                                                      >  1:1925773/1‑93 (MQ=255)
  aTCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCcccccc                                                                                      <  1:5113106/93‑1 (MQ=255)
   tCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGcccccccc                                                                                     <  2:4691228/93‑1 (MQ=255)
    cgcgAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCcccccccc                                                                                    <  2:1283251/93‑1 (MQ=255)
    cgcgAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCcccccccc                                                                                    >  2:4051953/1‑93 (MQ=255)
     gcgAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGcccccccccc                                                                                   <  1:1331778/93‑1 (MQ=255)
       gAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGtt                                                                                <  1:3271103/93‑1 (MQ=255)
        ggCAACCAGAAGGCATCGCCCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTTTAAACCCTGATTGCCCCCCCCCCGtt                                                                                <  1:567565/92‑1 (MQ=255)
          cAACCAGAGGCCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGa                                                                              <  1:4651840/93‑1 (MQ=255)
           aaCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGag                                                                             >  2:6196534/1‑92 (MQ=255)
            aCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATg                                                                            <  2:903617/93‑1 (MQ=255)
            aCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATg                                                                            >  1:6365072/1‑93 (MQ=255)
            aCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATg                                                                            <  1:5540902/93‑1 (MQ=255)
            aCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATg                                                                            <  1:4612717/93‑1 (MQ=255)
            aCCAGAAGGCATCGCGCCCACCTGGACCGATTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCCttgttg                                                                            >  2:114945/1‑93 (MQ=255)
             ccAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGTTgg                                                                           >  2:4680696/1‑93 (MQ=255)
               aGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGt                                                                         <  2:4376255/93‑1 (MQ=255)
               aGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGt                                                                         <  2:3291631/93‑1 (MQ=255)
               aGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTa                                                                        >  2:1597197/1‑93 (MQ=255)
                gAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCCGTGATgggga                                                                        >  1:3516608/1‑91 (MQ=255)
                  aGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTaaa                                                                      >  1:3912236/1‑93 (MQ=255)
                    gCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAAt                                                                    >  2:5387997/1‑93 (MQ=255)
                     cATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAg                                                                  <  1:719437/93‑1 (MQ=255)
                       acgcgCCCACCTGCACCGTTTGTTCATTGCCGGAACTGGAACAGTCCCGGTTGAAATTCCTATTCCCCCCCCCCCCGTTGATGGGTAAAATAg                                                                  <  2:1636222/92‑1 (MQ=255)
                        cgcgCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGcc                                                                <  2:5749197/93‑1 (MQ=255)
                        cgcgCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGcc                                                                <  2:455903/93‑1 (MQ=255)
                           gcCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGGAAAATAGCCGCa                                                             >  1:5858149/1‑93 (MQ=255)
                             ccACCTGGACCTTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAAt                                                           <  1:5586956/93‑1 (MQ=255)
                             ccACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAAt                                                           <  1:2787393/93‑1 (MQ=255)
                              cACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAAt                                                           <  1:4908910/93‑1 (MQ=255)
                               cccTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAAttt                                                         <  2:3592462/92‑1 (MQ=255)
                               aCCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGTCCCGCCCCCGTTTAAGGGTGAAATATCCTCAAtat                                                         >  1:4087108/1‑91 (MQ=255)
                               aCCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAAttt                                                         >  1:2869863/1‑93 (MQ=255)
                               aCCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAAAAGCCGCAAttt                                                         >  2:3597045/1‑93 (MQ=255)
                                 cTGGCCCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAAttttt                                                       <  2:267488/93‑1 (MQ=255)
                                 cTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAAAAGCCGAAAttttt                                                       >  1:3892755/1‑93 (MQ=255)
                                 cTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAAtttttt                                                      >  2:498265/1‑92 (MQ=255)
                                  tGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTc                                                      <  2:809615/93‑1 (MQ=255)
                                   ggACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGCTTGGTAAAATAGGCGCCAttttttg                                                     >  2:1702852/1‑91 (MQ=255)
                                   ggACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCg                                                     <  1:441553/93‑1 (MQ=255)
                                   ggACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGGAAAATAGCCGCAAttttttg                                                     >  2:6102293/1‑91 (MQ=255)
                                    gACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGt                                                    <  2:4325740/93‑1 (MQ=255)
                                    gACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTAGATGGTTAAATTACCCCCAATTTTTTGt                                                    >  2:3253904/1‑93 (MQ=255)
                                    gACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGtt                                                   <  2:5538314/93‑1 (MQ=255)
                                     aCCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGGGACATAGCCGCAATTTTTCGtt                                                   >  2:5314691/1‑93 (MQ=255)
                                      ccGTTTCTTCGTTGCACGATCTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGttt                                                  <  2:5842943/93‑1 (MQ=255)
                                       cGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGCTTATGGGGAAAATAACCCCAATTTTTTGgttt                                                 >  2:298980/1‑93 (MQ=255)
                                         tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACGGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTca                                               <  1:4991117/93‑1 (MQ=255)
                                          ttcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTcaac                                             <  1:1597197/93‑1 (MQ=255)
                                           tcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAACCACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTcaac                                             <  1:2150123/93‑1 (MQ=255)
                                           tcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGACCCAATTTTTTGtttttcaa                                             >  2:3844467/1‑89 (MQ=255)
                                           tcttcGTGCCACGACCTGGAGCAGCTCCTTTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTcaac                                             <  1:2591002/93‑1 (MQ=255)
                                            cttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTcaaca                                            <  1:1522369/93‑1 (MQ=255)
                                              tcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGCTGATGGGTAAAATAGCCGCACTTTTTCGGTTTCAACAAgc                                         >  1:3123434/1‑93 (MQ=255)
                                               cGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgc                                         <  2:4380502/93‑1 (MQ=255)
                                               cGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgc                                         <  2:1477301/93‑1 (MQ=255)
                                                gTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCCTTGAGGGGGTTAAAAGCCGAAAGTTTTGGTGTTTATCAAacg                                        >  1:2317543/1‑93 (MQ=255)
                                                gTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGGTGTGGGGGAAAAAAACCCCAAATTTTTGCTTTttacaaacag                                       >  1:5825436/1‑89 (MQ=255)
                                                 tttcacaAACTGGGGCCGCCCCTGTGTAAACCCTCATTGCCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcg                                        <  1:247592/90‑1 (MQ=255)
                                                 tttcACGAACTGGGGCATCCCGTGTTTAAACCCTGATTGCCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcg                                        <  1:5933224/90‑1 (MQ=255)
                                                 tggCACGAACGGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgc                                       <  1:3745687/91‑1 (MQ=255)
                                                  tGCCCGAACTGGAGCAGCTCCTGTGTAAACCCTGATTGCCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgc                                       <  1:1039902/93‑1 (MQ=255)
                                                   gCACGAACGGGAGCAGCTCCTGTGTAAACACGGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcgg                                     <  1:842254/93‑1 (MQ=255)
                                                    cACGAAGTGGAGCACCTCTTTTGTAAACCCTGATTCCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggc                                    <  2:5497659/93‑1 (MQ=255)
                                                     aCGAACTGGAGCAGCTGCTGTGTAAACACTGATTGCCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggc                                    <  1:3314618/93‑1 (MQ=255)
                                                     aCGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGTTGGGGAAAAAAGCCGCAAATTTTTGTTTTCAACAAccgcgggc                                   >  2:3800548/1‑91 (MQ=255)
                                                     aCGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTTGTTTTCAACAAacgcggcg                                   >  2:607653/1‑93 (MQ=255)
                                                      cGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggcgc                                  <  2:6081535/93‑1 (MQ=255)
                                                       gAACTGGAGCAGCTCCTGTGTAAACACTGATTCCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggcgcg                                 <  2:173472/93‑1 (MQ=255)
                                                        aaCTGGAGCATCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGa                                <  1:162657/93‑1 (MQ=255)
                                                         aCTGGAGGAGCTCCTGTGTAAACGCTGATTGCCCCCCCCACGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGAt                               <  1:5450758/93‑1 (MQ=255)
                                                         aCTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGAt                               <  2:649338/93‑1 (MQ=255)
                                                         aCTGGAGCAGCTCCTGTGTAAACACTGATGGCCGCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATg                              <  1:2556259/93‑1 (MQ=255)
                                                           tGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGc                             <  1:2071191/93‑1 (MQ=255)
                                                           tGGAGCAGCTCCTGTGTAAACACTGATTCCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGc                             <  2:1653467/93‑1 (MQ=255)
                                                            ggAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGcc                            <  2:188487/93‑1 (MQ=255)
                                                            ggAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGcc                            >  1:3164564/1‑93 (MQ=255)
                                                            ggAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCCTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTTAACAAAAGCGGCGGGATGcc                            >  2:1520568/1‑93 (MQ=255)
                                                               gcagcTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCt                         <  2:3973910/93‑1 (MQ=255)
                                                               gcagcTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCtt                        >  1:6163846/1‑93 (MQ=255)
                                                                 agcTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTa                       >  1:2204720/1‑93 (MQ=255)
                                                                  gcTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTAc                      <  1:6028237/93‑1 (MQ=255)
                                                                   cTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTc                    <  1:6210914/93‑1 (MQ=255)
                                                                     ccTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCa                   >  2:6427850/1‑93 (MQ=255)
                                                                       tgtgTAAACACTGATTGCCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGGGATACCGCTTACTCaa                  >  1:1034356/1‑93 (MQ=255)
                                                                        gtgtAAACACTGATTGCCCCCCCCCCGGTGATGGGGGAAAAAACCGCAATTTTTTGCTTTTTAAAAACGCGGCGCGATGCGGCCTTCTCaaga                >  2:5644808/1‑93 (MQ=255)
                                                                         tgtAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaagaa               <  1:3949937/93‑1 (MQ=255)
                                                                         tgtAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaagaa               <  1:1696878/93‑1 (MQ=255)
                                                                          gtAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaagaag              <  1:2534751/93‑1 (MQ=255)
                                                                           tAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAagaaga             <  2:5033947/93‑1 (MQ=255)
                                                                           tAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAagaaga             >  1:4318846/1‑93 (MQ=255)
                                                                              acacTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGaagaaaga         >  2:2845400/1‑93 (MQ=255)
                                                                                 cTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAAt       <  2:1973281/93‑1 (MQ=255)
                                                                                  tGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAAtt      <  2:3985483/93‑1 (MQ=255)
                                                                                  tGATTGCCCCCCCCCCGTTGATGGGGAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAAtt      >  1:718179/1‑93 (MQ=255)
                                                                                   gATTGCCCCCCCCCCGTTGATGGGGAAAAATACCGCAATTTTTTGTTTTCCACAAGCACGGCGCGACTCCGCATGCTCCAAAAGAAAGAAATa     >  2:4728126/1‑93 (MQ=255)
                                                                                   gATTGCCCCGCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTAt    <  2:2294124/93‑1 (MQ=255)
                                                                                    aTTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTAt    <  1:5343127/93‑1 (MQ=255)
                                                                                      tGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAATAATTATGa  <  1:6286169/93‑1 (MQ=255)
                                                                                            |                                                                                        
CCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGA  >  NC_000913/4606808‑4606986

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