Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,051,898 C→T 7.8% 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 7.8% 263.1 / 8.5 89R186H (CGT→CAT) PP_3571acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base C (37/45);  new base T (3/4);  total (40/49)
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.70e-01

CGAAGCGCTCGTCGTGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAG  >  NC_002947/4051755‑4052041
                                                                                                                                               |                                                                                                                                               
cGAAGCGCTCGTCGTGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCAtt                                                                                                                                            >  3:126019/1‑149 (MQ=255)
cGAAGCGCTCGTCGTGATGAACCGGCAGTGCCAGGGCCTCGCCGAATTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCAtt                                                                                                                                            >  1:288545/1‑149 (MQ=255)
 gAAGCGCTCTTCGTGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTTCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCTCCCTGCACCTGGCAACGCATTg                                                                                                                                           <  5:126850/149‑1 (MQ=255)
 gAAGCGCTCGTCGTGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCAtt                                                                                                                                            <  5:7378/148‑1 (MQ=255)
       cTCGTCGTGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCt                                                                                                                                     >  4:261083/1‑149 (MQ=255)
        tcgtcgTGATGTACCGGCAGTGCAAGGGCGTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCGCATACGGCGCGTACGCCGAATGGCACTCTGCGCCCTGCACCTGGCAACGCATTGCCAGCtt                                                                                                                                    <  2:369110/149‑1 (MQ=255)
          gtcgtGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCtt                                                                                                                                    <  3:30422/147‑1 (MQ=255)
          gtcgtGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCtt                                                                                                                                    <  3:147076/147‑1 (MQ=255)
          gtcgtGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGc                                                                                                                                  >  8:138429/1‑149 (MQ=255)
             gtGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGccc                                                                                                                                <  3:253656/148‑1 (MQ=255)
              tgatgaACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCtt                                                                                                                              <  3:316810/149‑1 (MQ=255)
                atgaACCGGCAGTGCCAGGGCCTCGCGGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTtgt                                                                                                                            >  3:279200/1‑149 (MQ=255)
                atgaACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTtg                                                                                                                             <  2:370143/148‑1 (MQ=255)
                 tgaACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTtgt                                                                                                                            <  2:205230/148‑1 (MQ=255)
                   aaCCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACg                                                                                                                                                <  4:106031/126‑1 (MQ=255)
                   aaCCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACg                                                                                                                                                >  3:106031/1‑126 (MQ=255)
                    aCCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGccc                                                                                                                        <  7:138429/149‑1 (MQ=255)
                     ccGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGccc                                                                                                                        >  2:73771/1‑148 (MQ=255)
                       ggCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAg                                                                                                                     <  4:126019/149‑1 (MQ=255)
                       ggCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAg                                                                                                                     >  3:51357/1‑149 (MQ=255)
                                  ggCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGa                                                                                                          <  4:143738/149‑1 (MQ=255)
                                  ggCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCAAGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGa                                                                                                          <  2:33088/149‑1 (MQ=255)
                                    ccTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGt                                                                                                        <  2:382760/149‑1 (MQ=255)
                                        gCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTc                                                                                                      <  3:261083/147‑1 (MQ=255)
                                                 aCCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCAc                                                                                           <  4:325331/149‑1 (MQ=255)
                                                 aCCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCAc                                                                                           >  7:137684/1‑149 (MQ=255)
                                                 aCCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCAc                                                                                           >  1:99368/1‑149 (MQ=255)
                                                 aCCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCAc                                                                                           >  2:46520/1‑149 (MQ=255)
                                                    cAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGa                                                                                        <  2:288545/149‑1 (MQ=255)
                                                       cTTGCCAATAAGTTTAGCAGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCa                                                                                     <  8:137684/149‑1 (MQ=255)
                                                        ttGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAg                                                                                    >  4:373516/1‑149 (MQ=255)
                                                          gCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTAcc                                                                                                                 <  3:216029/118‑1 (MQ=255)
                                                          gCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTAcc                                                                                                                 >  4:216029/1‑118 (MQ=255)
                                                           ccAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGca                                                                                 >  4:46949/1‑149 (MQ=255)
                                                            cAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGca                                                                                 >  5:99125/1‑148 (MQ=255)
                                                              aTCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAg                                                                              >  3:100812/1‑149 (MQ=255)
                                                              aTCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAg                                                                              >  6:43544/1‑149 (MQ=255)
                                                              aTCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAg                                                                              >  6:49487/1‑149 (MQ=255)
                                                               tCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCAc                                                                                           >  5:311793/1‑135 (MQ=255)
                                                               tCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCAc                                                                                           <  6:311793/135‑1 (MQ=255)
                                                                   tttCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTg                                                                         >  2:342777/1‑149 (MQ=255)
                                                                   tttCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTg                                                                         <  3:373516/149‑1 (MQ=255)
                                                                       gccgcATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGtt                                                                      >  3:20868/1‑148 (MQ=255)
                                                                        ccgcATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGATTGAGTTCGGAGGGTTCACCGATCAGCCACAGCCGTGGTTTa                                                                    >  6:262350/1‑149 (MQ=255)
                                                                             tACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAATGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACCCgg                                                               <  2:285599/149‑1 (MQ=255)
                                                                             tACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACgg                                                               <  3:46949/149‑1 (MQ=255)
                                                                             tACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACgg                                                               <  7:391043/149‑1 (MQ=255)
                                                                                tCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGgcg                                                            >  8:119190/1‑149 (MQ=255)
                                                                                  gATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCGGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGgcgc                                                           <  1:342777/148‑1 (MQ=255)
                                                                                      gCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTg                                                      <  4:100812/149‑1 (MQ=255)
                                                                                         ttCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAg                                                   <  7:143631/149‑1 (MQ=255)
                                                                                         ttCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAg                                                   <  5:262350/149‑1 (MQ=255)
                                                                                         ttCACACCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCTCAGCCGTGGTTTACACGGGCGCTGTTGGAg                                                   <  4:279200/149‑1 (MQ=255)
                                                                                              cccATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTggg                                                                                                <  3:187365/99‑1 (MQ=255)
                                                                                              cccATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTggg                                                                                                >  4:187365/1‑99 (MQ=255)
                                                                                              cccATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGt                                                                          <  2:386927/121‑1 (MQ=255)
                                                                                              cccATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcag                                              >  8:49368/1‑149 (MQ=255)
                                                                                                 aTACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGTCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcagcca                                           <  5:43544/149‑1 (MQ=255)
                                                                                                 aTACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGcagcca                                           <  5:49487/149‑1 (MQ=255)
                                                                                                     ggCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCAt                                       >  8:309750/1‑149 (MQ=255)
                                                                                                         cgTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAg                                   <  1:46520/149‑1 (MQ=255)
                                                                                                         cgTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAg                                   <  6:392163/149‑1 (MQ=255)
                                                                                                         cgTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAg                                   <  6:99125/149‑1 (MQ=255)
                                                                                                              gCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCAc                              >  6:230975/1‑149 (MQ=255)
                                                                                                                cGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCAc                                                                                           >  5:313717/1‑86 (MQ=255)
                                                                                                                cGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCAc                                                                                           <  6:313717/86‑1 (MQ=255)
                                                                                                                          gCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAAcc                  <  3:68144/149‑1 (MQ=255)
                                                                                                                           cTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAAcc                  >  8:375136/1‑148 (MQ=255)
                                                                                                                            tGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGa                >  4:131229/1‑149 (MQ=255)
                                                                                                                               gcCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAggg                                                 >  2:211293/1‑113 (MQ=255)
                                                                                                                               gcCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAggg                                                 <  1:211293/113‑1 (MQ=255)
                                                                                                                               gcCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGca                                               >  3:345919/1‑115 (MQ=255)
                                                                                                                               gcCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGca                                               <  4:345919/115‑1 (MQ=255)
                                                                                                                                cccTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAgg                                                  <  6:264425/111‑1 (MQ=255)
                                                                                                                                cccTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAgg                                                  >  5:264425/1‑111 (MQ=255)
                                                                                                                                 ccTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATc                                      >  5:227098/1‑122 (MQ=255)
                                                                                                                                 ccTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATc                                      <  6:227098/122‑1 (MQ=255)
                                                                                                                                    gCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCg        >  7:143075/1‑149 (MQ=255)
                                                                                                                                     aagcTGGCAATGCGTTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGa                                    <  6:60915/117‑1 (MQ=255)
                                                                                                                                     aagcTGGCAATGCGTTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGa                                    <  6:67109/117‑1 (MQ=255)
                                                                                                                                     aagcTGGCAATGCGTTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGa                                    >  5:67109/4‑120 (MQ=255)
                                                                                                                                     aagcTGGCAATGCGTTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGa                                    >  5:60915/4‑120 (MQ=255)
                                                                                                                                      agcTGGCAATGCGTTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGa                >  5:161519/3‑139 (MQ=255)
                                                                                                                                      agcTGGCAATGCGTTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGa                <  6:161519/137‑1 (MQ=255)
                                                                                                                                      aCCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAg                                   <  5:190456/120‑1 (MQ=255)
                                                                                                                                      aCCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAg                                   >  6:190456/1‑120 (MQ=255)
                                                                                                                                        cTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCa     >  3:73228/1‑148 (MQ=255)
                                                                                                                                           gCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAg  <  4:73228/148‑1 (MQ=255)
                                                                                                                                           gCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAg  <  7:375136/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
CGAAGCGCTCGTCGTGATGAACCGGCAGTGCCAGGGCCTCGCCGATTTCACCCAACTTGCCAATCAGTTTCGCCGCATACTCGATGGCGTTCACCCCATACGGCGCGTACGCCGAATGGCACGCTGCGCCCTGCACCTGGCAACGCATTGCCAGCTTGCCCTTGTGCCCCAGTACCGGTTTGAGTTCGGTGGGTTCACCGATCAGGCACAGCCGTGGTTTACACGGGCGCTGTTGGAGGGCAGCCAGCATCGAGCGCACGCCCAAACAACCGACTTCTTCGTCATAG  >  NC_002947/4051755‑4052041

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: