Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,317,372 C→T 5.7% A348A (GCG→GCA trpB ← tryptophan synthase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,317,3720CT 5.7% 2488.1 / 81.8 749A348A (GCG→GCAtrpBtryptophan synthase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (312/394);  new base T (18/25);  total (330/419)
Fisher's exact test for biased strand distribution p-value = 8.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 763 total aligned reads displayed.

AGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATG  >  NC_000913/1317224‑1317519
                                                                                                                                                    |                                                                                                                                                   
aGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAgc                                                                                                                                                     >  1:454050/1‑149 (MQ=255)
  aTTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAgcgc                                                                                                                                                   >  1:4306070/1‑149 (MQ=255)
       cccTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGga                                                                                                                                              <  4:4388880/149‑1 (MQ=255)
        ccTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGgat                                                                                                                                             <  4:667910/149‑1 (MQ=255)
        ccTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGgat                                                                                                                                             >  1:5800890/1‑149 (MQ=255)
          tCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCgg                                                                                                                                                <  1:2718107/144‑1 (MQ=255)
           cGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGgatgat                                                                                                                                          <  4:475687/149‑1 (MQ=255)
           cGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGgatgat                                                                                                                                          <  3:4899236/149‑1 (MQ=255)
           cGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGgatgat                                                                                                                                          <  2:3662767/149‑1 (MQ=255)
            gTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATc                                                                                                                                         <  2:1496285/149‑1 (MQ=255)
             tGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATcc                                                                                                                                        <  4:811555/149‑1 (MQ=255)
             tGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATcc                                                                                                                                        <  3:6272148/149‑1 (MQ=255)
               cTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCt                                                                                                                                      <  1:1207874/149‑1 (MQ=255)
                tttCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCtt                                                                                                                                     >  1:477094/1‑149 (MQ=255)
                tttCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCtt                                                                                                                                     <  4:4925244/149‑1 (MQ=255)
                tttCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCtt                                                                                                                                     >  3:2129985/1‑149 (MQ=255)
                   cAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGt                                                                                                                                  <  3:6951875/149‑1 (MQ=255)
                    aaaaTATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGga                                                                                                                                              <  3:2854477/136‑1 (MQ=255)
                     aaaTATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGc                                                                                                                                >  2:5456219/1‑149 (MQ=255)
                       atatCTTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAg                                                                                                                              <  4:5549549/149‑1 (MQ=255)
                                 cGGTGCAGCTGTGTTTATCGCCCTGACCGGAAAGGTTAACAACGAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGGGGATTCCAGCGCCGGGATGATCCCTTCGTCCAGGCGCAGCggt                                                                                                                    >  2:119437/1‑147 (MQ=255)
                                 cGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGtt                                                                                                                    <  2:7056005/149‑1 (MQ=255)
                                  ggTGAAGATGTCTTTATCGCCGCGACCGTAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCTGGATGATCCCTTCGTGCAGGCACAGCGttt                                                                                                                   <  3:3417235/149‑1 (MQ=255)
                                  ggTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGttt                                                                                                                   <  2:5797263/149‑1 (MQ=255)
                                  ggTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGttt                                                                                                                   <  3:3684861/149‑1 (MQ=255)
                                   gTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGtttt                                                                                                                  <  3:2512981/149‑1 (MQ=255)
                                    tGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTg                                                                                                                 >  2:1962759/1‑149 (MQ=255)
                                     gAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCACGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGa                                                                                                                >  2:2899146/1‑149 (MQ=255)
                                      aaGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGTGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGaa                                                                                                               <  1:915923/149‑1 (MQ=255)
                                      aaGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGaa                                                                                                               <  4:6888419/149‑1 (MQ=255)
                                      aaGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGaa                                                                                                               <  2:4223219/149‑1 (MQ=255)
                                      aaGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGaa                                                                                                               <  1:2543756/149‑1 (MQ=255)
                                      aaGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGaa                                                                                                               <  3:6785493/149‑1 (MQ=255)
                                        gATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAgg                                                                                                             <  2:4605779/149‑1 (MQ=255)
                                         aTGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGc                                                                                                            <  3:19208/149‑1 (MQ=255)
                                              tttATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGTAGGCTTCaa                                                                                                       >  3:1215415/1‑149 (MQ=255)
                                               ttATCGCCGCGACCGGAAAGGTTAACTACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAg                                                                                                      <  1:5487168/149‑1 (MQ=255)
                                               ttATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGGTTTGAAGGCTTCAAg                                                                                                      >  1:3539862/1‑149 (MQ=255)
                                                tATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTTGTGCAGGCACAGCGTTTTGAAGGCTTCAAgg                                                                                                     >  3:894937/1‑149 (MQ=255)
                                                  tCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGc                                                                                                   <  3:1023997/149‑1 (MQ=255)
                                                    gccgcGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATc                                                                                                                                         >  4:3011928/1‑109 (MQ=255)
                                                     ccgcGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAg                                                                                                      >  3:5293605/1‑143 (MQ=255)
                                                     ccgcGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTtc                                                                                                <  4:5177489/149‑1 (MQ=255)
                                                      cgcgACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTtca                                                                                               <  2:6534507/149‑1 (MQ=255)
                                                      cgcgACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTtca                                                                                               <  3:3181552/149‑1 (MQ=255)
                                                            cGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGAGCCCt                                                                                                                                      <  4:5053823/104‑1 (MQ=255)
                                                             ggAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACTCATGGGCCAGGGCGTGGGAGGATTCCAGCTCCGAGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCgg                                                                                        >  1:5189985/1‑149 (MQ=255)
                                                             ggAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGTGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAg                                                                                                                              <  1:1557478/111‑1 (MQ=255)
                                                             ggAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTtcatc                                                                                             <  4:5812960/144‑1 (MQ=255)
                                                             ggAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTtcatc                                                                                             >  3:4304233/1‑144 (MQ=255)
                                                             ggAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCgg                                                                                        >  3:6130406/1‑149 (MQ=255)
                                                              gAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGt                                                                                       >  3:4155980/1‑149 (MQ=255)
                                                               aaaGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAg                                                                                                                        <  1:2361437/115‑1 (MQ=255)
                                                                  ggTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATa                                                                                   >  1:7210177/1‑149 (MQ=255)
                                                                   gTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAg                                                                                  >  2:4760145/1‑149 (MQ=255)
                                                                        ccaccaGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACg                                                                             >  3:1687001/1‑149 (MQ=255)
                                                                         caccaGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGt                                                                            >  1:3869968/1‑149 (MQ=255)
                                                                         caccaGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGt                                                                            >  4:2080020/1‑149 (MQ=255)
                                                                         caccaGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGt                                                                            >  4:6987177/1‑149 (MQ=255)
                                                                             aGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGTTTCCAGTGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATa                                                                                   <  1:6633309/138‑1 (MQ=255)
                                                                                aGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAgc                                                                     >  1:4062781/1‑149 (MQ=255)
                                                                                    gctCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGTGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCgg                                                                                        <  1:6049971/126‑1 (MQ=255)
                                                                                    gctCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTaa                                                                          <  1:1993008/140‑1 (MQ=255)
                                                                                     ctctTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTcatcat                                                                                           <  4:4087453/122‑1 (MQ=255)
                                                                                     ctctTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTc                                                                >  4:2282755/1‑149 (MQ=255)
                                                                                        tttATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGTCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAg                                                             >  4:3363470/1‑149 (MQ=255)
                                                                                        tttATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTtcatc                                                                                             >  3:2982206/1‑117 (MQ=255)
                                                                                         ttATCCGGGTATTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCTTCGGTAATAGACACGTAATCAGCGCGTCCAGt                                                            <  4:6098730/149‑1 (MQ=255)
                                                                                           aTCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGc                                                          >  1:3328762/1‑149 (MQ=255)
                                                                                             ccGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTg                                                        <  2:2143712/149‑1 (MQ=255)
                                                                                              cGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGtt                                                                                                                    <  2:6043069/88‑1 (MQ=255)
                                                                                              cGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGt                                                       >  2:5123674/1‑149 (MQ=255)
                                                                                                ggTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTa                                                     <  4:3810330/149‑1 (MQ=255)
                                                                                                 gTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTaa                                                    <  3:3912099/149‑1 (MQ=255)
                                                                                                  ttttCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAg                                                   <  4:2677578/149‑1 (MQ=255)
                                                                                                    ttCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGTCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGAt                                                 <  2:4392727/149‑1 (MQ=255)
                                                                                                    ttCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGAt                                                 <  1:2899146/149‑1 (MQ=255)
                                                                                                    ttCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGAt                                                 <  2:2721355/149‑1 (MQ=255)
                                                                                                       gcgcATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATAcg                                              >  3:6504862/1‑149 (MQ=255)
                                                                                                          catcatTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAg                                                                                                              >  4:1234763/1‑82 (MQ=255)
                                                                                                            tcatTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAAt                                                                         <  2:6367900/117‑1 (MQ=255)
                                                                                                               ttttCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGtgttgt                                      >  1:59635/1‑149 (MQ=255)
                                                                                                                  tCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTggg                                   <  1:6870190/149‑1 (MQ=255)
                                                                                                                     aCGCATGGTGCAGGGCGTGGGAGGATTCTAGTGCCGGGCTGATACCTTCGGGCAGGCACAGCGTTTTGAACGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGAGCGTTCAGTGCTGATAAGATAAGCGTGGTGTGTGCCg                                >  2:3462655/1‑149 (MQ=255)
                                                                                                                           gggCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGAc                          <  4:1658699/149‑1 (MQ=255)
                                                                                                                             gCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACgg                        <  1:1000312/149‑1 (MQ=255)
                                                                                                                             gCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACgg                        >  4:636198/1‑149 (MQ=255)
                                                                                                                             gCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACgg                        <  3:4563934/149‑1 (MQ=255)
                                                                                                                             gCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACgg                        >  2:738589/1‑149 (MQ=255)
                                                                                                                              ccAGGGCGTGGGAGGATTCCAGTGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACggg                       <  2:3977230/149‑1 (MQ=255)
                                                                                                                               cAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCt                                                         >  2:2597343/1‑114 (MQ=255)
                                                                                                                               cAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGa                      >  2:6804693/1‑149 (MQ=255)
                                                                                                                                aGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCatcatc                                                                                          >  1:1719774/1‑80 (MQ=255)
                                                                                                                                aGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCatcatc                                                                                          <  2:1719774/80‑1 (MQ=255)
                                                                                                                                 gggCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGaaa                    <  4:4215993/149‑1 (MQ=255)
                                                                                                                                   gCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATc                  <  3:5698057/149‑1 (MQ=255)
                                                                                                                                    cGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGTGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGc                                  >  3:5110957/1‑132 (MQ=255)
                                                                                                                                       gggAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGt              >  1:6324708/1‑149 (MQ=255)
                                                                                                                                        ggaggaTTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTc             >  4:5342965/1‑149 (MQ=255)
                                                                                                                                                   gcgcCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATg  >  3:5443721/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCCAGCGCCGGGATGATCCCTTCGTGCAGGCACAGCGTTTTGAAGGCTTCAAGGGCTTCATCATCGGTAATAGACACGTAATCAGCGCGTCCAGTGCTGTTAAGATACGCGTGTTGTGGGCCGACAGACGGGAAATCCAGTCCGGCGGAGATG  >  NC_000913/1317224‑1317519

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