Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 798,070 A→C V121G (GTC→GGC)  ybhA ← pyridoxal phosphate (PLP) phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913798,0700AC98.7% 301.8 / ‑3.9 78V121G (GTC→GGC) ybhApyridoxal phosphate (PLP) phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (35/42);  total (35/43)
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.18e-01

GCTGCAATTGCGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGACATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCAT  >  NC_000913/797924‑798217
                                                                                                                                                  |                                                                                                                                                   
gctgcAATTGCGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCAt                                                                                                                                                  >  2:763203/1‑150 (MQ=255)
gctgcAATTGCGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCAt                                                                                                                                                  <  1:704677/150‑1 (MQ=255)
gctgcAATTGCGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATg                                                                                                                                                 >  1:73327/1‑151 (MQ=255)
 ctgcAATTGCGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGc                                                                                                                                                >  1:1025247/1‑151 (MQ=255)
   gcAATTGCGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGccc                                                                                                                                              >  2:1053858/1‑151 (MQ=255)
    cAATTGCGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCg                                                                                                                                             <  2:436270/151‑1 (MQ=255)
         gCGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCg                                                                                                                                         <  2:794133/150‑1 (MQ=255)
         gCGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGCCTCAGGCAGAGAAGCGACTTGGGTGAAAGTCGGACGCAGTTCCGCCAGCAGGGTTTGCGCCCAGTTAGATGAGCCAATGCCATGCCCGGta                                                                                                                                          <  2:665365/149‑2 (MQ=255)
          cGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCggg                                                                                                                                       >  1:101038/1‑151 (MQ=255)
           ggCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGt                                                                                                                                      <  2:89864/151‑1 (MQ=255)
                 tCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCAt                                                                                                                                >  2:711540/1‑151 (MQ=255)
                  cATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCAt                                                                                                                                >  1:39698/1‑150 (MQ=255)
                   aTCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCAt                                                                                                                                >  1:198535/1‑149 (MQ=255)
                    tCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCAt                                                                                                                                >  2:1185396/1‑148 (MQ=255)
                     cGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGAc                                                                                                                            <  1:151199/151‑1 (MQ=255)
                     cGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGAc                                                                                                                            >  1:707844/1‑151 (MQ=255)
                     cGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGAc                                                                                                                            >  2:298179/1‑151 (MQ=255)
                      gTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGAcc                                                                                                                           <  2:585775/151‑1 (MQ=255)
                       tGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGAcc                                                                                                                           >  2:1204666/1‑150 (MQ=255)
                         cGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCAtt                                                                                                                        >  1:176274/1‑151 (MQ=255)
                         cGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCAtt                                                                                                                        >  2:1058968/1‑151 (MQ=255)
                          gTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCTGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTg                                                                                                                       <  2:1027256/151‑1 (MQ=255)
                                 cGAACTTCCATACGGCGTTAACTTGTTGCGTCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCatcatc                                                                                                                <  1:691994/151‑1 (MQ=255)
                                 cGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCatcatc                                                                                                                <  1:1116293/151‑1 (MQ=255)
                                  gAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCg                                                                                                               >  2:1081045/1‑151 (MQ=255)
                                   aaCTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCg                                                                                                               <  1:1064341/150‑1 (MQ=255)
                                       tCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGacat                                                                                                           <  1:890628/150‑1 (MQ=255)
                                       tCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGacat                                                                                                           <  1:904400/150‑1 (MQ=255)
                                         cATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGacatac                                                                                                         >  1:770152/1‑150 (MQ=255)
                                          aTACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGacatacat                                                                                                       <  1:151910/151‑1 (MQ=255)
                                             cGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGATGCTGTTCCGGCGGCAGGGCTTGCGCCCAGTTAGATGTGCGACTGCCATGCCCGGTCGGGGGCTCATAGACCATTGCATCATCGaaagacat                                                                                                       >  2:762436/1‑148 (MQ=255)
                                                cGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGAcc                                                                                                 >  2:729816/1‑151 (MQ=255)
                                                 gTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGACCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCg                                                                                                <  2:242061/151‑1 (MQ=255)
                                                     aCTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGaa                                                                                            <  2:722502/151‑1 (MQ=255)
                                                              gcCGTTTCAGCCAGAGAAGCGACTTGTGTGCAAGCCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCAATGCGCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGgtgt                                                                                    >  2:565467/1‑150 (MQ=255)
                                                                cGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGttc                                                                                  >  2:573740/1‑150 (MQ=255)
                                                                 gTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGttcat                                                                                <  1:644049/151‑1 (MQ=255)
                                                                 gTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGttcat                                                                                <  2:453969/151‑1 (MQ=255)
                                                                   ttCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGttcattc                                                                              >  1:968377/1‑151 (MQ=255)
                                                                     cAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGttcattca                                                                             <  1:805452/150‑1 (MQ=255)
                                                                     cAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGttcattc                                                                              >  1:17836/1‑149 (MQ=255)
                                                                        ccAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCa                                                                          >  1:1124934/1‑150 (MQ=255)
                                                                           gagaAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCAtctc                                                                      <  1:1180626/151‑1 (MQ=255)
                                                                                   aCTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGtt                                                              <  2:256403/151‑1 (MQ=255)
                                                                                     ttGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGc                                                            >  1:1025471/1‑151 (MQ=255)
                                                                                       gtgttAAAGGCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAg                                                          <  2:891059/151‑1 (MQ=255)
                                                                                          tgAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGc                                                       >  1:960011/1‑151 (MQ=255)
                                                                                           gAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCt                                                      <  2:550260/151‑1 (MQ=255)
                                                                                           gAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCt                                                      <  1:503381/151‑1 (MQ=255)
                                                                                           gAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGACATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCt                                                      <  2:878946/151‑1 (MQ=255)
                                                                                            aaaGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCtt                                                     >  1:78330/1‑151 (MQ=255)
                                                                                              aGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCttt                                                    >  1:398647/1‑150 (MQ=255)
                                                                                                 cGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTtaata                                                <  2:26446/151‑1 (MQ=255)
                                                                                                  ggACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTtaataa                                               <  1:1013758/151‑1 (MQ=255)
                                                                                                   gACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAAc                                              <  1:242148/151‑1 (MQ=255)
                                                                                                      gCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACggg                                           >  1:1161740/1‑151 (MQ=255)
                                                                                                       cTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGc                                          >  1:865477/1‑151 (MQ=255)
                                                                                                           tCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCAtt                                       <  1:42843/150‑1 (MQ=255)
                                                                                                           tCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTg                                      >  2:1001157/1‑151 (MQ=255)
                                                                                                            ccGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTgg                                     >  2:815245/1‑151 (MQ=255)
                                                                                                            ccGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTg                                      <  2:1094299/150‑1 (MQ=255)
                                                                                                            ccGGCGGCAGGGTTTGCGCCCAGTGAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTgg                                     <  2:25077/151‑1 (MQ=255)
                                                                                                             cggcggCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTggg                                    <  2:682988/151‑1 (MQ=255)
                                                                                                               gcggcAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGGCCGTGAATGTGGTGTTCATTCTGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGt                                   <  2:95748/150‑1 (MQ=255)
                                                                                                               gcggcAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTc                                  <  1:350744/151‑1 (MQ=255)
                                                                                                               gcggcAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTc                                  <  1:1160574/151‑1 (MQ=255)
                                                                                                                 ggcAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTc                                  <  1:546710/149‑1 (MQ=255)
                                                                                                                    aGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGAGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCtt                             >  2:94731/1‑151 (MQ=255)
                                                                                                                            cgcCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACg                     <  2:300594/151‑1 (MQ=255)
                                                                                                                              cccAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGt                   >  2:1184856/1‑151 (MQ=255)
                                                                                                                                cAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGttt                 <  2:101033/151‑1 (MQ=255)
                                                                                                                                 aGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGtttt                <  2:176266/151‑1 (MQ=255)
                                                                                                                                 aGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGcttt                >  1:538220/1‑151 (MQ=255)
                                                                                                                                   ttAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGttttt               <  2:155059/150‑1 (MQ=255)
                                                                                                                                    tAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGc             >  1:226165/1‑151 (MQ=255)
                                                                                                                                     aGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGc             <  2:387240/150‑1 (MQ=255)
                                                                                                                                       aTGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATg          <  2:1050698/151‑1 (MQ=255)
                                                                                                                                       aTGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATg          <  1:1185472/151‑1 (MQ=255)
                                                                                                                                               atgccatgccCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCAt  >  2:443525/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCTGCAATTGCGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGACATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCAT  >  NC_000913/797924‑798217

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