Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,506,444 C→G 5.6% A63A (GCG→GCC ypdF ← Xaa‑Pro aminopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,506,4440CG 5.6% 291.6 / 4.7 86A63A (GCG→GCCypdFXaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand):  ref base C (39/42);  new base G (3/2);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGACTGCCAGCGGTGCGCGGTTTCCCAACTCACCTGCTGGCCCTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCACTGCATCCAGTTGT  >  NC_000913/2506301‑2506591
                                                                                                                                               |                                                                                                                                                   
tCAGACTGCCAGCGGTGCGCGGTTTCCCAACTCACCTGCTGGCCCTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCAc                                                                                                                                                 >  5:360975/1‑148 (MQ=255)
  aGACTGCCCGCGGTGCGCGGGTTCCCAACTCACCTGCTGGCCCTCAAACCCGAGCGTGTGCAACTGTTCATCGGCACTGATTTGCTTGACGATAGTGGTAAGCGTGTGCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACg                                                                                                                                                <  3:389885/147‑1 (MQ=255)
  aGACTGCCAGCGGTGCGCGGTTTCCCAACTCACCTGCTGGCCCTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACgg                                                                                                                                               >  1:276162/1‑148 (MQ=255)
        ccAGCGGTGCGCGGTTTCCCAACTCACCTGCTGGCCCTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTcc                                                                                                                                         <  3:159237/148‑1 (MQ=255)
        ccAGCGGTGCGCGGTTTCCCAACTCACCTGCTGGCCCTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTc                                                                                                                                          <  1:429749/147‑1 (MQ=255)
            cGGTGCGCGGTTTCCCAACTCACCTGCTGGCCCTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCt                                                                                                                                     >  4:76105/1‑148 (MQ=255)
                  gcgGTTTCCCAACTCACCTGCTGGCCGTCAAAACCGAGCGTTTGCAACTGTTCTGCGGAAATGATTTGATTGACGATTGTGGTACGCGGGTTGGTGGCGGCAAGCAAAGGCGGCGGGTAGCCTTGGGCCCGGGGTGCCACCTCACCg                                                                                                                                >  4:7143/1‑147 (MQ=255)
                  gcgGTTTCCCAACTCACCTGCTGGCCCTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACg                                                                                                                                >  5:63956/1‑147 (MQ=255)
                     gTTTCCCAACTCACCTGCTGGCCCTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAAt                                                                                                                            <  7:104486/148‑1 (MQ=255)
                           cAACTCACCTGCTGGCCCTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAgcgc                                                                                                                       >  3:349355/1‑147 (MQ=255)
                            aaCTCACCTGCTGGCCCTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGa                                                                                                                     <  3:110307/148‑1 (MQ=255)
                             aCTCACCTGCTGGCCCTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGaa                                                                                                                    >  4:261570/1‑148 (MQ=255)
                             aCTCACCTGCTGGCCCTCAAAACCGAGCGTTTGCAACTGATCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGGAATAGcgcgca                                                                                                                    >  6:122876/1‑146 (MQ=255)
                                   gtggtgtcccTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGAGTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATccacc                                                                                                              <  3:129746/141‑1 (MQ=255)
                                     gctgGCCCTCAAAACCGAGCGTTTGCACATGTTCATCGGCATTGATTTGATTGCCGATTGTGGTATGGGTGTTCGGCGCGTCAAGCAAATGCCGATGGTTGCGTTGCGCACGGGGTTCCACCGCCCGGTAATGGCGCGCTTccac                                                                                                               >  2:210316/1‑145 (MQ=255)
                                         gCCCTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAAtg                                                                                                        <  2:276162/148‑1 (MQ=255)
                                          cccTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGGGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAAtgt                                                                                                       >  7:157090/1‑148 (MQ=255)
                                          cccTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGGGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAAtgt                                                                                                       >  7:163121/1‑148 (MQ=255)
                                          cccTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAAtgt                                                                                                       >  2:400326/1‑148 (MQ=255)
                                           ccTGAAATCCGAGGGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCACGGAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAAtgtg                                                                                                      <  8:222651/148‑1 (MQ=255)
                                            cTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTgc                                                                                                     >  7:311268/1‑148 (MQ=255)
                                              cAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATcca                                                                                                                >  2:380655/1‑135 (MQ=255)
                                              cAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATcca                                                                                                                <  1:380655/135‑1 (MQ=255)
                                              cAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTgcg                                                                                                    >  8:26061/1‑147 (MQ=255)
                                              cAAAACCGAGAGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGa                                                                                                                     >  2:268103/1‑130 (MQ=255)
                                              cAAAACCGAGAGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGa                                                                                                                     <  1:268103/130‑1 (MQ=255)
                                               aaaaCCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCa                                                                                                  <  3:76105/148‑1 (MQ=255)
                                                aaaCCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCAc                                                                                                 >  2:17745/1‑148 (MQ=255)
                                                  aCCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAAtg                                                                                                        >  6:202380/1‑139 (MQ=255)
                                                  aCCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAAtg                                                                                                        <  5:202380/139‑1 (MQ=255)
                                                  aCCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACtt                                                                                               <  6:63956/148‑1 (MQ=255)
                                                    cGAGCGTTTGCAACTGTTCATCGGCAATGATTTGAGTGACGATAGTGGTAAGCTTGTTCGTCGCGTCAAGCATATGCAGCTGGTAGCCTTGCGGACTGGGTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCTCTTTc                                                                                             <  8:279371/148‑1 (MQ=255)
                                                     gAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTcc                                                                                                                                         <  8:391830/103‑1 (MQ=255)
                                                     gAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTcc                                                                                                                                         >  7:391830/1‑103 (MQ=255)
                                                       gCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACg                                                                                          <  6:183439/148‑1 (MQ=255)
                                                          tttGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGc                                                                                        >  2:435663/1‑147 (MQ=255)
                                                          tttGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCt                                                                                       <  6:360975/148‑1 (MQ=255)
                                                           ttGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCa                                                                                                  <  6:108863/136‑1 (MQ=255)
                                                           ttGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCa                                                                                                  >  5:108863/1‑136 (MQ=255)
                                                                aCTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCCCGGGATTACACCTCAACGTAATAGCGCGAATCCACCAGAATATGCGCACTTTCACGGCTAAtacc                                                                                 >  1:325772/1‑145 (MQ=255)
                                                                aCTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGGGCGCACTTTCACGGCTAATcac                                                                                 >  6:125466/1‑148 (MQ=255)
                                                                    ttCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCCCcaca                                                                             <  4:187770/148‑1 (MQ=255)
                                                                    ttCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACcaca                                                                             >  6:431704/1‑148 (MQ=255)
                                                                       aTCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATaa                                                                          <  1:400326/148‑1 (MQ=255)
                                                                          ggCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGAAGCTGGTAGCCTTGCGCACTGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAAcc                                                                        <  8:311268/147‑1 (MQ=255)
                                                                            cAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGa                                                                     >  3:116504/1‑148 (MQ=255)
                                                                             aaTGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGAt                                                                    <  5:431704/148‑1 (MQ=255)
                                                                              aTGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCCACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATc                                                                   >  5:187339/1‑148 (MQ=255)
                                                                              aTGATTTGAATGACGATAGTGGTAACCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAAGCCACCAGAATCTGCGCACTTTCACCGCAAATCACCACATAACCTGCTc                                                                   >  8:114734/1‑148 (MQ=255)
                                                                                gatttgattGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGa                                                                     <  4:279694/144‑1 (MQ=255)
                                                                                gatttgattGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGa                                                                     >  3:279694/1‑144 (MQ=255)
                                                                                 atttgattGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGTCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCg                                                                 <  2:205388/147‑1 (MQ=255)
                                                                                  tttgattGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGt                                                               <  3:392226/148‑1 (MQ=255)
                                                                                  tttgattGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGt                                                               <  7:26061/148‑1 (MQ=255)
                                                                                   ttgattgaCGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGt                                                               <  6:223348/147‑1 (MQ=255)
                                                                                   ttgattgaCGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTg                                                              >  3:193415/1‑148 (MQ=255)
                                                                                          aCGATAGTGGTCAGCGTGCTCGTCGCGTCAAGCAAAGGCAGCTGGCAGCCGTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATCACCTGATCCGGTGGAGATc                                                        <  4:349355/147‑1 (MQ=255)
                                                                                          aCGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATcc                                                       <  3:261570/148‑1 (MQ=255)
                                                                                           cGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATca                                                                                  >  8:319809/1‑120 (MQ=255)
                                                                                           cGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATca                                                                                  <  7:319809/120‑1 (MQ=255)
                                                                                            gATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATcccc                                                     <  8:177373/148‑1 (MQ=255)
                                                                                              tAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCgg                                                                >  1:188510/1‑135 (MQ=255)
                                                                                              tAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCgg                                                                <  2:188510/135‑1 (MQ=255)
                                                                                                 tGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGt                                                               <  4:172179/133‑1 (MQ=255)
                                                                                                 tGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGt                                                               >  3:172179/1‑133 (MQ=255)
                                                                                                     aaGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGc                                            >  6:225984/1‑148 (MQ=255)
                                                                                                          tgtTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGGTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCtgttt                                       <  1:8435/148‑1 (MQ=255)
                                                                                                                   cgTCAAGCAAATGCAGCTGGTAGCCTTGGGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCg                              <  8:163121/148‑1 (MQ=255)
                                                                                                                   cgTCAAGCAAATGCAGCTGGTAGCCTTGGGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCg                              <  8:157090/148‑1 (MQ=255)
                                                                                                                     tCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTg                            <  1:17745/148‑1 (MQ=255)
                                                                                                                     tCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTg                            <  1:435663/148‑1 (MQ=255)
                                                                                                                        aGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCa                                                                                            >  6:154793/1‑81 (MQ=255)
                                                                                                                        aGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCa                                                                                            <  5:154793/81‑1 (MQ=255)
                                                                                                                        aGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGGAATAGCGCGAATCCACCAGAATGTGCGCACGTTCACGGCTACTCACCACATAACCTGATCCGGTGGAGATCCCCATATGCGGCTGTTTGTTCTGCCGTgag                          >  5:117728/1‑147 (MQ=255)
                                                                                                                          cAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATGACCACATAACCTGATCCGGTGGAGATGCCCAGATGCGGCTGTTTGTTCTGCCGTgaggag                       <  7:114734/148‑1 (MQ=255)
                                                                                                                          cAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCATATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTgag                          >  8:288751/1‑145 (MQ=255)
                                                                                                                            aaTGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCtg                                          <  8:121/127‑1 (MQ=255)
                                                                                                                            aaTGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCtg                                          >  7:121/1‑127 (MQ=255)
                                                                                                                             aTGAAGCTGGTAGCCTTTCGCACTGGCTTCCACCTCAACGTAATCGCGCGTATCCATCAGAATGTGCGCAGTTTCACGGCTATTCACCACATAACCTGATCCGGTGGAGATCGCCAGTTGCGGCTGTTTGTTCTTCCGTGAGgagaga                    <  8:401358/148‑1 (MQ=255)
                                                                                                                                cAGCTGGTAGCCTTGCGCACGGGCTTCCACGTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGc                 <  3:7143/148‑1 (MQ=255)
                                                                                                                                cAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGc                 >  8:130871/1‑148 (MQ=255)
                                                                                                                                  gCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCAc               >  5:389472/1‑148 (MQ=255)
                                                                                                                                    tGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTgg                                                             >  4:396071/1‑100 (MQ=255)
                                                                                                                                    tGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTgg                                                             <  3:396071/100‑1 (MQ=255)
                                                                                                                                             tGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCACTGCTTCCAgtt    >  3:304910/1‑148 (MQ=255)
                                                                                                                                             tGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCACTGCATCCAgtt    >  7:197664/1‑148 (MQ=255)
                                                                                                                                               cgcACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCACTGCATCCAgttgt  <  8:244301/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TCAGACTGCCAGCGGTGCGCGGTTTCCCAACTCACCTGCTGGCCCTCAAAACCGAGCGTTTGCAACTGTTCATCGGCAATGATTTGATTGACGATAGTGGTAAGCGTGTTCGTCGCGTCAAGCAAATGCAGCTGGTAGCCTTGCGCACGGGCTTCCACCTCAACGTAATAGCGCGAATCCACCAGAATGTGCGCACTTTCACGGCTAATCACCACATAACCTGATCCGGTGGAGATCCCCAGATGCGGCTGTTTGTTCTGCCGTGAGGAGAGAAGCACTGCATCCAGTTGT  >  NC_000913/2506301‑2506591

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