Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,606,898 Δ1 bp 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,8980T.79.9% 638.2 / 236.2 277intergenic (+6/‑85)yjjG/prfCputative phosphatase/peptide chain release factor RF‑3
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base . (77/122);  total (111/166)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 285 total aligned reads displayed.

CCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGA  >  NC_000913/4606808‑4606986
                                                                                          |                                                                                        
ccATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCcccc                                                                                        >  1:367282/1‑93 (MQ=255)
  aTCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTcccc                                                                                      >  2:2692602/1‑93 (MQ=255)
   tCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGcccccccc                                                                                     >  2:1132506/1‑93 (MQ=255)
       gAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGtt                                                                                <  1:6425340/93‑1 (MQ=255)
       gAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGtt                                                                                <  1:7178925/93‑1 (MQ=255)
         gCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTg                                                                               >  2:5787617/1‑93 (MQ=255)
         gCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGa                                                                              <  1:1412380/93‑1 (MQ=255)
         gCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGa                                                                              <  2:1377203/93‑1 (MQ=255)
          cAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGAt                                                                             <  1:1477158/93‑1 (MQ=255)
            aCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATg                                                                            >  1:6032320/1‑93 (MQ=255)
            aCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATg                                                                            >  1:6761831/1‑93 (MQ=255)
            aCCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATgg                                                                           <  1:4021617/93‑1 (MQ=255)
             ccAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATggg                                                                          <  1:1555018/93‑1 (MQ=255)
              cAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATggg                                                                          <  1:749788/93‑1 (MQ=255)
              cAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGC‑‑CCCCCCCGTTGATGGGTa                                                                        <  1:235829/93‑1 (MQ=255)
               aGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTa                                                                        >  2:5428758/1‑93 (MQ=255)
                gAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTaa                                                                       <  2:297227/93‑1 (MQ=255)
                  aGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTaaa                                                                      >  1:4909966/1‑93 (MQ=255)
                  aGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTaaaa                                                                     <  2:3263962/93‑1 (MQ=255)
                   ggCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAAt                                                                    <  2:2519741/93‑1 (MQ=255)
                     cATCGCGCCCACCTGGACCGTTTCTTGGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAg                                                                  <  2:737053/93‑1 (MQ=255)
                      aTCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGc                                                                 <  2:4923615/93‑1 (MQ=255)
                      aTCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGc                                                                 >  2:447668/1‑93 (MQ=255)
                       tCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGcc                                                                <  2:1840658/93‑1 (MQ=255)
                            cccACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCaa                                                            <  1:1950634/93‑1 (MQ=255)
                             ccACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAAtt                                                          >  2:7039273/1‑93 (MQ=255)
                             ccACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGAAGGGTAAAATAGCCGCAAtt                                                          <  2:5478862/93‑1 (MQ=255)
                              cACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCCGTTTTGGGTGAAAAAACCGccaat                                                          >  2:1273366/1‑89 (MQ=255)
                              cACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGC‑‑CCCCCCCGTTGATGGGTAAAATAGCCGCAAtttt                                                        >  2:6829127/1‑93 (MQ=255)
                               cccTGGGCCCTTTCTTCGTTTGCCGAACCGGGGCCGCCTCCGTGTAAAACCTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAAttt                                                         <  1:3724657/92‑1 (MQ=255)
                               aCCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAAtttt                                                        >  2:908877/1‑93 (MQ=255)
                                 cTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTc                                                      <  2:6278008/93‑1 (MQ=255)
                                  tGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCg                                                     >  1:3104798/1‑93 (MQ=255)
                                   ggACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCg                                                     <  1:1649202/93‑1 (MQ=255)
                                   ggACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGGAAAATAGCCGCAATTTTTTGt                                                    >  1:2639894/1‑93 (MQ=255)
                                    gACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGt                                                    >  2:191847/1‑93 (MQ=255)
                                    gACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGtt                                                   <  2:310894/93‑1 (MQ=255)
                                    gACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGtt                                                   <  1:4406078/93‑1 (MQ=255)
                                    gACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATCGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGtt                                                   <  1:5122707/93‑1 (MQ=255)
                                     aCCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGttt                                                  >  2:3397634/1‑93 (MQ=255)
                                      ccGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGttt                                                  <  1:81369/93‑1 (MQ=255)
                                      ccGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGtttt                                                 <  1:919205/93‑1 (MQ=255)
                                       cGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGtttt                                                 >  2:1732309/1‑93 (MQ=255)
                                       cGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTc                                                <  2:4015124/93‑1 (MQ=255)
                                       cGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTc                                                >  1:5034021/1‑93 (MQ=255)
                                       cGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGC‑‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTca                                               <  2:7082078/93‑1 (MQ=255)
                                        gTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTc                                                <  1:5420510/93‑1 (MQ=255)
                                        gTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTca                                               <  2:6074616/93‑1 (MQ=255)
                                        gTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGGAAAATAGCCGCAATTTTTTGttttta                                               >  1:3889659/1‑91 (MQ=255)
                                         tttCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTcaa                                              >  1:6536598/1‑93 (MQ=255)
                                          ttcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTcaac                                             >  1:1070735/1‑93 (MQ=255)
                                           tcttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTcaaca                                            <  2:812297/93‑1 (MQ=255)
                                            cttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTcaacaa                                           <  2:5090986/93‑1 (MQ=255)
                                             ttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCCCTTGATGGGTAAAATAGCCCCAATTTTTCGTTTTtaaaaa                                           >  2:3179583/1‑93 (MQ=255)
                                             ttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAg                                          <  1:3232655/93‑1 (MQ=255)
                                             ttcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGC‑‑CCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgc                                         >  1:2764087/1‑93 (MQ=255)
                                              tcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgc                                         <  2:5220922/93‑1 (MQ=255)
                                              tcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgc                                         <  1:4358917/93‑1 (MQ=255)
                                              tcGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgc                                         <  1:705104/93‑1 (MQ=255)
                                               cGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcg                                        >  1:5563931/1‑93 (MQ=255)
                                                gTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgc                                       <  1:3451715/93‑1 (MQ=255)
                                                gTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgc                                       <  2:3526222/93‑1 (MQ=255)
                                                   gCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggc                                    <  2:213826/93‑1 (MQ=255)
                                                   gCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAgcgcggc                                    <  2:4321971/93‑1 (MQ=255)
                                                       gAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGa                                <  2:4710441/93‑1 (MQ=255)
                                                       gAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGa                                >  2:3711519/1‑93 (MQ=255)
                                                        aaCTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGAt                               <  1:452663/93‑1 (MQ=255)
                                                        aaCTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGAt                               <  2:3634954/93‑1 (MQ=255)
                                                          cTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGc                             <  2:5684536/93‑1 (MQ=255)
                                                          cTGGAGCAGCTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGc                             <  2:2443171/93‑1 (MQ=255)
                                                              agcagcTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCt                         <  1:3923382/93‑1 (MQ=255)
                                                              agcagcTCATGTGTAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGc                          <  2:166169/93‑1 (MQ=255)
                                                                 agcTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTAc                      <  2:5962085/93‑1 (MQ=255)
                                                                  gcTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACt                     >  1:5860334/1‑93 (MQ=255)
                                                                   cTCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCttcatc                    >  1:7094198/1‑89 (MQ=255)
                                                                    tCCTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCa                   <  2:4158940/93‑1 (MQ=255)
                                                                     ccTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCa                   <  1:4639817/93‑1 (MQ=255)
                                                                     ccTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaa                  <  2:1402197/93‑1 (MQ=255)
                                                                     ccTGTGTAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaa                  <  2:6748870/93‑1 (MQ=255)
                                                                       tgtgTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaag                 <  2:440669/93‑1 (MQ=255)
                                                                        gtgtAAACACTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaaga                <  2:4212025/93‑1 (MQ=255)
                                                                        gtgtAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaagaa               <  2:2842599/93‑1 (MQ=255)
                                                                        gtgtAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCaagaa               <  1:4785245/93‑1 (MQ=255)
                                                                          gtAAACACTGATTGCCCCCCCCCGTTGATGGGGAAAATAGCCGCAATTTTTTGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTTAagaaga             >  2:2549496/1‑93 (MQ=255)
                                                                           tAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAagaaga             <  2:6252139/93‑1 (MQ=255)
                                                                           tAAACACTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAgaagaa            >  2:7236319/1‑93 (MQ=255)
                                                                              acacTGATTGCCCCCCCTCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGaagaaaga         <  2:6491310/93‑1 (MQ=255)
                                                                                 cTGATTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAAt       <  1:1000412/93‑1 (MQ=255)
                                                                                 cTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAAtt      <  2:6621327/93‑1 (MQ=255)
                                                                                 cTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAAtt      <  1:1101684/93‑1 (MQ=255)
                                                                                 cTGATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAAtt      >  2:750652/1‑93 (MQ=255)
                                                                                 cTGATTCCCCCCCCCCGTTGATAGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAAtt      <  1:6339059/93‑1 (MQ=255)
                                                                                  tGATTGCCCCCCCCCGTTGATGAGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTa     <  1:54000/93‑1 (MQ=255)
                                                                                   taTTGCCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTa     <  1:3560202/92‑1 (MQ=255)
                                                                                   gATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTa     <  2:3931869/93‑1 (MQ=255)
                                                                                   gATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTAt    >  2:6445156/1‑93 (MQ=255)
                                                                                   gATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTAt    <  2:5564930/93‑1 (MQ=255)
                                                                                   gATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTAt    <  1:7400294/93‑1 (MQ=255)
                                                                                   gATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTAt    <  2:4354535/93‑1 (MQ=255)
                                                                                   gATTGCCCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTAt    <  2:4890469/93‑1 (MQ=255)
                                                                                          |                                                                                        
CCATCGCGAGCAACCAGAAGGCATCGCGCCCACCTGGACCGTTTCTTCGTTGCACGAACTGGAGCAGCTCCTGTGTAAACACTGATTGCCTCCCCCCCGTTGATGGGTAAAATAGCCGCAATTTTTCGTTTTCAACAAGCGCGGCGCGATGCCGCTTACTCAAGAAGAAAGAATTATGA  >  NC_000913/4606808‑4606986

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