Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,189,278 G→C 7.6% L458V (CTG→GTG)  PP_2799 ← class III aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,189,2780GC 7.6% 192.5 / 4.5 63L458V (CTG→GTG) PP_2799class III aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (30/28);  new base C (3/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

ACTTGGGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGT  >  NC_002947/3189131‑3189424
                                                                                                                                                   |                                                                                                                                                   
acTTGGGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAgg                                                                                                                                                    >  7:284329/1‑149 (MQ=255)
    gggCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTg                                                                                                                                                  <  3:63707/147‑1 (MQ=255)
        cATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGa                                                                                                                                            <  6:7596/149‑1 (MQ=255)
                  cccATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAACCCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCACCTGCtgt                                                                                                                                              >  6:12255/1‑137 (MQ=255)
                  cccATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCACCTGCtgt                                                                                                                                              <  5:12255/137‑1 (MQ=255)
                   ccATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCaa                                                                                                                                  >  5:123382/1‑148 (MQ=255)
                     aTGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATcc                                                                                                                               <  5:160608/149‑1 (MQ=255)
                     aTGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATcc                                                                                                                               >  4:197809/1‑149 (MQ=255)
                                  cGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTcagcca                                                                                                                  <  2:5237/149‑1 (MQ=255)
                                         ggAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCg                                                                                                           >  4:91101/1‑149 (MQ=255)
                                                               tcatcaCTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑c                                                                                     <  1:23533/148‑1 (MQ=255)
                                                                 atcaCTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTc                                                                                  >  6:85295/1‑149 (MQ=255)
                                                                    aCTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGaa                                                                                                     >  2:354660/1‑128 (MQ=255)
                                                                    aCTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGaa                                                                                                     <  1:354660/128‑1 (MQ=255)
                                                                           tCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTcc                                                                        >  3:313349/1‑149 (MQ=255)
                                                                             acgcacgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCGGAAGATCG‑CCTCGATATCGtgcgg                                                                      >  4:267166/1‑145 (MQ=255)
                                                                               gcacgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTggg                                                                    >  4:241513/1‑149 (MQ=255)
                                                                                 acgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGggtg                                                                  >  6:93935/1‑149 (MQ=255)
                                                                                  cgGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGggtgg                                                                 <  3:241513/149‑1 (MQ=255)
                                                                                     tCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATc                                                                            >  2:252140/1‑135 (MQ=255)
                                                                                     tCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATc                                                                            <  1:252140/135‑1 (MQ=255)
                                                                                       aaaCCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTgg                                                                     <  7:66355/140‑1 (MQ=255)
                                                                                       aaaCCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTgg                                                                     >  8:66355/1‑140 (MQ=255)
                                                                                         aCCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCaa                                                          <  5:85295/149‑1 (MQ=255)
                                                                                               aTGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACggag                                                     >  7:148523/1‑148 (MQ=255)
                                                                                               aTGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACggag                                                     <  3:267166/148‑1 (MQ=255)
                                                                                               aTGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACgga                                                      <  8:157423/147‑1 (MQ=255)
                                                                                               aTGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACgga                                                      <  3:197809/147‑1 (MQ=255)
                                                                                               aTGATTGTTGTTTTACGGTAGTAGCCATCAGATGCATCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACgga                                                      <  8:152672/147‑1 (MQ=255)
                                                                                                 gATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACggagga                                                   >  3:46464/1‑148 (MQ=255)
                                                                                                 gATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGATGCTCTGAGCGCGCCCAATACGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACggagga                                                   >  5:251643/1‑148 (MQ=255)
                                                                                                    tgttgtTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAg                                                >  2:266428/1‑148 (MQ=255)
                                                                                                     gttgttTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGACCTGGGTGGTAATCAACGGAGGACAg                                                >  7:254857/1‑147 (MQ=255)
                                                                                                      ttgttTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTAGCGTAAGCCAGCGTACGCTCGAAGGCGCTGAAGATCG‑CCTGGCTTTCGTCCAGGGTGGTAATCAACGGTGGACAGAAg                                             >  8:328000/1‑149 (MQ=255)
                                                                                                        gttTTACGTGAGTAGCCATGAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAgg                                            <  5:100850/148‑1 (MQ=255)
                                                                                                         ttttACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGc                                           >  8:182590/1‑148 (MQ=255)
                                                                                                             aCGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATc                                       <  7:182590/148‑1 (MQ=255)
                                                                                                               gTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCgtg                                    >  5:312630/1‑149 (MQ=255)
                                                                                                               gTGAGTAGCCATCAGATACAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCgtg                                    >  7:317227/1‑149 (MQ=255)
                                                                                                                tGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCGCCCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCgtg                                    <  1:266428/149‑1 (MQ=255)
                                                                                                                    tAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCg                                 >  6:296910/1‑147 (MQ=255)
                                                                                                                    tAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTa                               <  5:296910/149‑1 (MQ=255)
                                                                                                                        cATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACgg                                                       >  3:8147/1‑121 (MQ=255)
                                                                                                                         aTCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGggtg                                                                  >  4:340801/1‑109 (MQ=255)
                                                                                                                         aTCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGggtg                                                                  <  3:340801/109‑1 (MQ=255)
                                                                                                                         aTCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCgtgt                                   >  8:211580/1‑140 (MQ=255)
                                                                                                                         aTCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCgtgt                                   <  7:211580/140‑1 (MQ=255)
                                                                                                                             gATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGgcgc                       >  6:93891/1‑148 (MQ=255)
                                                                                                                               tGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGaa                    <  6:123382/149‑1 (MQ=255)
                                                                                                                                ccAGCCTTTGCTCACAGCAGGTGCTGTGTGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGaa                    <  7:337800/147‑1 (MQ=255)
                                                                                                                                 cAGCCTTTGCTCACAGCAGTTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAAt                   <  8:284329/148‑1 (MQ=255)
                                                                                                                                   gCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGa                 >  4:173597/1‑148 (MQ=255)
                                                                                                                                   gCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGACCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGAt                <  8:254857/149‑1 (MQ=255)
                                                                                                                                     cTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCa              <  8:317227/149‑1 (MQ=255)
                                                                                                                                     cTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCa              <  3:173597/149‑1 (MQ=255)
                                                                                                                                          cTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGccc          <  8:148523/148‑1 (MQ=255)
                                                                                                                                          cTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCa         >  6:104828/1‑149 (MQ=255)
                                                                                                                                           tCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCAt        >  8:203912/1‑149 (MQ=255)
                                                                                                                                           tCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAAcctt         >  5:347803/1‑146 (MQ=255)
                                                                                                                                           tCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTAGCGTCAGCCAGAGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCAt        >  6:87131/1‑149 (MQ=255)
                                                                                                                                           tCACAGCACCTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAAc                                                         <  7:98822/100‑1 (MQ=255)
                                                                                                                                           tCACAGCACCTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAAc                                                         >  8:98822/1‑100 (MQ=255)
                                                                                                                                           tCACAGCACCTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCAt        >  2:337429/1‑149 (MQ=255)
                                                                                                                                             acaGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATg       <  4:46464/148‑1 (MQ=255)
                                                                                                                                              cagcagGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATc               >  8:100360/1‑139 (MQ=255)
                                                                                                                                              cagcagGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATc               <  7:100360/139‑1 (MQ=255)
                                                                                                                                               agcagGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCtt    >  1:268793/1‑149 (MQ=255)
                                                                                                                                                 cagGTTCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTtgt  <  7:328000/149‑1 (MQ=255)
                                                                                                                                                 cagGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTtgt  <  5:93891/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACTTGGGCCATGATGGTTCCCATGGCACCGAGCCCGATCACGGAAACGTCGAACTCGAAGCTTTCATCACTGTGTTCACGCACGGTCAAACCCTCATGATTGTTGTTTTACGTGAGTAGCCATCAGATGCAGCCTTTGCTCACAGCAGGTGCTGTGAGCGCGCCCAATCCGTCGCGTCAGCCAGCGTCCGCTCGAAGGCGCTGAAGATCG‑CCTCGATATCGTCCTGGGTGGTAATCAACGGAGGACAGAAGGCGATCGTGTCGTAGATGGCGCGAATGATCAGCCCATGCTTGT  >  NC_002947/3189131‑3189424

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