Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,547,602 A→C 6.1% intergenic (‑257/‑87) 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,6020AC 6.1% 175.7 / 2.1 61intergenic (‑257/‑87)pydP/pydBNCS1 family transporter PydP/bifunctional D‑hydantoinase/dihydropyrimidinase
Reads supporting (aligned to +/- strand):  ref base A (31/26);  new base C (2/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

TCAAGTTAGTCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGTTACCC  >  NC_002947/4547454‑4547741
                                                                                                                                                     |                                                                                                                                             
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                                                                                                                                               <  4:273686/108‑1 (MQ=255)
                                         cgcgTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGa                                                                                                                                               >  3:273686/1‑108 (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           <  6:129588/149‑1 (MQ=255)
                                                                                                                                   gCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAaga           <  3:129524/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)
                                                                                                                                               gTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGTTAccc  <  8:89894/146‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
TCAAGTTAGTCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGAGAGTTACCC  >  NC_002947/4547454‑4547741

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: