Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,153,208 T→C 5.2% 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 5.2% 495.0 / 9.7 150intergenic (+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 (69/73);  new base C (4/4);  total (73/77)
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.69e-01

Only 100 of 151 total aligned reads displayed.

GTTGCGATCTGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCGGGCAACG  >  NC_000913/1153062‑1153356
                                                                                                                                                  |                                                                                                                                                     
gTTGCGATCTGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTcc                                                                                                                                                     <  1:137082/149‑1 (MQ=255)
gTTGCGATCTGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAAATCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTcc                                                                                                                                                     <  4:429327/149‑1 (MQ=255)
  tGCGATCTGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGc                                                                                                                                                   <  3:257760/149‑1 (MQ=255)
      aTCTGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGc                                                                                                                                               >  4:359301/1‑149 (MQ=255)
       tCTGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCg                                                                                                                                              >  3:56191/1‑149 (MQ=255)
           gATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGcc                                                                                                                                          >  3:107764/1‑149 (MQ=255)
           gATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGcc                                                                                                                                          <  1:521527/149‑1 (MQ=255)
            aTTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGCCCGCAAGCGGACC‑t                                                                                                                                        <  3:484669/149‑1 (MQ=255)
                   aCCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCt                                                                                                                                 >  2:1028280/1‑149 (MQ=255)
                   aCCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCt                                                                                                                                 <  3:22180/149‑1 (MQ=255)
                    ccGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCtt                                                                                                                                <  1:1001679/149‑1 (MQ=255)
                     cGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGCCCGCAAGCGGACC‑ttttataagg                                                                                                                               <  1:676603/149‑7 (MQ=255)
                        aCGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCt                                                                                                                            <  4:554131/149‑1 (MQ=255)
                         cGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCtt                                                                                                                           <  3:473823/149‑1 (MQ=255)
                          gAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCttt                                                                                                                          <  2:482044/149‑1 (MQ=255)
                            aGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTa                                                                                                                               <  1:243410/142‑1 (MQ=255)
                            aGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTa                                                                                                                               >  2:243410/1‑142 (MQ=255)
                                cgTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTcc                                                                                                                    <  3:975750/149‑1 (MQ=255)
                                 gTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCg                                                                                                                   >  4:856150/1‑149 (MQ=255)
                                  tCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTTCTTAGCTTTTATTCCg                                                                                                                   >  1:14837/1‑149 (MQ=255)
                                  tCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGa                                                                                                                  <  3:36920/149‑1 (MQ=255)
                                    aGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTGGCTTTTATTCCGACt                                                                                                                <  4:1034665/149‑1 (MQ=255)
                                      gTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTg                                                                                                              >  4:151342/1‑149 (MQ=255)
                                         aGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTc                                                                                                           <  4:817601/149‑1 (MQ=255)
                                         aGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTc                                                                                                           <  3:359301/149‑1 (MQ=255)
                                              aatggaatACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAg                                                                                                      >  2:590381/1‑149 (MQ=255)
                                              aatggaatACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGCCCGCAAGCGGACC‑tttta                                                                                                                                    <  1:909850/119‑2 (MQ=255)
                                                   aatACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAAc                                                                                                 >  3:1013882/1‑149 (MQ=255)
                                                   aatACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAAc                                                                                                 >  1:542559/1‑149 (MQ=255)
                                                    atACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACa                                                                                                <  2:646053/149‑1 (MQ=255)
                                                    atACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACa                                                                                                >  1:578259/1‑149 (MQ=255)
                                                       cacTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACAtgc                                                                                             <  1:775799/149‑1 (MQ=255)
                                                       cacTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACAtgc                                                                                             <  3:856150/149‑1 (MQ=255)
                                                        acTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACAtgct                                                                                            <  3:971663/149‑1 (MQ=255)
                                                         ctctGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCAACTTGTTCCGTAGTGAACAtgctg                                                                                           >  3:949709/1‑149 (MQ=255)
                                                          tctGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACa                                                                                                >  1:567345/1‑143 (MQ=255)
                                                          tctGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACa                                                                                                <  2:567345/143‑1 (MQ=255)
                                                          tctGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACAtgctgc                                                                                          >  2:33595/1‑149 (MQ=255)
                                                          tctGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACAtgc                                                                                             <  1:435474/146‑1 (MQ=255)
                                                          tctGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACAtgc                                                                                             >  2:435474/1‑146 (MQ=255)
                                                            tgtgTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCca                                                                                        >  4:519071/1‑149 (MQ=255)
                                                               gtAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCacac                                                                                     >  3:464630/1‑149 (MQ=255)
                                                                 aaCTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGc                                                                                                                             <  3:473941/107‑1 (MQ=255)
                                                                  aCTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTaa                                                                                  <  1:854201/149‑1 (MQ=255)
                                                                   cTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAAc                                                                                 <  3:50500/149‑1 (MQ=255)
                                                                   cTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAAc                                                                                 >  1:822088/1‑149 (MQ=255)
                                                                         cggcttcggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTTTTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATtct                                                                           <  3:616318/149‑1 (MQ=255)
                                                                          cgcttcggTTGCACTAATGGATCTTTGATCTTTAAAAAGATCTATGTTGTCATTTTCCACCCTTATAGAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTCAACATGCTTCCACACTAACAATtctc                                                                          <  2:449857/148‑1 (MQ=255)
                                                                           gcttcggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCAGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTctct                                                                         <  2:995738/149‑1 (MQ=255)
                                                                             ttcggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGa                                                                       >  3:880208/1‑149 (MQ=255)
                                                                                ggTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATaa                                                                    <  4:107731/149‑1 (MQ=255)
                                                                                  tGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAgg                                                                  >  1:31472/1‑149 (MQ=255)
                                                                                      aCTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGa                                                                       >  2:75487/1‑140 (MQ=255)
                                                                                      aCTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGa                                                                       <  1:75487/140‑1 (MQ=255)
                                                                                      aCTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGcc                                                              >  2:667342/1‑149 (MQ=255)
                                                                                      aCTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGcc                                                              >  1:370069/1‑149 (MQ=255)
                                                                                       cTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCg                                                             <  4:44830/149‑1 (MQ=255)
                                                                                          aTGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTa                                                          >  3:785125/1‑149 (MQ=255)
                                                                                          aTGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTa                                                          >  4:895471/1‑149 (MQ=255)
                                                                                           tGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCCTTTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTa                                                          <  2:14837/149‑1 (MQ=255)
                                                                                                 tttGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCtt                                                   <  3:734613/149‑1 (MQ=255)
                                                                                                 tttGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCtt                                                   >  2:765476/1‑149 (MQ=255)
                                                                                                   tGATCTTTAACAAGATGTAAGTTGTCATTTTCTTCCCTTAGAGAAGGTCTGCTTGCGGGCT‑TTTTTTCTTAGCGTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCATACTATCAATTCTCTGATATGGAGCCGTTATGTTCttaa                                                 <  2:85194/149‑1 (MQ=255)
                                                                                                     aTCTTTAAATAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATtctc                                                                          >  4:449710/1‑122 (MQ=255)
                                                                                                     aTCTTTAAATAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATtctc                                                                          <  3:449710/122‑1 (MQ=255)
                                                                                                     aTCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACtt                                                                                                               <  4:474969/85‑1 (MQ=255)
                                                                                                     aTCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGa                                                                       >  1:697618/1‑125 (MQ=255)
                                                                                                     aTCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGa                                                                       <  2:697618/125‑1 (MQ=255)
                                                                                                      tCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATtctc                                                                          >  3:513710/1‑121 (MQ=255)
                                                                                                            aaaaGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTc                                        <  1:1028280/149‑1 (MQ=255)
                                                                                                             aaaGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCa                                       >  2:26936/1‑149 (MQ=255)
                                                                                                                gATCTAAGTTGTCATTTTCCGCCCTTGTAAAAGGTCCGGTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACCATTCTCTAATAAGGAGCCGCTATGTTCTTAATTAACGGTCAttg                                    <  2:142839/149‑3 (MQ=255)
                                                                                                                  tCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGc                                  >  1:1070791/1‑149 (MQ=255)
                                                                                                                    tAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAg                                <  1:604842/149‑1 (MQ=255)
                                                                                                                    tAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAg                                >  3:267075/1‑149 (MQ=255)
                                                                                                                    tAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAg                                <  1:1084966/149‑1 (MQ=255)
                                                                                                                      aGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGa                              <  3:853132/149‑1 (MQ=255)
                                                                                                                       gTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGaa                             <  2:279561/149‑1 (MQ=255)
                                                                                                                          ttCATTTCCCACCTTTATAAAGGGTCCGCTTGCGGGCC‑TTTTTTCTTGGCTTTTTTTCCGACTGTTTCCTTGGTGACCATGCGCCCACACT‑ACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGc                         <  3:229485/148‑1 (MQ=255)
                                                                                                                           tCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAGCGGTCATAAGCAGGAATCGc                         <  4:500894/149‑1 (MQ=255)
                                                                                                                            cATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCt                        <  1:26936/149‑1 (MQ=255)
                                                                                                                             aTTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCttaa                                                 >  2:718374/1‑123 (MQ=255)
                                                                                                                             aTTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCttaa                                                 <  1:718374/123‑1 (MQ=255)
                                                                                                                             aTTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTg                       >  3:13347/1‑149 (MQ=255)
                                                                                                                                    aCCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTa                                                          <  3:16837/107‑1 (MQ=255)
                                                                                                                                    aCCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTa                                                          >  4:16837/1‑107 (MQ=255)
                                                                                                                                    aCCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTaa                >  1:391850/1‑149 (MQ=255)
                                                                                                                                    aCCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGAGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTaa                >  4:740548/1‑149 (MQ=255)
                                                                                                                                        ttATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGa            >  2:226616/1‑149 (MQ=255)
                                                                                                                                          ataAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCTACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATc          >  2:362131/1‑149 (MQ=255)
                                                                                                                                           taAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCg         <  3:519071/149‑1 (MQ=255)
                                                                                                                                            aaaaGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCgg        >  3:763953/1‑149 (MQ=255)
                                                                                                                                            aaaaGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAAAAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCgg        >  2:880194/1‑149 (MQ=255)
                                                                                                                                            aaaaGGCCCGCAAGCGGACC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCa                                 >  1:518404/1‑124 (MQ=255)
                                                                                                                                            aaaaGGCCCGCAAGCGGACC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCa                                 <  2:518404/124‑1 (MQ=255)
                                                                                                                                            aaaaGGCCCGCAAGCGGACC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCgg        >  3:202432/1‑149 (MQ=255)
                                                                                                                                             aaaGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCggg       <  1:669362/149‑1 (MQ=255)
                                                                                                                                             aaaGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCggg       <  3:476230/149‑1 (MQ=255)
                                                                                                                                                 gTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCAGATAAGGAGCTGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCGGGCAAc   <  2:868684/149‑1 (MQ=255)
                                                                                                                                                  tCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCGGGCAACg  >  4:754715/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTTGCGATCTGGATTTCGTACCGCACGAAGCGCGTCAGGTTAGCGGAATGGAATACACTCTGTGTAACTCCTTCGGCTTCGGTGGCACTAATGGTTCTTTGATCTTTAAAAAGATCTAAGTTGTCATTTTCCACCCTTATAAAAGGTCCGCTTGCGGGCC‑TTTTTTCTTAGCTTTTATTCCGACTTGTTCCGTAGTGAACATGCTGCCACACTAACAATTCTCTGATAAGGAGCCGGTATGTTCTTAATTAACGGTCATAAGCAGGAATCGCTGGCAGTAAGCGATCGGGCAACG  >  NC_000913/1153062‑1153356

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