Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,436,151 A→G 6.5% Y305H (TAT→CAT)  PP_1257 ← 1‑pyrroline‑4‑hydroxy‑2‑carboxylate deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,436,1510AG 6.5% 176.2 / 2.6 61Y305H (TAT→CAT) PP_12571‑pyrroline‑4‑hydroxy‑2‑carboxylate deaminase
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base G (2/2);  total (30/31)
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.21e-01

CGACCCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGTTGAGCACCATCAAGTGCTTGTAATACAACACCAGGTCGCAGCCTT  >  NC_002947/1436007‑1436297
                                                                                                                                                |                                                                                                                                                  
cGACCCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAg                                                                                                                                                >  6:352856/1‑149 (MQ=255)
    cccTGTAGGAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGTGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGc                                                                                                                                            >  6:269258/1‑149 (MQ=255)
    cccTGTAGGAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGc                                                                                                                                            >  6:64933/1‑149 (MQ=255)
       tGTAGGAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCgg                                                                                                                                          >  1:138331/1‑148 (MQ=255)
           ggAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACa                                                                                                                                      <  5:352856/148‑1 (MQ=255)
                  gTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTAc                                                                                                                              >  6:46144/1‑149 (MQ=255)
                        ccGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGt                                                                                                                        >  5:150443/1‑149 (MQ=255)
                          gcgAACACCGGCGAAGTCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGTATGCCAGc                                                                                                                                               <  1:213036/124‑1 (MQ=255)
                          gcgAACACCGGCGAAGTCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGTATGCCAGc                                                                                                                                               >  2:213036/1‑124 (MQ=255)
                               cacCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATGCCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCAt                                                                                                                 >  7:289625/1‑149 (MQ=255)
                                      gAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGcgcg                                                                                                                                 <  2:314683/126‑1 (MQ=255)
                                      gAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGcgcg                                                                                                                                 >  1:314683/1‑126 (MQ=255)
                                       aaGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTgcgc                                                                                                         <  5:64933/149‑1 (MQ=255)
                                         cccGGTGCATGGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGc                                                                                                                   >  3:207171/2‑137 (MQ=255)
                                         cccGGTGCATGGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGc                                                                                                                   <  4:207171/136‑1 (MQ=255)
                                          ccGGTGCATGGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCt                                                                                                        >  2:305112/1‑147 (MQ=255)
                                                     cACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATGCCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAg                                                                                           <  8:289625/149‑1 (MQ=255)
                                                          ggTCGCCTGATTCGCGGGCACGCCCGCTACCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTgcg                                                                                                          >  1:198510/1‑129 (MQ=255)
                                                          ggTCGCCTGATTCGCGGGCACGCCCGCTACCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTgcg                                                                                                          <  2:198510/129‑1 (MQ=255)
                                                              gCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGt                                                                                  >  5:238104/1‑149 (MQ=255)
                                                                 tGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTc                                                                                                 <  4:218820/131‑1 (MQ=255)
                                                                       gcgGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTg                                                                           >  8:134062/1‑147 (MQ=255)
                                                                       gcgGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGaa                                                                         <  6:150443/149‑1 (MQ=255)
                                                                       gcgGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGa                                                                          <  1:305112/148‑1 (MQ=255)
                                                                         gggCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCg                                                                                                                     >  5:347557/1‑103 (MQ=255)
                                                                         gggCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCg                                                                                                                     <  6:347557/103‑1 (MQ=255)
                                                                         gggCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGt                                                                       >  8:331869/1‑149 (MQ=255)
                                                                         gggCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGt                                                                       >  8:302895/1‑149 (MQ=255)
                                                                         gggCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGt                                                                       >  3:217853/1‑149 (MQ=255)
                                                                          ggCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTTCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGaa                                                                         >  1:83271/1‑146 (MQ=255)
                                                                          ggCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGaa                                                                         <  2:83271/146‑1 (MQ=255)
                                                                           gCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGc                                                                     <  6:238104/149‑1 (MQ=255)
                                                                              cGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAgg                                                                  <  7:302895/149‑1 (MQ=255)
                                                                                cccGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAgcg                                                                                                                                   <  8:298781/82‑1 (MQ=255)
                                                                                cccGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAgcg                                                                                                                                   >  7:298781/1‑82 (MQ=255)
                                                                                cccGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCt                                                                >  7:242139/1‑149 (MQ=255)
                                                                                       ccACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACt                                                           <  7:134062/147‑1 (MQ=255)
                                                                                       ccACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCg                                                         >  6:240670/1‑149 (MQ=255)
                                                                                         acaGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGgt                                                       <  5:46144/149‑1 (MQ=255)
                                                                                         acaGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGgt                                                       <  7:331869/149‑1 (MQ=255)
                                                                                            gggAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTc                                                    >  7:76124/1‑149 (MQ=255)
                                                                                                    aGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGtt                                              <  3:57001/147‑1 (MQ=255)
                                                                                                    aGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGtt                                              >  4:57001/1‑147 (MQ=255)
                                                                                                      tGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGTTGAg                                           >  4:127905/1‑148 (MQ=255)
                                                                                                             tCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGTTGAGc                                          >  4:328366/1‑142 (MQ=255)
                                                                                                             tCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGTTGAGc                                          <  3:328366/142‑1 (MQ=255)
                                                                                                               aGGTGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGTTGAGCACCATCaa                                  <  5:269258/148‑1 (MQ=255)
                                                                                                               aGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGt                                                                             <  8:32731/105‑1 (MQ=255)
                                                                                                               aGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGt                                                                             >  7:32731/1‑105 (MQ=255)
                                                                                                                        tttGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCGGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGTACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGTTGAGCACCATCAAGTGCTTGTaa                        <  7:249666/149‑1 (MQ=255)
                                                                                                                         ttGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGTTGAGCACCATCAAGTGCTTGTaa                        <  8:242139/148‑1 (MQ=255)
                                                                                                                            cTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGTTGAGCACCATCAAGTGCTTGTAATac                     <  8:76124/148‑1 (MQ=255)
                                                                                                                                    aCCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGTTGAGCACCATCAAGTGCTTGTaa                        <  7:242995/137‑1 (MQ=255)
                                                                                                                                    aCCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGTTGAGCACCATCAAGTGCTTGTaa                        >  8:242995/1‑137 (MQ=255)
                                                                                                                                     ccAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGTTGAGCACCATCAAGTGCTTGTAATACAacac                <  5:334493/144‑1 (MQ=255)
                                                                                                                                     ccAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGAAGGCAGTACTCGGAGTCGACTTTGAGCACCATCAAGTGCTTGTAATACAACACCAga            >  6:334493/1‑147 (MQ=255)
                                                                                                                                      cAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGTTGAGCACCATCAAGTGCTTGTAATACAACACCAGGTc          >  5:173842/1‑149 (MQ=255)
                                                                                                                                        gcagtcATACCTGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGATCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGTTGAGCACCATCAAGTGCTTGTAATACAACACCAGGTCGc        <  6:285709/144‑1 (MQ=255)
                                                                                                                                        gCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGTTGAGCACCATCAAGTGCTTGTAATACAACACCAGGTCGc        <  2:67183/149‑1 (MQ=255)
                                                                                                                                         cTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGTTGAGCACCATCAAGTGCTTGTAATACAACACCAGGTCGCa       <  6:173842/149‑1 (MQ=255)
                                                                                                                                              aTACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGTTGAGCACCATCAAGTGCTTGTAATACAACACCAGGTCGCAGCCtt  <  4:217853/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CGACCCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGGGTCGCCTGATTCGCGGGCACGCCCGCTCCCACAGGGAGAGCAGTGGCTGCTCAGGCGATGTTTTGCTCGGCCGACCAGCTGGCATACCAGCTGCGGAACAGCGCGTACTGCTGTTCGGCATAGTTGCGCTGAGCGTCGGTCAGGACGTCGGTCTCGTTGAAGTGCAGGCTGTACTCGGTGTCGCCGTTGAGCACCATCAAGTGCTTGTAATACAACACCAGGTCGCAGCCTT  >  NC_002947/1436007‑1436297

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