Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,051,901 A→G 6.5% M185T (ATG→ACG)  PP_3571 ← acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,051,9010AG 6.5% 223.9 / 4.7 77M185T (ATG→ACG) PP_3571acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base A (36/36);  new base G (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 = 9.15e-01

GAAGCGCTCGTCGTGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGC  >  NC_002947/4051756‑4052049
                                                                                                                                                 |                                                                                                                                                    
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                                                                                                                                            <  6:19500/99‑4 (MQ=255)
                                                         gCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAATGCGTTGCCAggtg                                                                                                                                            >  5:19500/1‑96 (MQ=255)
                                                         gCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCa                                                                                             <  4:180135/146‑1 (MQ=255)
                                                         gCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCa                                                                                             >  3:180135/1‑146 (MQ=255)
                                                         gCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAgg                                                                                           <  3:84398/148‑1 (MQ=255)
                                                         gCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAgg                                                                                           <  4:113288/148‑1 (MQ=255)
                                                         gCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAg                                                                                            <  8:123663/147‑1 (MQ=255)
                                                            aaTCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCAcc                                                                                                  <  8:142595/138‑1 (MQ=255)
                                                            aaTCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCAcc                                                                                                  >  7:142595/1‑138 (MQ=255)
                                                            aaTCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGca                                                                                         >  6:304915/1‑147 (MQ=255)
                                                             aTCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTg                                                                                                          <  8:111655/129‑1 (MQ=255)
                                                             aTCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTg                                                                                                          >  7:111655/1‑129 (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)
                                                                                                                                                 aTTGCCAGCTTGCCCTTGTGCCCCAGTACCTGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGATGGCATCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGc  <  2:80409/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GAAGCGCTCGTCGTGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAGGAAAACGC  >  NC_002947/4051756‑4052049

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