Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 235,330 T→C 8.2% intergenic (+27/+98) PP_0183 → / ← argH glutathione S‑transferase/argininosuccinate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947235,3300TC 8.2% 237.7 / 9.7 85intergenic (+27/+98)PP_0183/argHglutathione S‑transferase/argininosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base T (39/39);  new base C (5/2);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

CCGCAGGCCAAGGCGTTGCTGCAACTGATGGGGGAGAACCCGCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGATTAGCGAGACGCCAACAGGGCCTTGCCACGAACGACGGCAGCACGCACCT  >  NC_002947/235184‑235477
                                                                                                                                                   |                                                                                                                                                    
ccGCAGGCCAAGGCGTTGCTGCAACTGATGGGGGAGAACCCGCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGc                                                                                                                                                    >  8:6318/1‑149 (MQ=255)
  gCAGGCCAAGGCGTTGCTGCAACTGATGGGGGAGAACCCGCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGccc                                                                                                                                                  >  7:24836/1‑149 (MQ=255)
       ccAAGGCGTTGCTGCAACTGATGGGGGAGAACCCGCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGcgaa                                                                                                                                             <  2:262752/149‑4 (MQ=255)
       ccAAGGCGTTGCTGCAACTGATGGGGGAGAACCCGCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTcc                                                                                                                                             <  6:26477/149‑1 (MQ=255)
           ggCGTTGCTGCAACTGATGGGGGAGAACCCGCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCACGGGCATGCCCGCTCCca                                                                                                                                           >  3:68758/1‑147 (MQ=255)
           gaaggctcggtgaatgcgcgtgaccacattggtggtacggcgccggcgcaggtgcgtgctgccgtcgttcgtggcaaggccctgttggcgtctcgctaatcccccaaggcgagggcttcgccctcgTTCGCGGGCGtgcccgcgaacga                                                                                                                                          >  5:8651‑M1/127‑143 (MQ=255)
           gaaggctcggtgaatgcgcgtgaccacattggtggtacggcgccggcgcaggtgcgtgctgccgtcgttcgtggcaaggccctgttggcgtctcgctaatcccccaaggcgagggcttcgccctcgTTCGCGGGCGtgcccgcgaacga                                                                                                                                          >  2:140152‑M1/127‑143 (MQ=255)
              gTTGCTGCAACTGATGGGGGAGAACCCGCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGt                                                                                                                                      <  7:261081/149‑1 (MQ=255)
                   tgcAACTGATGGGGGAGAACCCGCACAGTGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTg                                                                                                                                  >  1:206582/1‑148 (MQ=255)
                           aTGGGGGAGAACCCGCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGttttt                                                                                                                        <  7:6318/149‑1 (MQ=255)
                              ggggAGAACCCGCACATGCCGCGGATCTTTGCTGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGa                                                                                                                      <  8:76042/148‑1 (MQ=255)
                               gggAGAACCCGCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGaaa                                                                                                                    >  5:88620/1‑149 (MQ=255)
                                ggAGAACCCGCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGGTGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAAc                                                                                                                   >  1:81630/1‑149 (MQ=255)
                                  agaACCCGCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCACGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACtt                                                                                                                 <  4:68758/149‑1 (MQ=255)
                                   cacattggtggtacggcgccgtcgcaggtgcgtgctgccgtcgttcgtgtcaaggcccttttggcttctcgctaatcccccaaggcgtgggcttcgccctcgTTCTCTGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACtt                                                                                                                 >  7:272115‑M1/103‑149 (MQ=255)
                                   cacattggtggtacggcgccggcgcaggtgcgtgctgccgtcgttcgtggcaaggccctgttggcgtctcgctaatcccccaaggcgagggcttcgccctcgTTCGCGGGCACGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACtt                                                                                                                 >  3:142945‑M1/103‑149 (MQ=255)
                                    acattggtggtacggcgccggcgcaggtgcgtgctgccgtcgttcgtggcaaggccctgttggcgtctcgctaatcccccaaggcgagggcttcgccctcgTTCGCGGGCACGCCCGCTCCCACAGGTAATGCACCGTTTTTGAAACtt                                                                                                                 >  2:90299‑M1/102‑149 (MQ=255)
                                      cccGCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGc                                                                                                             <  8:24836/149‑1 (MQ=255)
                                       ccGCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCa                                                                                                            >  4:116720/1‑149 (MQ=255)
                                       ccGCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCa                                                                                                            >  2:254687/1‑149 (MQ=255)
                                         gCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCAtt                                                                                                          <  5:301600/149‑1 (MQ=255)
                                         gCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCAt                                                                                                           >  4:317563/1‑148 (MQ=255)
                                         gCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCAt                                                                                                           >  3:18348/1‑148 (MQ=255)
                                         gCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCAt                                                                                                           >  3:13267/1‑148 (MQ=255)
                                               acggcgccggcgcaggtgcgtgctgccgtcgttcgtggcaaggccctgttggcgtctcgctaatcccccaaggcgagggcttcgccctcgTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCtg                                                                                                     <  4:227016‑M1/59‑1 (MQ=255)
                                                 ggcgccggcgcaggtgcgtgctgccgtcgttcgtggcaaggccctgttggcgtctcgctaatcccccaaggcgagggcttcgccctcgTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCtgtg                                                                                                   >  8:110684‑M1/89‑149 (MQ=255)
                                                 ggcgccggcgcaggtgcgtgctgccgtcgttcgtggcaaggccctgttggcgtctcgctaatcccccaaggcgagggcttcgccctcgTTCGCGGGCACGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACtt                                                                                                                 >  3:245354‑M1/89‑135 (MQ=255)
                                                 ggcgccggcgcaggtgcgtgctgccgtcgttcgtggcaaggccctgttggcgtctcgctaatcccccaaggcgagggcttcgccctcgTTCGCGGGCACGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACtt                                                                                                                 <  4:245354‑M1/47‑1 (MQ=255)
                                                  gcgcgggcgcaggtgcgtgctgccgtcgttcgtggcaaggccctgttggcgtctcgctaatcgcccactgcgagggcttcgccctcgTTCGCGGGCACGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTgg                                                                                                  >  3:365364‑M1/88‑149 (MQ=255)
                                                  gcgccggcgcaggtgcgtgccgccgtcgttcgtggcaaggccctgttggcgtctcgctaatcccccacggcgagggcttcgccctcgTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTgg                                                                                                  >  2:280115‑M1/88‑149 (MQ=255)
                                                   gcgGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGcgaa                                                                                                                                             <  3:162871/106‑4 (MQ=255)
                                                   gcgGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGcgaa                                                                                                                                             >  4:162871/1‑103 (MQ=255)
                                                    cgGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTtgt                                                                                                               >  4:153143/1‑134 (MQ=255)
                                                    cgGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTtgt                                                                                                               <  3:153143/134‑1 (MQ=255)
                                                              cGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGc                                                                                      <  3:48571/149‑1 (MQ=255)
                                                               ggACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGcc                                                                                     >  3:47024/1‑149 (MQ=255)
                                                                gACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGccc                                                                                    >  7:106406/1‑149 (MQ=255)
                                                                        ccgtcgttcgtggcaaggccctgttggcgtctcgctaatcccccaaggcgagggcttcgccctcgTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAAc                                                                              >  3:227016‑M1/66‑147 (MQ=255)
                                                                        aGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAAc                                                                              >  8:73849/1‑147 (MQ=255)
                                                                         cgtcgttcgtggcaaggccctgttggcgtctcgctaatcccccaaggcgagggcttcgccctcgTTCGCGGGCACGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTgg                                                                                                  <  7:262076‑M1/62‑1 (MQ=255)
                                                                         cgtcgttcgtggcaaggccctgttggcgtctcgctaatcccccaaggcgagggcttcgccctcgTTCGCGGGCACGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTgg                                                                                                  >  8:262076‑M1/65‑126 (MQ=255)
                                                                          catccaTGCCGCCCTTCAGGCATAGCACCCGCAGCGCCACCCGCTGTGTCTGGTGGTCACCCCTCCCCGCCCATCCCCGCTCCCCCAGGTCTGCCCCGTTTTTTAAACTTGTCCTTTACCTGTGGTAGCGGGCGTGCCCGCGAACGAgg                                                                          <  8:137412/144‑1 (MQ=255)
                                                                          cGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAgg                                                                          <  6:88620/149‑1 (MQ=255)
                                                                           gTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCCCCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAggg                                                                         <  5:357555/149‑1 (MQ=255)
                                                                           gTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCAGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAggg                                                                         >  2:274013/1‑149 (MQ=255)
                                                                             cGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGct                                                                       >  2:299615/1‑148 (MQ=255)
                                                                               aTGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCAGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGATACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGaa                                                                     <  1:274013/149‑1 (MQ=255)
                                                                                tGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAg                                                                    >  2:266379/1‑149 (MQ=255)
                                                                                 gCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGc                                                                   >  7:336153/1‑149 (MQ=255)
                                                                                        ttCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGcc                                                            >  5:278279/1‑149 (MQ=255)
                                                                                          cATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCCCCGTTTTTGAAACTTGTGCATTACCGGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGTAGCCCTCGCCtt                                                          >  1:156354/1‑149 (MQ=255)
                                                                                          cATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCTTCGCCCTCgttc                                                           <  1:299615/148‑4 (MQ=255)
                                                                                           aTGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTTGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTg                                                         >  5:233101/1‑149 (MQ=255)
                                                                                            tGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTgg                                                                                                  <  8:16103/107‑1 (MQ=255)
                                                                                            tGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTgg                                                                                                  >  7:16103/1‑107 (MQ=255)
                                                                                            tGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCTTCGCCCTCgttcgcg                                                        >  3:295376/1‑143 (MQ=38)
                                                                                            tGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTgg                                                        <  2:206582/149‑1 (MQ=255)
                                                                                             ggAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTTGATTAACTGTGGGAGCGGGCGTGCTCTCGAACGAGGGCGAAGCCCTCGGCGTggg                                                       >  4:5177/1‑149 (MQ=39)
                                                                                             ggAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCTTCGCCCTCgttcgcg                                                        <  4:295376/148‑7 (MQ=38)
                                                                                              gagaTGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTggg                                                       <  1:246434/148‑1 (MQ=255)
                                                                                              gagaTGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTggg                                                       <  3:317563/148‑1 (MQ=255)
                                                                                                 aggatCCGCAGCGGCATGCGCTGAGGCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGa                                                    <  4:80519/146‑1 (MQ=255)
                                                                                                  tgatCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGAtt                                                  <  5:175728/149‑1 (MQ=255)
                                                                                                   gatCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGAt                                                   >  1:313352/1‑147 (MQ=255)
                                                                                                   gatCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGAt                                                   >  8:311157/1‑147 (MQ=255)
                                                                                                     tCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGATTAGc                                               <  2:81630/149‑1 (MQ=255)
                                                                                                     gccGCTGCGGCACGCGCGGAGTCGGGGGTTGCCCCCTGCGCGGGCATGCCCGCCCCCCCAGGGCTGCACCGTTTTTCAAACTTGTGCATTACCTGGGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGATTAGc                                               <  5:43427/148‑1 (MQ=255)
                                                                                                      ccGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGc                                                                                                             <  6:20968/86‑1 (MQ=255)
                                                                                                      ccGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGc                                                                                                             >  5:20968/1‑86 (MQ=255)
                                                                                                      ccGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGATTAGCg                                              >  2:260477/1‑149 (MQ=255)
                                                                                                       cGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGATTAGc                                               >  5:377866/1‑147 (MQ=255)
                                                                                                          aGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGATTAGCGAGAc                                          <  3:77022/149‑1 (MQ=255)
                                                                                                          aGCGGCAAGCGCTGAGGCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGGGGGAGCGGGCGTGCTCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGATTAGCGAGAc                                          <  3:5177/149‑1 (MQ=255)
                                                                                                                aaGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTgg                                                                                                  >  7:209101/1‑87 (MQ=255)
                                                                                                                aaGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTgg                                                                                                  <  8:209101/87‑1 (MQ=255)
                                                                                                                 aGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGATTAGCGAGACGCCAACa                                   <  7:73849/149‑1 (MQ=255)
                                                                                                                    gcTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGATTAGCGAGACGCCAACAggg                                <  1:254687/149‑1 (MQ=255)
                                                                                                                            tGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGGCCCGTTTTTGAAACTTGTGGATTACCGGTGGGAGCGGGCGGGGCCGCGCACGCGGGCGAAGCCCGCGCCGTGGGGGATTAGCGAGACGCCAACAGGGCCTTGccc                        >  6:252008/1‑148 (MQ=255)
                                                                                                                            tGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAgg                                                                          >  3:247732/1‑99 (MQ=255)
                                                                                                                            tGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAgg                                                                          <  4:247732/99‑1 (MQ=255)
                                                                                                                            tGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGATTAGCGAGACGCCAACAGGGCCTTGCCa                        >  2:58449/1‑149 (MQ=255)
                                                                                                                                    aGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGATTAGCGAGACGCCAACAGGGCCTTGCCACGAacgac                <  3:116720/149‑1 (MQ=255)
                                                                                                                                    aGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGATTAGCGAGACGCCAACAGGGCCTTGCCACGAacgac                <  7:311157/149‑1 (MQ=255)
                                                                                                                                     gCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGATTAGCGAGACGCCAACAGGGCCTTGCCACGAacgacg               <  8:106406/149‑1 (MQ=255)
                                                                                                                                      cccTTCGGGGGCATGCCCGAGCCCACAGGTGTCCACCTTTTTTTAACCTTGTGCATTACCTGTGGTTGGGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGATTAGCGAGACGCCAACAGGGCCTTGCCACGAACGACgg              <  6:233101/149‑1 (MQ=255)
                                                                                                                                       ccTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGGGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGATTAGCGAGACGCCAACAGGGCCTTGCCACGAACGACGgc             <  2:313352/149‑1 (MQ=255)
                                                                                                                                                  ATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGATTAGCGAGACGCCAATAGGGCCTTGCCACGAACGACGGCAGCACGCACCt  <  2:230497‑M1/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCGCAGGCCAAGGCGTTGCTGCAACTGATGGGGGAGAACCCGCACATGCCGCGGATCTTGGCGGACAAGGAAGCGTCGATGCCGGCCTTCATGGAGATGATCCGCAGCGGCAAGCGCTGAGTCTGGTGTTGAGCCCTTCGCGGGCATGCCCGCTCCCACAGGTCTGCACCGTTTTTGAAACTTGTGCATTACCTGTGGGAGCGGGCGTGCCCGCGAACGAGGGCGAAGCCCTCGCCTTGGGGGATTAGCGAGACGCCAACAGGGCCTTGCCACGAACGACGGCAGCACGCACCT  >  NC_002947/235184‑235477

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