Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 691,865 C→A 6.6% intergenic (‑81/+9) ybeX ← / ← ybeY putative ion transport/ssRNA‑specific endoribonuclease; 16S rRNA 3' end maturation and quality control co‑endoribonuclease working with RNase R; rRNA transcription antitermination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913691,8650CA 6.6% 247.1 / 4.4 75intergenic (‑81/+9)ybeX/ybeYputative ion transport/ssRNA‑specific endoribonuclease; 16S rRNA 3' end maturation and quality control co‑endoribonuclease working with RNase R; rRNA transcription antitermination factor
Reads supporting (aligned to +/- strand):  ref base C (37/33);  new base A (2/3);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

GTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCA  >  NC_000913/691718‑692010
                                                                                                                                                   |                                                                                                                                                 
gTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGcc                                                                                                                                                  <  4:285814/149‑1 (MQ=255)
   aCAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTgg                                                                                                                                               >  2:366190/1‑149 (MQ=255)
      gggAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCa                                                                                                                                             >  2:344558/1‑148 (MQ=255)
        gagaAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAat                                                                                                                                          <  7:78815/149‑1 (MQ=255)
        gagaAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGGCATTTGACTGGCagac                                                                                                                                          <  5:375370/149‑4 (MQ=255)
            aaaaTCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAattatt                                                                                                                                      <  1:366190/149‑1 (MQ=255)
              aaTCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTAttct                                                                                                                                    <  3:158123/149‑1 (MQ=255)
               aTCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCt                                                                                                                                                 >  5:279813/1‑135 (MQ=255)
               aTCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCt                                                                                                                                                 <  6:279813/135‑1 (MQ=255)
                  ccTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTAttctttct                                                                                                                                <  1:344558/149‑1 (MQ=255)
                      cttGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCgg                                                                                                                             >  3:388178/1‑148 (MQ=255)
                         gttgCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCg                                                                                                                              <  5:319833/144‑1 (MQ=255)
                         gttgCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCg                                                                                                                              >  6:319833/1‑144 (MQ=255)
                            gCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTa                                                                                                                      >  4:395301/1‑149 (MQ=255)
                             cTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTa                                                                                                                      <  1:349849/148‑1 (MQ=255)
                                aTCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACgg                                                                                                                   >  7:239519/1‑148 (MQ=255)
                                 tCGTGTCACTACTGTGTAAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGAt                                                                                                                 <  8:149386/149‑1 (MQ=255)
                                   gtgtCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGCCGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATcc                                                                                                               >  5:417552/1‑149 (MQ=255)
                                   gtgtCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATcc                                                                                                               <  4:388178/149‑1 (MQ=255)
                                    tgtCACTACTGTTTGAATTGTCGTCGCTCATGGCGTCTTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCt                                                                                                              <  6:252014/149‑1 (MQ=255)
                                    tgtCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCt                                                                                                              <  2:366621/149‑1 (MQ=255)
                                        actactGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCa                                                                                                            >  2:122224/1‑147 (MQ=255)
                                               tgttAATTGTCGTCGCTCATGGCGTCGTATGGGTGCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTAttctttc                                                                                                                                 <  3:257214/119‑1 (MQ=255)
                                               tgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTAttctttc                                                                                                                                 >  4:257214/1‑119 (MQ=255)
                                               tgtgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCa                                                                                                     <  1:342312/147‑1 (MQ=255)
                                                 tgAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCagag                                                                                                  <  2:260140/148‑1 (MQ=255)
                                                        tcgtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATcc                                                                                                               >  5:332943/1‑128 (MQ=255)
                                                        tcgtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATcc                                                                                                               <  6:332943/128‑1 (MQ=255)
                                                        tcgtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCAt                                                                                           >  8:146894/1‑148 (MQ=255)
                                                        tcgtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATa                                                                                          <  8:239519/149‑1 (MQ=255)
                                                           tcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTGTCTCGGCAATGTACGGATACTCATCGCACAGAGCAAGCATaagc                                                                                       >  6:139841/1‑147 (MQ=255)
                                                            cgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAAtct                                                                                      <  3:395301/149‑1 (MQ=255)
                                                              cTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCAAAATctct                                                                                    <  2:183703/149‑1 (MQ=255)
                                                               tCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTg                                                                                   >  8:83048/1‑149 (MQ=255)
                                                                cATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGt                                                                                  <  8:143873/149‑1 (MQ=255)
                                                                 aTGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTTTTCTTTCTCGGAAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGtt                                                                                 <  7:415340/149‑1 (MQ=255)
                                                                    gCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTc                                                                               <  1:122224/148‑1 (MQ=255)
                                                                      gtcgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAg                                                                            <  5:153150/149‑1 (MQ=255)
                                                                              ggTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCa                                                                    >  7:458870/1‑149 (MQ=255)
                                                                               gTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCAt                                                                   >  5:233149/1‑149 (MQ=255)
                                                                               gTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCAAg                                                                            <  8:190953/140‑1 (MQ=255)
                                                                               gTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCAAg                                                                            >  7:190953/1‑140 (MQ=255)
                                                                                     aTGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAggg                                                                          >  7:106041/1‑136 (MQ=255)
                                                                                     aTGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAggg                                                                          <  8:106041/136‑1 (MQ=255)
                                                                                     aTGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATttcttc                                                             >  1:125027/1‑149 (MQ=255)
                                                                                     aTGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATttctt                                                              <  8:458870/148‑1 (MQ=255)
                                                                                     aTGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATttct                                                               <  7:83048/147‑1 (MQ=255)
                                                                                         tagttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCt                                                         >  8:56681/1‑149 (MQ=255)
                                                                                          agttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCtt                                                        >  2:104713/1‑149 (MQ=255)
                                                                                            ttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGagcaagca                                                                                            <  7:60125/111‑1 (MQ=255)
                                                                                            ttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGagcaagca                                                                                            >  8:60125/1‑111 (MQ=255)
                                                                                            ttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCg                                                      >  4:377156/1‑149 (MQ=255)
                                                                                            ttaATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCg                                                      >  2:239329/1‑149 (MQ=255)
                                                                                             taATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGt                                                     >  5:136992/1‑149 (MQ=255)
                                                                                                       gCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGAt                                           <  3:377156/149‑1 (MQ=255)
                                                                                                                 aTCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTc                                                                      <  4:193589/112‑1 (MQ=255)
                                                                                                                 aTCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTATGTTTCGAGGGATTc                                                                      >  3:193589/1‑112 (MQ=255)
                                                                                                                      cAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTAc                                                                                                                     <  6:360180/60‑1 (MQ=255)
                                                                                                                      cAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTAc                                                                                                                     >  5:360180/1‑60 (MQ=255)
                                                                                                                        aCGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAAc                           >  4:345403/1‑148 (MQ=255)
                                                                                                                         cggcggGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACa                          >  4:356957/1‑148 (MQ=255)
                                                                                                                           gcggGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAAt                       >  1:220437/1‑149 (MQ=255)
                                                                                                                              gggACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCa                    <  2:125027/149‑1 (MQ=255)
                                                                                                                                 taatttGCCAGGCATTTGACTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGAt                                           >  3:139558/6‑123 (MQ=255)
                                                                                                                                 taatttGCCAGGCATTTGACTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGAt                                           <  4:139558/118‑1 (MQ=255)
                                                                                                                                   gTCTGCCAGTCAAATGCCTGGCAAATTATTCTATCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCAATTTCAACTGCTTCGTCATATTCGATGTGATCGTAACCTAACAAATGCAGACTg               >  7:426742/1‑149 (MQ=255)
                                                                                                                                    tttGCCAGGCATTTGACTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTc                                                       <  6:345808/106‑1 (MQ=255)
                                                                                                                                    tttGCCAGGCATTTGACTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTc                                                       >  5:345808/3‑108 (MQ=255)
                                                                                                                                    tCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATCGCCCAGCGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGc              >  8:235148/1‑149 (MQ=255)
                                                                                                                                       gCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGt           >  5:282416/1‑149 (MQ=255)
                                                                                                                                       gCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGt           >  6:22906/1‑149 (MQ=255)
                                                                                                                                        ccAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAAcc                               >  7:5796/1‑128 (MQ=255)
                                                                                                                                        ccAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAAcc                               <  8:5796/128‑1 (MQ=255)
                                                                                                                                            tCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTAcacc      <  1:406914/149‑1 (MQ=255)
                                                                                                                                              aaaTGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCAATTTATTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCTTGcaccac    >  5:108167/1‑149 (MQ=255)
                                                                                                                                                aTGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCAATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCAcaacca  >  6:347231/1‑149 (MQ=255)
                                                                                                                                                   ccTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCa                    >  1:120713/1‑128 (MQ=255)
                                                                                                                                                   ccTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTACATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCa                    <  2:120713/128‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACGGATCCTCATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCA  >  NC_000913/691718‑692010

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