Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,547,594 T→C 6.2% intergenic (‑249/‑95) pydP ← / → pydB NCS1 family transporter PydP/bifunctional D‑hydantoinase/dihydropyrimidinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,547,5940TC 6.2% 188.7 / 2.1 64intergenic (‑249/‑95)pydP/pydBNCS1 family transporter PydP/bifunctional D‑hydantoinase/dihydropyrimidinase
Reads supporting (aligned to +/- strand):  ref base T (32/28);  new base C (2/2);  total (34/30)
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.67e-01

ATCTTGTCAAGTTAGTCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGT  >  NC_002947/4547448‑4547736
                                                                                                                                                   |                                                                                                                                                
aTCTTGTCAAGTTAGTCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCtt                                                                                                                                                 <  2:360283/148‑1 (MQ=255)
     gTCAAGTTAGTCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCTTATTGATATCTACAACTTAAAAACCAGAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCTGCATCTAGCCTCGAATCTTGTCAg                                                                                                                                            <  6:32011/148‑1 (MQ=255)
      tCAAGTTAGTCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGa                                                                                                                                          <  4:52662/149‑1 (MQ=255)
                 aaaTACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGAtt                                                                                                                               >  3:347456/1‑149 (MQ=255)
                 aaaTACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCCTGTCAGTCCTGACAAGAtt                                                                                                                               >  5:383667/1‑149 (MQ=255)
                    tACCCGGACAACCCCGGCGCCCACTCGCCCGCTAATTTATTCCTATTGATATCAACACCTTAAAAACCACAAATTTAATCGGAAAATTTTCTTGCTGACTCCAGAATCCGCATCTAGCCTCGAATCCTGTCAGTCCTGACAAGAttcga                                                                                                                            <  8:112826/149‑4 (MQ=255)
                    tACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCCTGTCAGTCCTGACAAGAttcga                                                                                                                            <  6:383667/149‑4 (MQ=255)
                       cAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGa                                                                                                                                      >  5:306198/1‑136 (MQ=255)
                       cAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGa                                                                                                                                      <  6:306198/136‑1 (MQ=255)
                       cAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgc                                                                                                                         >  2:168418/1‑149 (MQ=255)
                       cAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgc                                                                                                                         >  7:310375/1‑149 (MQ=255)
                         ggAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgccc                                                                                                                       >  2:273093/1‑149 (MQ=255)
                         gcaaacCCCCGCCCCCCATCCCCCGGTAAATTTTTCCCAATTTATACAAAAACCTTAAAAACACAAAATTAAATGGCAAATTTTCCTGCCTTATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgccct                                                                                                                      <  4:360733/147‑1 (MQ=255)
                                         cATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATcc                                                                                                      <  1:347800/149‑1 (MQ=255)
                                              gcgggTTAAATTTTTTCTTTTTATAATAACAACTTTAAAACCACAAATTTTATACATTAAATTAGATATCTTATTCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGa                                                                                                  <  1:320723/149‑1 (MQ=255)
                                               cgcgTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGa                                                                                                                                          >  3:273686/1‑108 (MQ=255)
                                               cgcgTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGa                                                                                                                                          <  4:273686/108‑1 (MQ=255)
                                                 cgTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAAt                                                                                              <  2:172396/149‑1 (MQ=255)
                                                  gTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATg                                                                                             >  4:129524/1‑149 (MQ=255)
                                                  gTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAACGCCCTGCCCGCCCTGCATCCATGACAATg                                                                                             >  4:45556/1‑149 (MQ=255)
                                                       ttttttCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCGCTGCATCCATGACAATGCa                                                                                           >  2:267906/1‑146 (MQ=255)
                                                             ccTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCg                                                                                   <  1:273093/148‑1 (MQ=255)
                                                              cTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCgg                                                                                  <  8:310375/148‑1 (MQ=255)
                                                                   gATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTTTAAGGACTGACATGATTGACCGCCCTGCCCGCCCTGCATCAATGACAATGCAAAGAACCGGCCAAGa                                                                            >  3:209975/1‑149 (MQ=255)
                                                                    atatCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGAc                                                                           <  6:176439/149‑1 (MQ=255)
                                                                     tatCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGAcc                                                                          >  7:89894/1‑149 (MQ=255)
                                                                       tCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCg                                                                         <  2:513/148‑1 (MQ=255)
                                                                         aacaacTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGt                                                                       <  3:45556/148‑1 (MQ=255)
                                                                         aacaacTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGt                                                                       >  1:196398/1‑148 (MQ=255)
                                                                         aacaacTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTc                                                                      <  2:196398/149‑1 (MQ=255)
                                                                         aacaacTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATAAGCATCTACACACGAATCTTGTCAGGACTGACAGGATTGACCGCCCTTCACCACCTTCATCCATGACAAACCAAAGAACCGGACAAGACCgga                                                                       >  1:236673/1‑147 (MQ=255)
                                                                          acaacTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCa                                                                     <  8:126548/149‑1 (MQ=255)
                                                                            aacTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAg                                                                    >  8:332272/1‑148 (MQ=255)
                                                                               ttAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTTACAGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTg                                                                >  7:143659/1‑149 (MQ=255)
                                                                               ttAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTg                                                                >  2:140672/1‑149 (MQ=255)
                                                                                  aaaaCCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGc                                                               >  6:282399/1‑147 (MQ=255)
                                                                                     aCCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAgg                                                           >  7:317961/1‑148 (MQ=255)
                                                                                       cacaAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCGAAGAACCGGCCAAGACCGGTCAGCCTGCGa                                                             >  3:287775/1‑144 (MQ=255)
                                                                                       cacaAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGaaa                                                       <  4:347456/149‑1 (MQ=255)
                                                                                       cacaAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGaa                                                        >  7:98400/1‑148 (MQ=255)
                                                                                        acaAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGaaaa                                                      >  3:232388/1‑149 (MQ=255)
                                                                                            aTTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGc                                                  <  7:154616/149‑1 (MQ=255)
                                                                                                  tGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTcc                                            <  7:332272/149‑1 (MQ=255)
                                                                                                     ggAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGGCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCt                                          >  2:244351/1‑149 (MQ=255)
                                                                                                      gAAAATTTTCTTGCTGAATCCAGAATCATCATCTAGCCTCGAATCTTGTCAGGACTGACATGATTGACCGCCCTGCCCGCCCTGCTTCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCGCt                                          <  2:236673/148‑1 (MQ=255)
                                                                                                      gAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGACCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTg                                         >  2:250902/1‑149 (MQ=255)
                                                                                                             ttCTTGCTGAATGCAGAATCAGCATCTAGCCTCGAATGTTGTCAGGACTTACAGGATTGAACGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGACAAGACCGGTCAGCCTGCTAGGAAAAGGGCATGTCCCTGTTGATc                                   <  8:143659/148‑1 (MQ=255)
                                                                                                             ttCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATc                                   <  1:168418/148‑1 (MQ=255)
                                                                                                                  gCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGgc                             >  5:129588/1‑149 (MQ=255)
                                                                                                                   cTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGgcg                            >  6:107889/1‑149 (MQ=255)
                                                                                                                               tCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCGTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTa                <  1:162048/149‑1 (MQ=255)
                                                                                                                                cAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGTCCCGCCCTGCATCCATGACAATGCAAAGACCCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTa                >  1:91539/1‑149 (MQ=255)
                                                                                                                                cAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTCACCGCCCTGACCGCCCTGCATGAATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCACTGTTGATCCGTGCCGCCACCGTGGTTAc               >  4:274717/1‑149 (MQ=255)
                                                                                                                                       tAGGCTCGAATCTTGTCAGGAATGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGGAAAAGGAAAGAACCGGCCAAGACCGGTCAGCCGGGGAGGAAAACGGCATGGGCCGGTTGCTCCGTGGCGCCACCGTGGTTACCCACg          >  8:192287/1‑147 (MQ=255)
                                                                                                                                         gCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAaga      <  4:232388/149‑1 (MQ=255)
                                                                                                                                         gCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAaga      <  3:129524/149‑1 (MQ=255)
                                                                                                                                         gCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAaga      <  6:129588/149‑1 (MQ=255)
                                                                                                                                         gCCTCGAATCTTGTCACGACTGACAGGTATGACCGCCCAGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGTAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGaa        <  1:291339/148‑1 (MQ=255)
                                                                                                                                          ccTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCACTGTTGATCCGTGGCGCCACCGTGGTTACCCACGaa        >  4:184285/1‑146 (MQ=255)
                                                                                                                                           cTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAagaga    <  3:84705/149‑1 (MQ=255)
                                                                                                                                           cTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAagaga    <  2:328851/149‑1 (MQ=255)
                                                                                                                                           cTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAaga      >  3:264986/1‑147 (MQ=255)
                                                                                                                                            tCGAATCTTGTCAGGACTGACAGGATTAACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAagagag   >  2:44579/1‑149 (MQ=255)
                                                                                                                                            gtcaaTCCTGTCAGTCCTGACAAGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAagagag   >  4:344517/4‑149 (MQ=255)
                                                                                                                                              gAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGt  >  5:299004/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ATCTTGTCAAGTTAGTCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGT  >  NC_002947/4547448‑4547736

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