Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,051,898 C→T 6.7% R186H (CGT→CAT)  PP_3571 ← acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,051,8980CT 6.7% 229.0 / 3.9 77R186H (CGT→CAT) PP_3571acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base C (35/36);  new base T (2/3);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

TGGATCGAAGCGCTCGTCGTGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAA  >  NC_002947/4051750‑4052044
                                                                                                                                                    |                                                                                                                                                  
tGGATCGAAGCGCTCGTCGTGACGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAAc                                                                                                                                                    <  6:74131/149‑1 (MQ=255)
      gAAGCGCTCGTCGTGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCAtt                                                                                                                                               <  1:278841/148‑1 (MQ=255)
                  gtGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCTTGCACCTGGCAATGCGTTGCCAggtgcagg                                                                                                                                   <  1:215555/148‑8 (MQ=255)
                    gatgaACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTtg                                                                                                                                >  3:113288/1‑149 (MQ=255)
                    gatgaACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTtg                                                                                                                                >  6:126806/1‑149 (MQ=255)
                       gaACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGc                                                                                                                             >  1:350363/1‑149 (MQ=255)
                       gaACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGc                                                                                                                             >  7:11251/1‑149 (MQ=255)
                        aaCCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCAtt                                                                                                                                               >  5:298700/1‑130 (MQ=255)
                        aaCCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCAtt                                                                                                                                               <  6:298700/130‑1 (MQ=255)
                        aaCCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGc                                                                                                                             <  4:72727/148‑1 (MQ=255)
                         aCCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGccc                                                                                                                           <  3:260818/149‑1 (MQ=255)
                          ccGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAATGCGTTGCCAgg                                                                                                                                         >  8:81615/1‑133 (MQ=255)
                          ccGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAATGCGTTGCCAgg                                                                                                                                         <  7:81615/134‑2 (MQ=255)
                          ccGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGccc                                                                                                                           >  4:304467/1‑148 (MQ=255)
                           cGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGccc                                                                                                                           >  4:296458/1‑147 (MQ=255)
                             gCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAg                                                                                                                        >  4:84398/1‑148 (MQ=255)
                             gCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGACCTGCACCTGGCAACGCATTGCCAGCTAGCCCTTGTGCACCAg                                                                                                                        >  8:313598/1‑148 (MQ=255)
                                     agggcCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGttt                                                                                                               <  2:350363‑M2/144‑1 (MQ=255)
                                     agggcCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGttt                                                                                                               <  5:247160‑M2/144‑1 (MQ=255)
                                     agggcCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTAGTGCCCCAGTACCGGttt                                                                                                               <  5:227561‑M2/144‑1 (MQ=255)
                                           tCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCACCGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTc                                                                                                         >  7:286560/1‑149 (MQ=255)
                                                gATTTCACCCAACTTGCCAATCAGTTTCGCCGCATTCTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTgg                                                                                                    <  3:48940/149‑1 (MQ=255)
                                                      aCCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCACTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCAc                                                                                              >  4:146614/1‑149 (MQ=255)
                                                           aCTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATc                                                                                         <  1:74412/149‑1 (MQ=255)
                                                               gCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAATGCGTTGCCAggtg                                                                                                                                       >  5:19500/1‑96 (MQ=255)
                                                               gCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAATGCGTTGCCAggtg                                                                                                                                       <  6:19500/99‑4 (MQ=255)
                                                               gCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCa                                                                                        >  3:180135/1‑146 (MQ=255)
                                                               gCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCa                                                                                        <  4:180135/146‑1 (MQ=255)
                                                               gCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAgg                                                                                      <  4:113288/148‑1 (MQ=255)
                                                               gCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAgg                                                                                      <  3:84398/148‑1 (MQ=255)
                                                               gCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAg                                                                                       <  8:123663/147‑1 (MQ=255)
                                                                  aaTCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCAcc                                                                                             >  7:142595/1‑138 (MQ=255)
                                                                  aaTCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCAcc                                                                                             <  8:142595/138‑1 (MQ=255)
                                                                  aaTCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGca                                                                                    >  6:304915/1‑147 (MQ=255)
                                                                   aTCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTg                                                                                                     >  7:111655/1‑129 (MQ=255)
                                                                   aTCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTg                                                                                                     <  8:111655/129‑1 (MQ=255)
                                                                            gccgcATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGtt                                                                         <  7:313598/148‑1 (MQ=255)
                                                                                cATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTacac                                                                    >  2:352129/1‑149 (MQ=255)
                                                                                  tACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAAT‑GCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACggg                                                                 >  8:235851/1‑149 (MQ=255)
                                                                                     tCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTac                                                                      <  4:128591/142‑1 (MQ=255)
                                                                                     tCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTac                                                                      >  3:128591/1‑142 (MQ=255)
                                                                                     tCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGgcg                                                               >  1:109728/1‑149 (MQ=255)
                                                                                       gATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGAt                                                                                          <  1:91295/120‑1 (MQ=255)
                                                                                       gATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGAt                                                                                          >  2:91295/1‑120 (MQ=255)
                                                                                         tGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGt                                                           <  3:146614/149‑1 (MQ=255)
                                                                                           gCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGtt                                                          <  2:103218/148‑1 (MQ=255)
                                                                                             gTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGa                                                       <  8:11251/149‑1 (MQ=255)
                                                                                              ttCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAg                                                      <  5:126806/149‑1 (MQ=255)
                                                                                                 aCCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAggg                                                    >  4:108310/1‑148 (MQ=255)
                                                                                                 aCCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAggg                                                    >  2:245179/1‑148 (MQ=255)
                                                                                                 aCCCCATACGGCGCGTACGCCGAAT‑GCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGc                                                   <  7:235851/148‑1 (MQ=255)
                                                                                                  ccccATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGca                                                  <  2:109728/149‑1 (MQ=255)
                                                                                                   cccATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGt                                                                                                                 <  7:252586/85‑1 (MQ=255)
                                                                                                   cccATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGt                                                                                                                 >  8:252586/1‑85 (MQ=255)
                                                                                                   cccATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcag                                                 >  5:92904/1‑149 (MQ=255)
                                                                                                    ccATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGa                                                       <  4:205376/142‑1 (MQ=255)
                                                                                                    ccATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTGTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGa                                                       >  3:205376/1‑142 (MQ=255)
                                                                                                           gcgcGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATc                                         <  5:212984/149‑1 (MQ=255)
                                                                                                           gcgcGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAAAGCATTGCCAGCTTGCCATTGTGCCCCAGTACCGGTTTGAGTTAGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGTGCGCTGTTGGAGGGCAGACAGCATc                                         >  1:80409/1‑149 (MQ=255)
                                                                                                                   gCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCAc                                 >  5:131852/1‑149 (MQ=255)
                                                                                                                     cGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGc                               >  6:333223/1‑149 (MQ=255)
                                                                                                                     cGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGc                               <  5:333223/149‑1 (MQ=255)
                                                                                                                      gAATGGCACGCTGCGCCCTGCACCTGTCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcagcca                                              <  6:218083/133‑1 (MQ=255)
                                                                                                                      gAATGGCACGCTGCGCCCTGCACCTGCCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcagcca                                              >  5:218083/1‑133 (MQ=255)
                                                                                                                         tGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCa                            >  6:154215/1‑148 (MQ=255)
                                                                                                                           gCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCaaa                          <  6:92904/148‑1 (MQ=255)
                                                                                                                            cACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTCCACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAaac                         >  3:226326/1‑148 (MQ=255)
                                                                                                                                   cgcCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGa                   <  1:352129/147‑1 (MQ=255)
                                                                                                                                      ccTGCACCTGGCGACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGActtc               >  8:346759/1‑148 (MQ=255)
                                                                                                                                      ccTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGActtct              >  1:126975/1‑149 (MQ=255)
                                                                                                                                          cACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACttcttc            <  5:304915/147‑1 (MQ=255)
                                                                                                                                            ccTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCa        >  5:327555/1‑149 (MQ=255)
                                                                                                                                             cTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCa        >  4:85990/1‑148 (MQ=255)
                                                                                                                                             cTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCAt       <  6:327555/149‑1 (MQ=255)
                                                                                                                                                ggaaCGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACAGGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATa      <  2:116146/145‑1 (MQ=255)
                                                                                                                                                gcgaCGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATa      <  7:346759/144‑1 (MQ=255)
                                                                                                                                                gCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATa      <  3:304467/147‑1 (MQ=255)
                                                                                                                                                   aCGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGaa  >  8:283346/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGGATCGAAGCGCTCGTCGTGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAA  >  NC_002947/4051750‑4052044

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