Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,484,494 T→C 6.1% V486V (GTA→GTG)  pepA ← multifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,484,4940TC 6.1% 210.6 / 4.0 65V486V (GTA→GTG) pepAmultifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base C (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.44e-01

TGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGC  >  NC_000913/4484347‑4484640
                                                                                                                                                   |                                                                                                                                                  
tgttgtGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTa                                                                                                                                                   >  2:593921/1‑149 (MQ=255)
  ttgtGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCCAcc                                                                                                                                                 >  3:562152/1‑149 (MQ=255)
       tCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGgacg                                                                                                                                            >  1:224873/1‑149 (MQ=255)
         ttATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGgacgac                                                                                                                                          >  2:148578/1‑149 (MQ=255)
           aTCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCg                                                                                                                                        >  2:407412/1‑149 (MQ=255)
                   tACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGcc                                                                                                                                <  3:6655/149‑1 (MQ=255)
                   tACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGcc                                                                                                                                >  2:530902/1‑149 (MQ=255)
                   tACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCCACCGGTCGTCCGGTAgcgtt                                                                                                                                <  4:562152/149‑3 (MQ=39)
                   tACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCCACCGGTCGTCCGGTAgcg                                                                                                                                  <  2:451415/147‑1 (MQ=39)
                   tACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCCACCGGTCGTCCGGTAgcg                                                                                                                                  >  1:451415/1‑147 (MQ=39)
                        gTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTg                                                                                                                           >  3:499509/1‑149 (MQ=255)
                          tCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCt                                                                                                                         <  1:148578/149‑1 (MQ=255)
                               gCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTAcc                                                                                                                    <  2:284390/149‑1 (MQ=255)
                                cTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCtt                                                                                                                        >  3:360473/1‑144 (MQ=255)
                                cTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCtt                                                                                                                        <  4:360473/144‑1 (MQ=255)
                                   tgtAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGaa                                                                                                                >  4:397048/1‑149 (MQ=255)
                                        gTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCa                                                                                                           <  3:558097/149‑1 (MQ=255)
                                              tAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAGCAGGTACTGTGCCGGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGt                                                                                                     >  3:552922/1‑149 (MQ=255)
                                                aGGCGTTCACGCCGCATCCGGGAATAACGGCCTTGCCTGACGCAGTTACTCTTCGACGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTAc                                                                                                   >  2:489640/1‑149 (MQ=255)
                                                    gTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGc                                                                                               >  4:583299/1‑149 (MQ=255)
                                                    gTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGc                                                                                               >  1:155633/1‑149 (MQ=255)
                                                    gTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGc                                                                                               >  2:146477/1‑149 (MQ=255)
                                                       cACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGat                                                                                            <  3:47725/149‑1 (MQ=255)
                                                        aCGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGata                                                                                           >  4:65458/1‑149 (MQ=255)
                                                           ccgcATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATcc                                                                                        <  3:150830/149‑1 (MQ=255)
                                                            cgcATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCa                                                                                       <  1:407412/149‑1 (MQ=255)
                                                             gcATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAg                                                                                      <  4:167667/149‑1 (MQ=255)
                                                             gcATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAg                                                                                      <  1:220460/149‑1 (MQ=255)
                                                             gcATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAg                                                                                      <  3:549005/149‑1 (MQ=255)
                                                                 ccGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTgc                                                                                  <  1:510431/149‑1 (MQ=255)
                                                                   ggCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTgcgc                                                                                >  3:508772/1‑149 (MQ=255)
                                                                   ggCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTgcgc                                                                                >  2:568396/1‑149 (MQ=255)
                                                                   ggCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTgcgc                                                                                <  3:331802/149‑1 (MQ=255)
                                                                   ggCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCACCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTgcgc                                                                                >  1:230577/1‑149 (MQ=255)
                                                                     cAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGccc                                                                              >  1:349186/1‑149 (MQ=255)
                                                                      aaTAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCa                                                                             <  2:347065/149‑1 (MQ=255)
                                                                      aaTAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCa                                                                             >  1:517593/1‑149 (MQ=255)
                                                                      aaTAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCa                                                                             <  4:508772/149‑1 (MQ=255)
                                                                                 tGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTAc                                                                  <  2:471706/149‑1 (MQ=255)
                                                                                  gCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTa                                                                   <  3:305554/147‑1 (MQ=255)
                                                                                  gCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTa                                                                   >  4:305554/1‑147 (MQ=255)
                                                                                   ccTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACgg                                                                <  3:397048/149‑1 (MQ=255)
                                                                                      gACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTa                                                             >  1:261017/1‑149 (MQ=255)
                                                                                        cGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTaaa                                                           <  1:435224/149‑1 (MQ=255)
                                                                                         gCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAg                                                          >  1:398459/1‑149 (MQ=255)
                                                                                             ttACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTg                                                      >  4:482626/1‑149 (MQ=255)
                                                                                             ttACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTg                                                      >  3:297066/1‑149 (MQ=255)
                                                                                                ctctTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACa                                                   <  3:450903/149‑1 (MQ=255)
                                                                                                ctctTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACa                                                   <  4:330780/149‑1 (MQ=255)
                                                                                                 tctTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAg                                                  >  4:156684/1‑149 (MQ=255)
                                                                                                   ttCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCC‑TTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGaa                                               >  1:41816/1‑149 (MQ=255)
                                                                                                    tCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGaa                                               <  4:54243/149‑1 (MQ=255)
                                                                                                    tCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGaa                                               <  1:530902/149‑1 (MQ=255)
                                                                                                        cGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGc                                                         <  1:431112/135‑1 (MQ=255)
                                                                                                        cGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAAAGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGc                                                         >  2:431112/1‑135 (MQ=255)
                                                                                                                             aCAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCccacca                      <  1:433139/149‑1 (MQ=255)
                                                                                                                              cAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAg                     <  1:96384/149‑1 (MQ=255)
                                                                                                                               aGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTAc                                                                  >  2:229098/1‑103 (MQ=255)
                                                                                                                               aGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTAc                                                                  <  1:229098/103‑1 (MQ=255)
                                                                                                                                 gAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGgac                  >  1:441619/1‑149 (MQ=255)
                                                                                                                                    cTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGgacgac               <  2:398459/149‑1 (MQ=255)
                                                                                                                                        gCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGcc           <  1:568396/149‑1 (MQ=255)
                                                                                                                                           aGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAAt        >  4:542457/1‑149 (MQ=255)
                                                                                                                                               aCGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTc    <  3:65458/149‑1 (MQ=255)
                                                                                                                                                 gCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTAACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGc  >  3:339504/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGTTGTGTCTTATCAGGCCTACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGC  >  NC_000913/4484347‑4484640

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