Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 689,975 C→T 5.4% R221H (CGC→CAC)  PP_0591 ← adenosine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947689,9750CT 5.4% 284.4 / 2.5 92R221H (CGC→CAC) PP_0591adenosine deaminase
Reads supporting (aligned to +/- strand):  ref base C (43/44);  new base T (3/2);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTCACACCACGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGCTGGGTGGGT  >  NC_002947/689829‑690122
                                                                                                                                                  |                                                                                                                                                    
cTTTCACACCACGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGgcgcgc                                                                                                                                                    <  6:238067/149‑1 (MQ=255)
 tttCACACCACGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCa                                                                                                                                                   >  1:136205/1‑149 (MQ=255)
   tCACACCACGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCTAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCTGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCAc                                                                                                                                                  <  1:42537/148‑1 (MQ=255)
     acacCACGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGc                                                                                                                                                >  8:340727/1‑148 (MQ=255)
     acacCACGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGc                                                                                                                                                >  6:372938/1‑148 (MQ=255)
       accacGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGTGCGCGCCAt                                                                                                                                             <  5:141696/149‑1 (MQ=255)
       accacGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCAt                                                                                                                                             >  6:101438/1‑149 (MQ=255)
       accacGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCAt                                                                                                                                             >  2:32395/1‑149 (MQ=255)
         tacGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATgg                                                                                                                                           >  4:208124/2‑149 (MQ=255)
         cacGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATgg                                                                                                                                           >  6:244538/1‑149 (MQ=255)
                gAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATAc                                                                                                                                    <  1:338890/149‑1 (MQ=255)
                  gCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGATAGGTGTTGGACAGAGGGCACACAGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTAGTCCTCCATGGCGCGCACGCCATTGTCGATAAGt                                                                                                                                  >  7:234296/1‑149 (MQ=255)
                   cATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGtt                                                                                                                                 <  3:351560/149‑1 (MQ=255)
                    aTGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGttt                                                                                                                                <  2:40638/149‑1 (MQ=255)
                      gTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTg                                                                                                                               <  2:136205/148‑1 (MQ=255)
                      gTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTAGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTg                                                                                                                               >  3:287664/1‑148 (MQ=255)
                       tCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGa                                                                                                                              <  3:208124/148‑1 (MQ=255)
                       tCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGATCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGa                                                                                                                              >  6:227273/1‑148 (MQ=255)
                          aGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCtt                                                                                                                          <  7:92063/149‑1 (MQ=255)
                             aTGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTcag                                                                                                                       <  7:109167/149‑1 (MQ=255)
                              tgttgtGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTcagc                                                                                                                      <  4:72889/149‑1 (MQ=255)
                               gttgtGCTGGGTCATTCGGTCGAAGACGCAGAGCTTGGTGTTGGACAGCTGGCACTCCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGGTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTcagc                                                                                                                      <  7:277219/148‑1 (MQ=255)
                               gttgtGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGAt                                                                                                                                      >  6:59440/1‑132 (MQ=255)
                               gttgtGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGAt                                                                                                                                      <  5:59440/132‑1 (MQ=255)
                                  gtgCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGt                                                                                                                  >  5:112701/1‑149 (MQ=255)
                                     cTGGCTCATATGGTCGAACACGCAGAGCTTGGAGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTc                                                                                                                 >  3:146545/1‑147 (MQ=255)
                                        gCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAg                                                                                                              <  3:365119/147‑1 (MQ=255)
                                          tCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCt                                                                                                          >  5:63836/1‑149 (MQ=255)
                                             tatGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGc                                                                                                           <  3:230036/145‑1 (MQ=255)
                                             tatGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGc                                                                                                           >  4:230036/1‑145 (MQ=255)
                                             tatGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTcc                                                                                                       <  7:158743/149‑1 (MQ=255)
                                                 gTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAg                                                                                                    <  5:227273/148‑1 (MQ=255)
                                                 gTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGa                                                                                                   <  6:63836/149‑1 (MQ=255)
                                                     aaCACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTa                                                                                               >  3:127002/1‑149 (MQ=255)
                                                        acGCAGAGCTTGGTGTTGGACAGCGGGCTCACCGGGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTc                                                                                            <  8:188396/149‑1 (MQ=255)
                                                         cgctgttcTTGGTGTTGTTGAGCGGTCAGAGCGTTAGTGGGATCTGCTCGTCTTTGATGCGCTGCTTCAGGTGGTCGTTGTTGATGTCGCTTACGCCATTGTCTATATGTTTGATCTTCAGCAGGTCGAGGGCTTCCCTGATGTACTc                                                                                            <  8:234296/141‑1 (MQ=255)
                                                          gCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCCGATGTACTcgg                                                                                          >  7:355937/1‑149 (MQ=255)
                                                          gCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTcgg                                                                                          >  8:288583/1‑149 (MQ=255)
                                                          gCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTcgg                                                                                          <  2:182424/149‑1 (MQ=255)
                                                           cAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTcggg                                                                                         <  5:101438/149‑1 (MQ=255)
                                                            agagCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCGATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTcgggc                                                                                        >  5:342446/1‑149 (MQ=255)
                                                              agCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTcgggcgg                                                                                      <  2:372228/149‑1 (MQ=255)
                                                                    gtgtTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCcttc                                                                                >  1:56080/1‑149 (MQ=255)
                                                                    gtgtTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCcttc                                                                                <  2:56080/149‑1 (MQ=255)
                                                                         ggACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCcttc                                                                                <  7:178959/144‑1 (MQ=255)
                                                                         ggACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCcttc                                                                                >  8:178959/1‑144 (MQ=255)
                                                                         ggACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGc                                                                          >  1:62157/1‑149 (MQ=255)
                                                                          gACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCTGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGc                                                                          <  4:146545/148‑1 (MQ=255)
                                                                          gACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGc                                                                          >  2:216728/1‑148 (MQ=255)
                                                                                 ggCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATgg                                                                  <  5:372938/149‑1 (MQ=255)
                                                                                  gCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATggg                                                                 >  4:278440/1‑149 (MQ=255)
                                                                                  gCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATggg                                                                 <  7:288583/149‑1 (MQ=255)
                                                                                  gCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATggg                                                                 <  7:340727/149‑1 (MQ=255)
                                                                                      acCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCacc                                                             >  7:256428/1‑147 (MQ=255)
                                                                                      acCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAAc                                                             <  3:278440/149‑1 (MQ=255)
                                                                                      acCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAAc                                                             >  8:110921/1‑149 (MQ=255)
                                                                                      acCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAAc                                                             >  8:163420/1‑149 (MQ=255)
                                                                                          tGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCAGGCGGGCCTTCTTCGCCAGCATGGGCAACGGcg                                                         >  1:136818/1‑149 (MQ=255)
                                                                                           gAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGcga                                                        >  2:361979/1‑149 (MQ=255)
                                                                                            aGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGcgac                                                       <  6:112701/149‑1 (MQ=255)
                                                                                            aGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGcgac                                                       <  8:252826/149‑1 (MQ=255)
                                                                                                 taTCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGc                                                  <  1:309435/148‑1 (MQ=255)
                                                                                                       cTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCTGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCg                                             <  2:136818/148‑1 (MQ=255)
                                                                                                       cTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTc                                              <  4:249609/147‑1 (MQ=255)
                                                                                                       cTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTc                                              <  3:96938/147‑1 (MQ=255)
                                                                                                       cTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCACACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTCTCGCCAGCATGGGCAACGGCGACGAAGCCTTCg                                             <  1:50668/149‑1 (MQ=255)
                                                                                                        tcgtcgATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGcc                                                                                   >  8:346106/1‑110 (MQ=255)
                                                                                                        tcgtcgATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGcc                                                                                   <  7:346106/110‑1 (MQ=255)
                                                                                                          gtcgATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTAc                                         >  5:303982/1‑149 (MQ=255)
                                                                                                           tcgATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATCCGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATTGGCACAGGCGACGAAGCCTTCGCTACg                                        >  4:59304/1‑149 (MQ=255)
                                                                                                           tcgATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACg                                        >  8:275048/1‑149 (MQ=255)
                                                                                                              atgatgCGTTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGgc                                     >  3:9274/1‑149 (MQ=255)
                                                                                                              atgatgCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCGACGGgc                                     >  8:298323/1‑149 (MQ=255)
                                                                                                                gatgCGCTGCATCAGGCTCTTGTCCTCGATGGCGCGCACGCCATGGTAGATACGTTTGATCTTCAGCAGGTCGTGGGCTTCCGAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGTAACGGCGACGAAGCCTTCGCTACGGgcg                                    <  6:342446/148‑1 (MQ=255)
                                                                                                                gatgCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGgcgc                                   >  6:345972/1‑149 (MQ=255)
                                                                                                                  tgCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCgg                                 <  5:89736/149‑1 (MQ=255)
                                                                                                                   gcgcTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGt                                >  1:47630/1‑149 (MQ=255)
                                                                                                                    cgcTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCttcttc                                                                            >  3:2735/1‑104 (MQ=255)
                                                                                                                    cgcTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCttcttc                                                                            <  4:2735/104‑1 (MQ=255)
                                                                                                                    cgcTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTc                               <  1:216728/149‑1 (MQ=255)
                                                                                                                      cTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGa                             <  3:167356/149‑1 (MQ=255)
                                                                                                                              ggCGGTCGTCCTCGATGGCGCGCACTCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGc                     <  4:287664/149‑1 (MQ=255)
                                                                                                                               gcgcTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCt                    >  8:240240/1‑149 (MQ=255)
                                                                                                                               gcgcTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCt                    >  8:224805/1‑149 (MQ=255)
                                                                                                                                cgcTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTg                   <  7:275048/149‑1 (MQ=255)
                                                                                                                                 gcTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGAt                                                                                                  >  8:71597/1‑70 (MQ=255)
                                                                                                                                  cTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTgg                  >  2:70366/1‑148 (MQ=255)
                                                                                                                                     atcgaccaTGGCGTGCGCGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGaa                >  6:315551/8‑147 (MQ=255)
                                                                                                                                     atcgaccaTGGCGTGCGCGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGaa                <  5:315551/140‑1 (MQ=255)
                                                                                                                                           gATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGCtgg        >  2:60205/1‑149 (MQ=255)
                                                                                                                                           caTGGCGTGCGCGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGCtgg        >  3:372594/2‑149 (MQ=255)
                                                                                                                                           caTGGCGTGCGCGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGAAGTACTCGGGCGGGCCTTCTTCGCCAGCATCGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGCtgg        >  3:213027/2‑149 (MQ=255)
                                                                                                                                             tGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGCtgggt      <  1:70366/149‑1 (MQ=255)
                                                                                                                                                 gcgcACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGCTgggtgggt  >  1:199227/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTTTCACACCACGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGCTGGGTGGGT  >  NC_002947/689829‑690122

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