Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,484,509 G→A 6.1% A481A (GCC→GCT)  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,5090GA 6.1% 208.6 / 3.5 65A481A (GCC→GCT) pepAmultifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base G (30/31);  new base A (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

TACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATT  >  NC_000913/4484366‑4484655
                                                                                                                                               |                                                                                                                                                  
tACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGcc                                                                                                                                               <  3:6655/149‑1 (MQ=255)
tACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGcc                                                                                                                                               >  2:530902/1‑149 (MQ=255)
tACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCCACCGGTCGTCCGGTAgcgtt                                                                                                                                               <  4:562152/149‑3 (MQ=39)
tACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCCACCGGTCGTCCGGTAgcg                                                                                                                                                 >  1:451415/1‑147 (MQ=39)
tACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCCACCGGTCGTCCGGTAgcg                                                                                                                                                 <  2:451415/147‑1 (MQ=39)
     gTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTg                                                                                                                                          >  3:499509/1‑149 (MQ=255)
       tCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCt                                                                                                                                        <  1:148578/149‑1 (MQ=255)
            gCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTAcc                                                                                                                                   <  2:284390/149‑1 (MQ=255)
             cTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCtt                                                                                                                                       <  4:360473/144‑1 (MQ=255)
             cTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCtt                                                                                                                                       >  3:360473/1‑144 (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                                                                                                              >  1:155633/1‑149 (MQ=255)
                                 gTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGc                                                                                                              >  4:583299/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                                                                                                     <  1:220460/149‑1 (MQ=255)
                                          gcATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAg                                                                                                     <  4:167667/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                                                                                               <  3:331802/149‑1 (MQ=255)
                                                ggCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTgcgc                                                                                               >  2:568396/1‑149 (MQ=255)
                                                ggCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCACCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTgcgc                                                                                               >  1:230577/1‑149 (MQ=255)
                                                  cAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGccc                                                                                             >  1:349186/1‑149 (MQ=255)
                                                   aaTAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCa                                                                                            >  1:517593/1‑149 (MQ=255)
                                                   aaTAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCa                                                                                            <  4:508772/149‑1 (MQ=255)
                                                   aaTAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCa                                                                                            <  2:347065/149‑1 (MQ=255)
                                                              tGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTAc                                                                                 <  2:471706/149‑1 (MQ=255)
                                                               gCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTa                                                                                  >  4:305554/1‑147 (MQ=255)
                                                               gCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTa                                                                                  <  3:305554/147‑1 (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)
                                                                                                                                   cGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCatat            >  2:275317/1‑149 (MQ=255)
                                                                                                                                     gtcgtccGGTAGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCg          >  3:303690/6‑149 (MQ=255)
                                                                                                                                         acCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCaa      >  4:410377/1‑149 (MQ=255)
                                                                                                                                          ccGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCaaa     >  1:478950/1‑149 (MQ=255)
                                                                                                                                             gTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAAtt  <  1:224512/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TACGAGTTCAGTGCTGTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGGCAATAACAGCCTTGCCTGACGCAATTACTCTTCGCCGTTAAACCCAGCGCGGTTTAACAGGAACTGTGCCAGCAACGCTACCGGACGACCGGTGGCGCCTTTTGCTTTACCAGAACGCCAGGCGGTACCGGCGATATCCAGGTGCGCCCAGTTGTACTTACGGGTAAAGCGTGACAGGAAGCAACCTGCGGTAATCGCCCCACCAGGACGACCGCCAATGTTCGCCATATCGGCAAAATT  >  NC_000913/4484366‑4484655

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: ‑