Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,153,208 T→C 7.7% intergenic (+28/‑92) fabF → / → pabC 3‑oxoacyl‑[acyl‑carrier‑protein] synthase II/4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,153,2080TC 7.7% 286.7 / 12.5 104intergenic (+28/‑92)fabF/pabC3‑oxoacyl‑[acyl‑carrier‑protein] synthase II/4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base T (46/50);  new base C (4/4);  total (50/54)
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.97e-01

Only 100 of 107 total aligned reads displayed.

AAGGTTGCGATCTGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCGGGCAACGC  >  NC_000913/1153059‑1153357
                                                                                                                                                     |                                                                                                                                                     
aaGGTTGCGATCTGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGt                                                                                                                                                       <  2:897364/150‑1 (MQ=255)
   gTTGCGATCTGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCg                                                                                                                                                    <  2:419568/150‑1 (MQ=255)
   gTTGCGATCTGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCCTATAAAAGGTCCg                                                                                                                                                    <  2:547899/150‑1 (MQ=255)
     tGCGATCTGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCt                                                                                                                                                  >  1:748589/1‑150 (MQ=255)
             ggATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGcc                                                                                                                                          >  1:1073839/1‑150 (MQ=255)
               aTTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGAGGGCCtt                                                                                                                                        >  1:1396738/1‑150 (MQ=255)
                tttCGTACCGCACGAAGCGCGGCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCttt                                                                                                                                       <  2:1252775/150‑1 (MQ=255)
                 ttCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGCCCGCAAGCGGACCtttt                                                                                                                                      <  2:1364986/150‑1 (MQ=255)
                  tCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGCCCGCAAGCGGACCtttta                                                                                                                                     <  2:1411679/150‑2 (MQ=255)
                   cGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCtttttt                                                                                                                                    >  2:883799/1‑150 (MQ=255)
                   cGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCtttttt                                                                                                                                    >  1:906303/1‑150 (MQ=255)
                   cGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGCCCGCAAGCGGACCttttat                                                                                                                                    <  1:1422778/150‑3 (MQ=255)
                    gTACCGCACGAAGCGCGTCAGGTTAGCGGAGTGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCTCTTGCGGGCCTTTTTTc                                                                                                                                   >  2:1377799/1‑150 (MQ=255)
                    gTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGCCCGCAAGCGGACCttttata                                                                                                                                   <  2:752333/150‑4 (MQ=255)
                      aCCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATCCACTCCGTGTAACTCCTACGGCTTCGGGGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTCTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCtt                                                                                                                                 <  1:290409/150‑1 (MQ=255)
                      aCCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCtt                                                                                                                                 >  1:889613/1‑150 (MQ=255)
                           aCGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCtt                                                                                                                            >  2:1465174/1‑150 (MQ=255)
                            cGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCttt                                                                                                                           >  2:1076797/1‑150 (MQ=255)
                            cGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCAACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCttt                                                                                                                           >  1:703816/1‑150 (MQ=255)
                             gAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCtttt                                                                                                                          >  2:1669885/1‑150 (MQ=255)
                               aGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTAt                                                                                                                        <  2:496219/150‑1 (MQ=255)
                                  gcgTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCTGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTcc                                                                                                                     <  1:607604/150‑1 (MQ=255)
                                  gcgTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTcc                                                                                                                     <  2:642910/150‑1 (MQ=255)
                                   ggTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCg                                                                                                                    <  2:1225950/149‑1 (MQ=255)
                                   cgTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCg                                                                                                                    >  1:302832/1‑150 (MQ=255)
                                   agTCAGGTTAGAGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGAGCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCg                                                                                                                    <  2:1099748/149‑1 (MQ=255)
                                       aGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTAGCGTGCCTTTTTTCTTATCTTTTATTCCGACtt                                                                                                                >  2:257939/1‑150 (MQ=255)
                                        ggTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTg                                                                                                               >  2:826446/1‑150 (MQ=255)
                                        ggTTAGCGGAATGGAATACACACTGTGAAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGATTGCGGGCCTTTTTTATTAGCTTTTATTCCGACtag                                                                                                               >  1:1503960/1‑148 (MQ=255)
                                         gTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGt                                                                                                              >  1:498139/1‑150 (MQ=255)
                                           tAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTc                                                                                                            <  1:682697/150‑1 (MQ=255)
                                            aGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTcc                                                                                                           <  1:1001880/150‑1 (MQ=255)
                                            aGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTcc                                                                                                           >  2:380797/1‑150 (MQ=255)
                                              cGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGt                                                                                                         <  2:669499/150‑1 (MQ=255)
                                              cGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGt                                                                                                         <  2:1121320/150‑1 (MQ=255)
                                              cGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGt                                                                                                         <  2:1451768/150‑1 (MQ=255)
                                               ggaatggaatACACTCTGTGTAACTCCTTCGGCTTCGGTGGCCCTAATGGTTCTTTGCTCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTGTTAGCTTTTATTCCGACTTGTTCCGtt                                                                                                        >  1:81088/1‑149 (MQ=255)
                                                gaatggaatACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAg                                                                                                       >  1:86844/1‑150 (MQ=255)
                                                  atggaatACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTg                                                                                                     <  1:970454/150‑1 (MQ=255)
                                                    ggaatACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTAAGACTTGTTCCGTAGTGaa                                                                                                   >  1:1187810/1‑150 (MQ=255)
                                                     gaatACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAAc                                                                                                  <  1:386966/150‑1 (MQ=255)
                                                      aatACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACa                                                                                                 <  2:574017/150‑1 (MQ=255)
                                                      aatACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACa                                                                                                 <  1:941033/150‑1 (MQ=255)
                                                       atACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACAt                                                                                                >  1:613804/1‑150 (MQ=255)
                                                        tACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACAtg                                                                                               <  2:636883/150‑1 (MQ=255)
                                                            ctctGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGAGGGCCTTATTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACAtgctgc                                                                                           >  2:320911/1‑150 (MQ=255)
                                                            cactGTTTAACTCCTTCGGCTTCGGAGGCACTAATGGTTCTTTGATCTTTACAAAGATCTATGTTGTCATTTTCCACCCAAATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACAtgctgc                                                                                           <  2:1630280/148‑1 (MQ=255)
                                                             tctGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGcc                                                                                          <  1:380797/150‑1 (MQ=255)
                                                               tgtgTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTGCGACTAGTTCCGTAGTGAACATGCTGCcac                                                                                        >  1:811746/1‑150 (MQ=255)
                                                               tgtgTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCcac                                                                                        <  2:748589/150‑1 (MQ=255)
                                                                gtgtAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCcaca                                                                                       <  1:826446/150‑1 (MQ=255)
                                                                gtgtAACTCCTTCGGCTTCGGTGGCACTAATGGTTATTTGTTCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCcaca                                                                                       <  2:743425/150‑1 (MQ=255)
                                                                  gtAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCTACTTGTTCCGTAGTGATCATGCGGCCacacg                                                                                     >  1:886582/1‑149 (MQ=255)
                                                                  gtAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACt                                                                                     >  1:596614/1‑150 (MQ=255)
                                                                     aCTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGTTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAAc                                                                                  >  1:1091315/1‑150 (MQ=255)
                                                                        ccTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAAt                                                                               >  1:829637/1‑150 (MQ=255)
                                                                         cttcggcttcggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAAtt                                                                              >  2:1350014/1‑150 (MQ=255)
                                                                                ttcggTTGCACTAATGGTTCTTTGATCTTTAAATAGATATAAGTTGTCATTTTCCTCCCTTATAAAAGGTCCGCTTGAGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTGCGTAGTGAACATGCTGCCACACTAACAATTCTCTGAt                                                                       <  1:273183/150‑1 (MQ=255)
                                                                                 tcggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATa                                                                      >  2:560466/1‑150 (MQ=255)
                                                                                     tGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGa                                                                  <  1:582826/150‑1 (MQ=255)
                                                                                     tGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGa                                                                  <  2:1603578/150‑1 (MQ=255)
                                                                                        cACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGcc                                                               <  1:318542/150‑1 (MQ=255)
                                                                                         aCTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCTACCATTAAAAAAGGTCCTCTTGCGGGCCTTTTTTCTTAGCTTATATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCg                                                              <  1:79910/150‑1 (MQ=255)
                                                                                         aCTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCg                                                              <  1:1377799/150‑1 (MQ=255)
                                                                                          cTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCgg                                                             <  2:1468002/150‑1 (MQ=255)
                                                                                          atAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCgg                                                             <  2:845049/149‑1 (MQ=255)
                                                                                           tAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACAtgctg                                                                                            >  1:1321856/1‑118 (MQ=255)
                                                                                           tAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACAtgctg                                                                                            <  2:1321856/118‑1 (MQ=255)
                                                                                              tGGTTCTTTGATATTTATAAAGATCTTAGTTGTCATTTTCCACCCTTATATAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATg                                                         <  2:1396738/150‑1 (MQ=255)
                                                                                                gTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTATCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGtt                                                       >  2:778849/1‑150 (MQ=255)
                                                                                                  tCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCt                                                     <  1:883799/150‑1 (MQ=255)
                                                                                                   cTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCtt                                                    >  2:485523/1‑150 (MQ=255)
                                                                                                   cTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCtt                                                    >  1:758521/1‑150 (MQ=255)
                                                                                                    tttGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCtta                                                   <  2:1073839/150‑1 (MQ=255)
                                                                                                    tttGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCtta                                                   >  2:1428206/1‑150 (MQ=255)
                                                                                                      tGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCttaat                                                 <  2:703816/150‑1 (MQ=255)
                                                                                                      tGAGCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCttaat                                                 <  1:1499368/150‑1 (MQ=255)
                                                                                                        aTCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCttaatta                                               <  1:1346555/150‑1 (MQ=255)
                                                                                                         tCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCttaattaa                                              <  2:1091315/150‑1 (MQ=255)
                                                                                                          atttAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGCGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAAc                                             <  1:381854/149‑1 (MQ=255)
                                                                                                           tttAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACg                                            <  1:1350014/150‑1 (MQ=255)
                                                                                                            ttAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACgg                                           >  1:1293351/1‑150 (MQ=255)
                                                                                                            ttAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACgg                                           <  1:1465174/150‑1 (MQ=255)
                                                                                                                        aaGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGa                               >  1:964904/1‑150 (MQ=255)
                                                                                                                          gTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTTAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAAt                             <  1:257939/150‑1 (MQ=255)
                                                                                                                          gTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAAt                             >  2:416124/1‑150 (MQ=255)
                                                                                                                            tGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGAAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCg                           >  1:402583/1‑150 (MQ=255)
                                                                                                                            tGTCATTTTCCACCCATATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTAATCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCg                           <  2:313638/150‑1 (MQ=255)
                                                                                                                                      cACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTaa                 >  2:1363122/1‑150 (MQ=255)
                                                                                                                                       aCCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAg                >  2:678934/1‑150 (MQ=255)
                                                                                                                                        cccTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGc               >  1:368872/1‑150 (MQ=255)
                                                                                                                                         taagaaaaaaGGCCCGCAAGCGGACCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCg              >  1:103969/7‑150 (MQ=255)
                                                                                                                                          aagaaaaaaGGCCCGCAAGCGGACCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGa             >  1:58040/6‑150 (MQ=255)
                                                                                                                                          aagaaaaaaGGCCCGCAAGCGGACCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGa             >  1:401761/6‑150 (MQ=255)
                                                                                                                                          aagaaaaaaGGCCCGCAAGCGGACCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGa             >  1:1409538/6‑150 (MQ=255)
                                                                                                                                            tataAAAGGTCCGCATGCGGGCCGTTTTTCTTAGCTTTTATTCCGACTTGTTCGGTTGTGAACATGCTCCGACACTAACAATTATATGATAAGGATCCTATAAGATCTTAATTATCGGTCAATAGCAGGAATCGCTGGCTGTAAGCGagc           >  2:260816/1‑148 (MQ=255)
                                                                                                                                             ataAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCg          >  2:286400/1‑150 (MQ=255)
                                                                                                                                             ataAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCg          <  2:889613/150‑1 (MQ=255)
                                                                                                                                                 aaGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCGGGCaa     <  2:134977/150‑1 (MQ=255)
                                                                                                                                                     tCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCGGGCAACGc  <  1:1428206/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGGTTGCGATCTGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCGGGCAACGC  >  NC_000913/1153059‑1153357

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: