Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,261,0550AC44.8% 136.4 / 20.3 87V455G (GTG→GGG) tdcEpyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base A (41/7);  new base C (0/39);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 6.49e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTT  >  NC_000913/3260907‑3261205
                                                                                                                                                     |                                                                                                                                                       
gTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCAcc                                                                                                                                                       >  1:1119/1‑151 (MQ=255)
 tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCAccc                                                                                                                                                     >  1:409766/1‑151 (MQ=255)
  aGTGAATGATATTCAGCGCGCTGATTTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCTTCCAcccc                                                                                                                                                    <  2:931720/151‑1 (MQ=255)
  aGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCAcccc                                                                                                                                                    <  1:995624/151‑1 (MQ=255)
    tGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccgc                                                                                                                                                  >  2:970038/1‑151 (MQ=255)
     gAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccgcc                                                                                                                                                 >  2:837283/1‑151 (MQ=255)
       aTGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGt                                                                                                                                               >  2:1107872/1‑151 (MQ=255)
         gATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGtt                                                                                                                                              >  2:1079674/1‑150 (MQ=255)
              tCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTg                                                                                                                                         >  1:1187573/1‑150 (MQ=255)
                aGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCg                                                                                                                                       >  2:218325/1‑150 (MQ=255)
                aGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCg                                                                                                                                       >  1:1404697/1‑150 (MQ=255)
                   gcgcTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGt                                                                                                                                      >  1:983334/1‑148 (MQ=255)
                     gcTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGagc                                                                                                                                 >  1:1213469/1‑151 (MQ=255)
                            tACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGCTAATTGCGTAGAGCAGCGtt                                                                                                                           >  2:97531/1‑150 (MQ=255)
                              cTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTc                                                                                                                        >  1:663036/1‑151 (MQ=255)
                                 ccccGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGCCCTAAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGc                                                                                                                      <  1:150764/148‑1 (MQ=255)
                                    cGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                                                                                                  >  2:530369/1‑151 (MQ=255)
                                      ccagccggACAATGAAGTGATCGAGTCTGCCCATCACTTTTTCGTTAACCAGCCCGTCTTCCCTCCGGGGTGTTGTTTCCGCCACGACCAGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                                                                                                  <  1:893310/149‑1 (MQ=255)
                                       cgtcctccccATGAGGTGACCGGGGTTGCCCATCGCTTTGCCGTAACCCAACCCGCCGTCCTTCAGCGGTGCTGTGTGCGGCCCGCCTGGATTCTTCATCTCCTGTCCCCCGCCGCCGTTAAGTGCGTAGAGCAGCGTTTTCGCCAGGTTa                                                                                                               <  2:430256/143‑1 (MQ=255)
                                          ctcttccgttctgTGATCGAGGCTGTCCACCCTTTTGTCGTAATCCAGCACGTCGTCCACCAGCGGTGCTTTTTTCGCCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgc                                                                                                             <  1:76717/138‑1 (MQ=255)
                                          cgtggccATGAAGTGAGCGGGGCTGCCCCTCCTTTTGTCGTAACCCGGCCCGCGGTCCACCGGCGGTGGTTTTTTCGCCCCGCCGTGAGTCGTCAGTTTCCCGTCCCCGCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgc                                                                                                             <  2:789851/145‑1 (MQ=255)
                                          ccagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCCCGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgc                                                                                                             <  1:394621/150‑1 (MQ=255)
                                          ccagTCCATGAAGTGATCGAGGCTCTCCCTCCCTTTGTCGTACTCCACCACGGCGTCCACCAGCGGTGCTTTTTGCGGCCCGACCTGAATCTTCAGCTTTCCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgcg                                                                                                            <  2:547066/151‑1 (MQ=255)
                                          ccagTCCATGAAGTGACCGAGGCTGTCCATCCCTTTGTCGTAACCCAGCACGTCGCCCATCAGCGGTGCTTTTTTCGGCCCGACCTAAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTa                                                                                                               <  1:394549/148‑1 (MQ=255)
                                          ccagTCCAGGAATGGATCGAGTCTGCCCTTCTTTTTGTCGAACCCCAGCGCTTCGTCCACCAGGGTTGCTGTTTTCGGCCCGCCCTGAATCTTCAGTTTCCCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTa                                                                                                               <  2:644786/148‑1 (MQ=255)
                                           cagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgcgc                                                                                                           >  1:495826/1‑151 (MQ=255)
                                               ccATGAAGTGATCGGGGCTGTCCATCTTTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTCCGGCCCGACCTGATTTTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGt                                                                                                         <  1:951605/149‑1 (MQ=255)
                                                    aaGTGATCGGGGCTCTCCATCTCTTTGTCGTAACCCGGCCCGTCGTCCATCAGCGGTGTTGTTTTCGCCCCGCCCTTAATCTTCAGCTTTTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCa                                                                                                   <  1:1166564/150‑1 (MQ=255)
                                                    aaGTGATCGAGGCTGTCCATCCCTTTGTCGTAATCCAGCACGTCGTCCATCAGGGGTGCTTTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAc                                                                                                     <  1:195707/148‑1 (MQ=255)
                                                      ggaagcgggacgttCCATCGTTTTGTCGGACTCCACGCCGCCGCCCCTCAGGGGTGGGTTTGTCGGCCCGCCTGGAAGCTCCTTCTCCGCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTGAGCGCGTGCGCCCAAg                                                                                                <  2:79734/138‑1 (MQ=255)
                                                      gTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgc                                                                                                             >  2:76703/1‑138 (MQ=255)
                                                        gagcGAGGCTGTCCATCCCTTTGTCGTAACCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTAAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGaa                                                                                              <  1:1049326/148‑1 (MQ=255)
                                                        gATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGa                                                                                               >  2:1422174/1‑150 (MQ=255)
                                                        gATCGAGGCTGTCAATCACTTTGTCATAATCCACCTCGTCGTCCATCAGCGTTGCTTTTTTCGCCCCGACCTAAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAg                                                                                                <  1:412770/149‑1 (MQ=255)
                                                         aTCGAGGCTGTCCATCACTCTGACGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTACAGCAGCGTTTTCGCCCGGCTAGCGCGCGCACCAAAGAAc                                                                                             >  2:1084564/1‑151 (MQ=255)
                                                            ggggCTGTCCATCCCTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACt                                                                                            <  1:422820/147‑1 (MQ=255)
                                                            gAGGCTGTCCATCTCTTTGTCGTAACCCAGCACGTCGTCCATCAGCGGTGTTTTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTg                                                                                           <  1:1030678/150‑1 (MQ=255)
                                                             aGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGc                                                                                          >  1:1293248/1‑150 (MQ=255)
                                                              ggCTGTCCTTCCTTTTGTCGTAACCCACCCCGTCGTCCCTCAGCGGTGCTTTTTGCGCCCCGCCCTGAATCTGCGGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACt                                                                                            <  1:463957/147‑1 (MQ=255)
                                                              ggCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                                                                        >  2:311006/1‑151 (MQ=255)
                                                               gCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt                                                                                       >  1:855248/1‑151 (MQ=255)
                                                                cTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAttt                                                                                      >  1:321734/1‑151 (MQ=255)
                                                                cTGCCCATCATTTTGTCGTAACCCAGCCCGTCGTCCATCGGCGGTGCTGTTTTCGCCCGGACCGTATTCTTCAGCTTCGCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                                                                        <  2:1119/149‑1 (MQ=255)
                                                                 tGTCCACCTCTTTGTCGTAATCCGGCCCGTCGCCCCCCGGGGGTGCTTTTTGCGGCCCGCCTTGATTCTTCATCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt                                                                                       <  2:213336/149‑1 (MQ=255)
                                                                   tCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCt                                                                                   >  2:987927/1‑151 (MQ=255)
                                                                     cATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCtt                                                                                  >  2:744122/1‑150 (MQ=255)
                                                                      aTCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                                                                                 >  2:701714/1‑150 (MQ=255)
                                                                        cACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                                                                               >  1:562172/1‑150 (MQ=255)
                                                                               tCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCAcc                                                                        >  2:15854/1‑150 (MQ=255)
                                                                                cGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAt                                                                      >  1:1288277/1‑151 (MQ=255)
                                                                                  gaacccaccgcctcgccCATGAGGGTTCCTGTTTCCGCCCCGACATAATCCTGAAGTTCCTAGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAGAGAACTGCATTTGCTTACCAATCACCATTg                                                                    <  2:45540/139‑1 (MQ=255)
                                                                                  aaacccaccaCTTCGCCCATCAGGGGTGCTGTTTTCGCCCCGCCTGAATTCTCCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAt                                                                      <  2:623324/145‑1 (MQ=255)
                                                                                   agtCCAGCGCGTCCCCCACCCTCGGTGCTTTTTGCTGCCCGACCTGAATCTTCAGCTTCTCTCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTgg                                                                   <  2:1012339/149‑1 (MQ=255)
                                                                                    aTCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTggg                                                                  >  1:744659/1‑151 (MQ=255)
                                                                                     cccAGCCCGTCGTCCATCAGCGTTCCTTTTTTCGCCCCGACCTGAATCTCCATTTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTggg                                                                  <  2:452315/149‑1 (MQ=255)
                                                                                      ccAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGc                                                                 >  2:1099989/1‑150 (MQ=255)
                                                                                       cAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTg                                                               >  2:514078/1‑151 (MQ=255)
                                                                                         gCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGAc                                                             >  1:472083/1‑151 (MQ=255)
                                                                                            cgtcgtCCATCAGCGGTGTTTTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgca                                                          <  1:1423105/151‑1 (MQ=255)
                                                                                              tcccccATCGGCGGTGTTGTTTTCGCCCCGACCTGATTCTTCGTCTTCTCTTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgca                                                          <  2:566870/145‑1 (MQ=255)
                                                                                                  ccATCGGCGGTGTTGTTTCCGCCCCGACCTAAATCTTCTGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagca                                                       <  2:572247/148‑1 (MQ=255)
                                                                                                    aTCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCa                                                   >  2:256706/1‑150 (MQ=255)
                                                                                                        gCGGTGCTGTTTCCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgc                                              <  1:1079990/151‑1 (MQ=255)
                                                                                                          ggTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcg                                               <  1:610449/148‑1 (MQ=255)
                                                                                                          ggTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGt                                            >  2:129715/1‑151 (MQ=255)
                                                                                                            tgttGTTTTCGGCCCGACCTAAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTa                                           <  1:982774/147‑1 (MQ=255)
                                                                                                               tGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcg                                       >  1:28222/1‑151 (MQ=255)
                                                                                                                 ttttCGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcg                                       <  1:204103/149‑1 (MQ=255)
                                                                                                                 tttcCGGCCCGACCTGAATCTTCAGCTCCTCGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgt                                      <  1:1087425/150‑1 (MQ=255)
                                                                                                                  ttgcggcccggccTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtc                                     <  2:835627/147‑1 (MQ=255)
                                                                                                                   ttCGGCCCGCCCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcg                                    <  1:397059/150‑1 (MQ=255)
                                                                                                                   ttCGCCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcg                                    <  1:1037256/150‑1 (MQ=255)
                                                                                                                     cGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTg                                 >  1:233570/1‑151 (MQ=255)
                                                                                                                      ggCCCGCCCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCt                                  <  1:530544/149‑1 (MQ=255)
                                                                                                                      ggCCCGCCATGAATTTTCTTCTTCCCTTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTg                                 <  2:1038129/150‑1 (MQ=255)
                                                                                                                       gCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGtt                               >  1:587178/1‑151 (MQ=255)
                                                                                                                        cccGACCTGAATCTTCAGCTTCTCGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGtt                               <  1:837533/150‑1 (MQ=255)
                                                                                                                         ccGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTg                              >  2:933636/1‑150 (MQ=255)
                                                                                                                          cagactttactCTTCACCTTCGCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTg                              <  2:178454/140‑1 (MQ=255)
                                                                                                                              cTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc                         >  1:479092/1‑150 (MQ=255)
                                                                                                                                  aatttgcAGCTTCCCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc                         <  2:1023734/141‑1 (MQ=255)
                                                                                                                                     cTTCGGCTTCTCGTCCACCCCGCCGTTAATTGCGTATAGCAGCGTTTTCGCCAGGGTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCGTTGGGCTGACGCAGCAGGCAATCGCGTAATCATCGCTGTTGAAGTCAGTACGc                  >  2:1390020/1‑151 (MQ=255)
                                                                                                                                        cAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACg                   <  1:72050/147‑1 (MQ=255)
                                                                                                                                         aGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACg                   <  1:878331/146‑1 (MQ=255)
                                                                                                                                            ttCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCa                 <  2:364569/145‑1 (MQ=255)
                                                                                                                                                cGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCg         <  1:311509/149‑1 (MQ=255)
                                                                                                                                                   ccACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAtt    >  1:479778/1‑151 (MQ=255)
                                                                                                                                                     aCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAtttt  >  2:1020575/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
GTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTT  >  NC_000913/3260907‑3261205

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