Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,547,594 T→C 8.5% 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 8.5% 191.2 / 6.6 70intergenic (‑249/‑95)pydP/pydBNCS1 family transporter PydP/bifunctional D‑hydantoinase/dihydropyrimidinase
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base C (3/3);  total (34/36)
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.76e-01

ATCTTGTCAAGTTAGTCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGA  >  NC_002947/4547448‑4547732
                                                                                                                                                  |                                                                                                                                           
aTCTTGTCAAGTTAGTCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCtt                                                                                                                                            <  1:5075/148‑1 (MQ=255)
 tCTTGTCAAGTTAGTCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCCTg                                                                                                                                           >  3:76700/1‑148 (MQ=255)
                 aaaTACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGAtt                                                                                                                          >  1:77157/1‑149 (MQ=255)
                    tACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCCTGTCAGTCCTGACAAGAttcga                                                                                                                       <  4:76700/149‑4 (MQ=255)
                      ccAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATAAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCCTGTCAGTCCTGACAAGAttcgag                                                                                                                      >  1:41232/1‑144 (MQ=255)
                      ccAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATAAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCCTGTCAGTCCTGACAAGAttcgag                                                                                                                      <  2:41232/148‑5 (MQ=255)
                       cAGGAAAACCCCGGCGCCCATTCGCGAGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTATTGCTGAATCCAGAATCAGCTTCTAGCCTAGAATCTTGTAAGGTATGACATGATTTAccgc                                                                                                                    >  8:279580/1‑149 (MQ=255)
                       cAGGAAAACCCCGGCGCCCATACGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgc                                                                                                                    >  5:42825/1‑149 (MQ=255)
                        aGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgcc                                                                                                                   <  2:71300/149‑1 (MQ=255)
                        aGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCCTGTCAGTCCTGACAAGAttcga                                                                                                                       >  8:48862/1‑142 (MQ=255)
                        aGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCCTGTCAGTCCTGACAAGAttcga                                                                                                                       <  7:48862/145‑4 (MQ=255)
                                cccGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACt                                                                                                                                   >  6:88312/1‑125 (MQ=255)
                                cccGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACt                                                                                                                                   <  5:88312/125‑1 (MQ=255)
                                cccGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgccctgccc                                                                                                             <  4:122066/147‑1 (MQ=255)
                                cccGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgccctgccc                                                                                                             >  3:122066/1‑147 (MQ=255)
                                   ggCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgccctgcccgccct                                                                                                       <  6:58359/149‑1 (MQ=255)
                                     cgcCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGa                                                                                                                                     <  6:67528/118‑1 (MQ=255)
                                      gcCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCa                                                                                                    >  3:299840/1‑149 (MQ=255)
                                        ccATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCATCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCAt                                                                                                   >  1:81635/1‑148 (MQ=255)
                                              gcgcgTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGa                                                                                             <  6:42825/148‑1 (MQ=255)
                                               cgcgTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGAc                                                                                                                                    <  6:325111/109‑1 (MQ=255)
                                                              cTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGc                                                                            >  5:178360/1‑149 (MQ=255)
                                                               tATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGc                                                                            >  2:120086/1‑148 (MQ=255)
                                                               tATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCATGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGc                                                                            >  7:259115/1‑148 (MQ=255)
                                                                   gATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgcc                                                                                                                   >  2:187602/1‑106 (MQ=255)
                                                                   gATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgcc                                                                                                                   <  1:187602/106‑1 (MQ=255)
                                                                    atatCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGaa                                                                                 <  2:150735/138‑1 (MQ=255)
                                                                    atatCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGaa                                                                                 >  1:150735/1‑138 (MQ=255)
                                                                    atatCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGAc                                                                      >  1:204106/1‑149 (MQ=255)
                                                                       tCAACAACTTAAAAACCACAAATTAAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCgg                                                                   >  7:219260/1‑149 (MQ=255)
                                                                         aacaacTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgccctgcccgccct                                                                                                       <  2:265835/111‑1 (MQ=255)
                                                                         aacaacTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGt                                                                  <  2:204106/148‑1 (MQ=255)
                                                                         aacaacTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAAGCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGAccgccctgcccgccct                                                                                                       >  1:265835/1‑111 (MQ=255)
                                                                          acaacTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGt                                                                  >  5:211010/1‑147 (MQ=255)
                                                                              cTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAAt                                                                                         <  2:119530/120‑1 (MQ=255)
                                                                              cTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCt                                                            >  7:199465/1‑149 (MQ=255)
                                                                                 aaaaaCCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCATGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGc                                                          >  7:260712/1‑148 (MQ=255)
                                                                                 aaaaaCCACAAATTTAATGGGAAAATTTTCTTGCTGAATACAGAATCACCATCTAGCCTCGAATCTTGTCAGGACTGACACGATTGCCCGCCCTGACCGCCCTCCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCCGCCTGc                                                          >  4:194959/1‑148 (MQ=255)
                                                                                  aaaaCCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGa                                                        <  4:299840/149‑1 (MQ=255)
                                                                                  aaaaCCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGa                                                        <  6:178360/149‑1 (MQ=255)
                                                                                    aaCCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAgg                                                      <  2:81635/149‑1 (MQ=255)
                                                                                         caAATTTAAAGGGAAAATTTTCTTGCTGAATCCAGAATCATCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAAc                                                 <  1:135226/149‑1 (MQ=255)
                                                                                           aaTTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACgt                                               <  7:279580/149‑2 (MQ=255)
                                                                                            aTTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGc                                              <  8:199465/149‑1 (MQ=255)
                                                                                             attaaatggAAAATTTTCTTGCTGAATCCATAATCAGCATCTAGCCTCGAATCTTGTCAGGACTTACATGATTGAACGCCCTGCCCGCCAATGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGc                                              <  1:274786/142‑1 (MQ=255)
                                                                                               tAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCAt                                            <  6:211010/148‑1 (MQ=255)
                                                                                                   gggAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGa                                                        <  8:130004/132‑1 (MQ=255)
                                                                                                   gggAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGa                                                        >  7:130004/1‑132 (MQ=255)
                                                                                                     gAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTg                                                           >  1:208021/1‑127 (MQ=255)
                                                                                                     gAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTg                                                           <  2:208021/127‑1 (MQ=255)
                                                                                                     gAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTg                                     >  4:86943/1‑149 (MQ=255)
                                                                                                      aaatTTTCCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGt                                    <  1:120086/149‑1 (MQ=255)
                                                                                                       aaaTTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGtt                                   >  2:157849/1‑149 (MQ=255)
                                                                                                              cTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCg                             <  1:157849/148‑1 (MQ=255)
                                                                                                               ttGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTg                           >  1:40781/1‑149 (MQ=255)
                                                                                                                tGCTGATTCTGGATTCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCTCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTgg                          <  8:267628/149‑1 (MQ=255)
                                                                                                                tGCTGATTCTGGATTCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTgg                          >  7:267628/1‑149 (MQ=255)
                                                                                                                 gCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTgg                          <  3:86943/148‑1 (MQ=255)
                                                                                                                 gCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGAAAATGCAAAGAACCGGCCAAGACAGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTgg                          <  3:194959/148‑1 (MQ=255)
                                                                                                                  cTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGgcg                        >  7:330025/1‑149 (MQ=255)
                                                                                                                    gAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGt                            >  2:192030/1‑143 (MQ=255)
                                                                                                                    gAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGt                            <  1:192030/143‑1 (MQ=255)
                                                                                                                         cAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGa                                                        >  2:316539/1‑110 (MQ=255)
                                                                                                                                aGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTAcc          <  2:77157/149‑1 (MQ=255)
                                                                                                                                 gCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTAccc         <  2:40781/149‑1 (MQ=255)
                                                                                                                                     cTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTgg                          <  7:322935/128‑1 (MQ=255)
                                                                                                                                     cTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTgg                          >  8:322935/1‑128 (MQ=255)
                                                                                                                                        gCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAaga  <  8:259115/149‑1 (MQ=255)
                                                                                                                                          cTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAaga  >  4:68783/1‑147 (MQ=255)
                                                                                                                                           tCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGa                                                        >  5:73818/1‑92 (MQ=255)
                                                                                                                                           tCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGa                                                        <  6:73818/92‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
ATCTTGTCAAGTTAGTCAAATACCAGGAAAACCCCGGCGCCCATTCGCGCGTTAATTTTTTCCTATTGATATCAACAACTTAAAAACCACAAATTTAATGGGAAAATTTTCTTGCTGAATCCAGAATCAGCATCTAGCCTCGAATCTTGTCAGGACTGACAGGATTGACCGCCCTGCCCGCCCTGCATCCATGACAATGCAAAGAACCGGCCAAGACCGGTCAGCCTGCGAGGAAAACGGCATGTCCCTGTTGATCCGTGGCGCCACCGTGGTTACCCACGAAGA  >  NC_002947/4547448‑4547732

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: