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.6% 265.2 / ‑3.4 69V121G (GTC→GGC) ybhApyridoxal phosphate (PLP) phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (30/38);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01

GCTGCAATTGCGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGACATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCA  >  NC_000913/797924‑798216
                                                                                                                                                  |                                                                                                                                                  
gctgcAATTGCGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCAt                                                                                                                                                 <  1:375650/150‑1 (MQ=255)
  tgcAATTGCGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAatgccatgcc                                                                                                                                              >  2:955490/1‑151 (MQ=255)
   gcAATTGCGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAatgccatgcca                                                                                                                                             <  2:412666/151‑2 (MQ=255)
     aaTTGCGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCgg                                                                                                                                           >  2:222821/1‑151 (MQ=255)
        tGCGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTc                                                                                                                                         <  2:630293/150‑1 (MQ=255)
          cGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCggg                                                                                                                                      <  2:1049644/151‑1 (MQ=255)
           ggCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCgg                                                                                                                                       >  2:1119369/1‑149 (MQ=255)
              aGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGAGCCGTTGCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCt                                                                                                                                  >  1:685059/1‑151 (MQ=255)
               ggTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTc                                                                                                                                 <  2:426421/151‑1 (MQ=255)
                     cGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGAc                                                                                                                           <  1:930799/151‑1 (MQ=255)
                       tGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGAcc                                                                                                                          >  2:961248/1‑150 (MQ=255)
                                 cGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCatcatc                                                                                                               <  1:156994/151‑1 (MQ=255)
                                   aaCTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGa                                                                                                             <  1:428520/151‑1 (MQ=255)
                                   aaCTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGa                                                                                                             <  2:1242341/151‑1 (MQ=255)
                                        ccATACGGCGATAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGacata                                                                                                         >  1:764036/1‑150 (MQ=255)
                                         cATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGacataca                                                                                                       <  2:394591/151‑1 (MQ=255)
                                         cATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGacataca                                                                                                       <  1:1270545/151‑1 (MQ=255)
                                          aTACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGacatacat                                                                                                      >  2:932659/1‑151 (MQ=255)
                                            aCGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATc                                                                                                     <  1:300990/150‑1 (MQ=255)
                                            aCGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCa                                                                                                    <  1:297193/151‑1 (MQ=255)
                                                cGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGAcc                                                                                                >  1:822149/1‑151 (MQ=255)
                                                     aCTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGa                                                                                            >  2:721510/1‑150 (MQ=255)
                                                        tgttgCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAAt                                                                                          <  1:977200/149‑1 (MQ=255)
                                                              gcCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGtt                                                                                  >  2:959336/1‑151 (MQ=255)
                                                                cGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGttc                                                                                 >  1:342540/1‑150 (MQ=255)
                                                                     cAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGttcattc                                                                             >  1:535730/1‑149 (MQ=255)
                                                                         cAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCAtc                                                                       >  1:871993/1‑151 (MQ=255)
                                                                          agagAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCAtct                                                                      <  1:113645/151‑1 (MQ=255)
                                                                          agagAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGCGCGAATGCCATGCCCGGTCGGGGGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCAtct                                                                      >  2:143125/1‑151 (MQ=255)
                                                                                  gACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGt                                                              >  1:385395/1‑151 (MQ=255)
                                                                                   aCTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGtt                                                             <  1:727529/151‑1 (MQ=255)
                                                                                   aCTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGtt                                                             <  2:387366/151‑1 (MQ=255)
                                                                                      tgtgtgAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATTAGTTGc                                                           <  2:1212632/150‑1 (MQ=255)
                                                                                      tgtgtgAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGc                                                           >  1:760595/1‑150 (MQ=255)
                                                                                      tgtgtgAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGCGCGAATGCCATGCCCGGCCGGGGGCTCATAGACCATTGCATCATCGCCATACATCAGACCGTGAATGTGGCGTTCATTCAGCATCTCAATGAGTTGc                                                           >  1:355040/1‑150 (MQ=255)
                                                                                       gtgtgAAAGTCGGACGCTGTTCCGGCGGCATGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCa                                                          <  1:1117404/150‑1 (MQ=255)
                                                                                       gtgtgAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCa                                                          <  1:908109/150‑1 (MQ=255)
                                                                                       gtgtgAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAg                                                         <  1:462166/151‑1 (MQ=255)
                                                                                         gtgAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAggg                                                       >  2:361422/1‑151 (MQ=255)
                                                                                         gtgAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGACATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAgg                                                        <  1:214885/150‑1 (MQ=255)
                                                                                             aaGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCttt                                                   <  2:292771/151‑1 (MQ=255)
                                                                                                tCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTtaat                                                >  1:740157/1‑151 (MQ=255)
                                                                                                 cGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTtaat                                                <  2:691459/150‑1 (MQ=255)
                                                                                                  ggACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTtaata                                               >  2:435642/1‑150 (MQ=255)
                                                                                                   gACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAAc                                             >  1:1111680/1‑151 (MQ=255)
                                                                                                        tGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGc                                         >  1:689678/1‑150 (MQ=255)
                                                                                                        tGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGc                                         >  1:350116/1‑150 (MQ=255)
                                                                                                           tCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCAtt                                      <  1:271052/150‑1 (MQ=255)
                                                                                                                cggcAGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTcc                                <  1:955432/151‑1 (MQ=255)
                                                                                                                    aGGGTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCtt                            >  2:86018/1‑151 (MQ=255)
                                                                                                                       gTTTGCGCCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCa                         <  1:341261/151‑1 (MQ=255)
                                                                                                                            cgcCCAGTTAGATGTGCGAATGCCATGCCCGTTCGGGTGATCATAGACCATTGCATCATCGACATAGATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACg                    <  1:51748/151‑1 (MQ=255)
                                                                                                                            cgcCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACg                    >  1:229897/1‑151 (MQ=255)
                                                                                                                             gcCCAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACgg                   >  2:991158/1‑151 (MQ=255)
                                                                                                                              cccAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGt                  >  1:65935/1‑151 (MQ=255)
                                                                                                                              cccAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGt                  >  2:425008/1‑151 (MQ=255)
                                                                                                                              cccAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGt                  <  2:428817/151‑1 (MQ=255)
                                                                                                                                cAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGttt                <  1:806446/151‑1 (MQ=255)
                                                                                                                                cAGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGttt                <  2:277254/151‑1 (MQ=255)
                                                                                                                                 aGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGtttt               <  1:732666/151‑1 (MQ=255)
                                                                                                                                 aGTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGtttt               <  2:693326/151‑1 (MQ=255)
                                                                                                                                  gTTAGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGttttt              >  1:1200324/1‑151 (MQ=255)
                                                                                                                                     aGATGTGTGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCa           <  2:140369/151‑1 (MQ=255)
                                                                                                                                     aGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGc            >  2:682376/1‑150 (MQ=255)
                                                                                                                                     aGATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGc            <  2:1146427/150‑1 (MQ=255)
                                                                                                                                      gATGTGCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCa           <  1:952942/150‑1 (MQ=255)
                                                                                                                                         gtgCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGAt       <  2:110051/151‑1 (MQ=255)
                                                                                                                                           gCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGCGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATa      <  1:278274/150‑1 (MQ=255)
                                                                                                                                           gCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATaa     <  2:350124/151‑1 (MQ=255)
                                                                                                                                           gCGAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATa      >  2:237116/1‑150 (MQ=255)
                                                                                                                                             gAATGCCATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATc   <  2:503944/151‑1 (MQ=255)
                                                                                                                                               atgccatgccCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCa  <  2:1129980/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCTGCAATTGCGGCAGGTCATCGTGCGTCAGGGCGAACTTCCATACGGCGTTAACTTGTTGCGCCGTTTCAGCCAGAGAAGCGACTTGTGTGAAAGTCGGACGCTGTTCCGGCGGCAGGGTTTGCGCCCAGTTAGATGTGCGAATGACATGCCCGGTCGGGTGCTCATAGACCATTGCATCATCGACATACATCAGACCGTGAATGTGGTGTTCATTCAGCATCTCAATGAGTTGCAGGGCTTTAATAACGGGCATTGGGTCCGCTTCCAGCACGGTTTTTGCATGATAATCA  >  NC_000913/797924‑798216

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