Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 473,813 C→A 14.7% A283E (GCG→GAG)  amtB → ammonium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913473,8130CA14.7% 188.4 / 17.7 67A283E (GCG→GAG) amtBammonium transporter
Reads supporting (aligned to +/- strand):  ref base C (28/29);  new base A (5/5);  total (33/34)
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.95e-01

GTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCTTGCTGCGGGTGGATGA  >  NC_000913/473667‑473960
                                                                                                                                                  |                                                                                                                                                   
gttTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCgt                                                                                                                                                   <  1:200158/149‑1 (MQ=255)
    ggCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGAGTGtt                                                                                                                                                <  1:181250/148‑1 (MQ=255)
      cTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTg                                                                                                                                             >  3:10273/1‑149 (MQ=255)
       tttAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTgg                                                                                                                                            >  3:209369/1‑149 (MQ=255)
       tttAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTgg                                                                                                                                            >  3:340360/1‑149 (MQ=255)
             gCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGAt                                                                                                                                      <  3:296021/149‑1 (MQ=255)
                gggTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCCCTGCTGGGGGCGTGTTCTGGCGCGATTGc                                                                                                                                   >  4:113700/1‑149 (MQ=255)
                gggTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGc                                                                                                                                   >  4:380160/1‑149 (MQ=255)
                  gTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCg                                                                                                                                 <  4:191635/149‑1 (MQ=255)
                    cAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGt                                                                                                                               <  4:209369/149‑1 (MQ=255)
                     aGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTc                                                                                                                              >  3:301758/1‑149 (MQ=255)
                       cGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTg                                                                                                                            <  4:340360/149‑1 (MQ=255)
                         ggCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGt                                                                                                                          <  8:151546/149‑1 (MQ=255)
                          gCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTc                                                                                                                         >  7:7673/1‑149 (MQ=255)
                              ggCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCg                                                                                                                     >  7:123056/1‑149 (MQ=255)
                               gCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGt                                                                                                                    <  6:163973/149‑1 (MQ=255)
                                 gAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGa                                                                                                                  >  8:77089/1‑149 (MQ=255)
                                       aaTCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCa                                                                                                             <  4:10273/148‑1 (MQ=255)
                                               cACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGAtt                                                                                                                                     >  1:134844/1‑116 (MQ=255)
                                               cACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGAtt                                                                                                                                     <  2:134844/116‑1 (MQ=255)
                                                    gCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACAt                                                                                               >  8:243207/1‑149 (MQ=255)
                                                      aTTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTg                                                                                             <  4:154011/149‑1 (MQ=255)
                                                       tttGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTgg                                                                                            >  6:168111/1‑149 (MQ=255)
                                                       tttGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATATTCGGTGAATGGGCGCTGCGTGGTAAGCCATCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTTACTCCAGCCTGCGGCTCCATTgg                                                                                            >  5:12184/1‑149 (MQ=255)
                                                         tgtgAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTggg                                                                                           <  6:56534/148‑1 (MQ=255)
                                                          gtgAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGAGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGt                                                                                         >  3:356981/1‑149 (MQ=255)
                                                             aaTACTTGGGGCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGGCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTgg                                                                                      <  3:203551/149‑1 (MQ=255)
                                                             aaTACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCTCTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTgg                                                                                      >  1:41040/1‑149 (MQ=255)
                                                             aaTACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGAGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTgg                                                                                      >  1:258017/1‑149 (MQ=255)
                                                               tACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTggcg                                                                                    >  6:153337/1‑149 (MQ=255)
                                                                 cTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGTGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTggcggc                                                                                  <  8:7673/149‑1 (MQ=255)
                                                                  tgtgGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGgcggcg                                                                                 >  8:174851/1‑149 (MQ=255)
                                                                    tgGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCGGCGGGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGgcggcg                                                                                 <  3:54229/147‑1 (MQ=255)
                                                                      gTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCATGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGcgcg                                                                               >  6:57425/1‑147 (MQ=255)
                                                                      gTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCATGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGcgcg                                                                               <  5:57425/147‑1 (MQ=255)
                                                                         gCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGGCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGa                                                                           <  2:109445/148‑1 (MQ=255)
                                                                         gCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGAGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCg                                                                                                      >  6:42647/1‑121 (MQ=255)
                                                                            aCGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTAt                                                                       <  8:123056/149‑1 (MQ=255)
                                                                              ggcggcggcAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCg                                                                     <  4:301758/149‑1 (MQ=255)
                                                                                       aaTTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGTAGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAg                                                            >  5:129256/1‑149 (MQ=255)
                                                                                            ttGGCTGGAACTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGt                                                       <  7:206359/149‑1 (MQ=255)
                                                                                               gCTGGATCTTCGGTGCATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGGGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCt                                                     <  7:77089/148‑1 (MQ=255)
                                                                                                   gATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGAGTGTTCTGGCGCGATTGCCGGTCTGGt                                                                                                                          >  8:121612/1‑75 (MQ=255)
                                                                                                   gATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGAGTGTTCTGGCGCGATTGCCGGTCTGGt                                                                                                                          <  7:121612/75‑1 (MQ=255)
                                                                                                       ttCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCtt                                            >  2:218285/1‑149 (MQ=255)
                                                                                                        tCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTtg                                           >  4:340938/1‑149 (MQ=255)
                                                                                                          ggTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTtgtg                                         <  5:54038/149‑1 (MQ=255)
                                                                                                             gAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTgggg                                      >  1:222736/1‑149 (MQ=255)
                                                                                                              actGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGAGTGTTCTGGCGCGATTGCCGGTATGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGc                                     <  3:294393/147‑1 (MQ=255)
                                                                                                                   gcgcTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTAc                                <  5:153337/149‑1 (MQ=255)
                                                                                                                       tGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATg                            >  1:195337/1‑149 (MQ=255)
                                                                                                                       tGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTATGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATg                            >  7:170540/1‑149 (MQ=255)
                                                                                                                             gTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTgggg                                                                                          <  3:172106/81‑1 (MQ=255)
                                                                                                                             gTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTgggg                                                                                          >  4:172106/1‑81 (MQ=255)
                                                                                                                             gTAAGCCTTCACTGCTGGGGGCGTGTGCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCaaa                      <  1:372190/149‑1 (MQ=255)
                                                                                                                              tAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAAc                     <  2:41040/149‑1 (MQ=255)
                                                                                                                              tAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAAc                     >  8:278886/1‑149 (MQ=255)
                                                                                                                                aGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGc                   >  2:282217/1‑149 (MQ=255)
                                                                                                                                aGCCTTCACTGCTGGGGGAGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACg                    <  4:356981/148‑1 (MQ=255)
                                                                                                                                  cattCACTGCGGGGGGAGTGTTCTGGCGCGATTGACGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCtt                 <  6:129256/147‑1 (MQ=255)
                                                                                                                                    ttCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCTtgc               >  3:126872/1‑149 (MQ=255)
                                                                                                                                       aCTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCTtgctgc            >  7:102726/1‑149 (MQ=255)
                                                                                                                                           ctgGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCTTGCTGCggg         <  1:282217/148‑1 (MQ=255)
                                                                                                                                           ctgGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCTTGCTGCggg         <  7:243207/148‑1 (MQ=255)
                                                                                                                                           ctgGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGGGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCTTGCTGCgg          <  3:380160/147‑1 (MQ=255)
                                                                                                                                              ggggCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCTTGCTGCGGGTGga     >  2:399358/1‑149 (MQ=255)
                                                                                                                                                 gCGTGTGATGGCTGGATTGCCTGTAGGGTCGGCGGGTCGCCAGCCTGTGGTTAGATTGGGGGTGGAGGCGCGTTTATTATAGGCGTGGTATCTGGGCTGGCGGGCTTGTGGGGCGGTACCATGCTGAAACGCTTGCTGCGGGTGgaaga  <  6:12184/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTTTGGCTTTAACGCCGGGTCAGCGGGCACGGCGAATGAAATCGCGGCACTGGCATTTGTGAATACTGTGGTCGCAACGGCGGCGGCAATTCTTGGCTGGATCTTCGGTGAATGGGCGCTGCGTGGTAAGCCTTCACTGCTGGGGGCGTGTTCTGGCGCGATTGCCGGTCTGGTCGGCGTGACGCCAGCCTGCGGCTACATTGGGGTTGGCGGCGCGTTGATTATCGGCGTGGTAGCTGGTCTGGCGGGCTTGTGGGGCGTTACCATGCTCAAACGCTTGCTGCGGGTGGATGA  >  NC_000913/473667‑473960

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: