Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,501 T→A 15.0% intergenic (+56/+31) ilvC → / ← ppiC ketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,959,5010TA15.0% 347.7 / 15.6 140intergenic (+56/+31)ilvC/ppiCketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Reads supporting (aligned to +/- strand):  ref base T (29/83);  new base A (20/2);  total (51/89)
Fisher's exact test for biased strand distribution p-value = 1.16e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

GAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACC  >  NC_000913/3959352‑3959648
                                                                                                                                                     |                                                                                                                                                   
gAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACt                                                                                                                                                     >  2:2900572/1‑150 (MQ=255)
  aCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCACACCCCGCCCCACTCCCACAGGGGGGGGGAGAAAAAAcaca                                                                                                                                                   >  1:384982/1‑150 (MQ=255)
     aaGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                                                                                                                <  1:1952011/150‑1 (MQ=255)
      aGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGc                                                                                                                                               <  1:2743857/150‑1 (MQ=255)
         gATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCtt                                                                                                                                            <  1:1639864/150‑1 (MQ=255)
          aTTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                                                                                                                           <  1:2747981/150‑1 (MQ=255)
            tCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                                                                                                                         <  2:1818580/150‑1 (MQ=255)
            tCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACAACCCCACGGAGATGGAGTAAAAACACAACGCCAtctc                                                                                                                                         >  2:577189/1‑150 (MQ=255)
             cGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                                                                                                                        <  1:1049306/150‑1 (MQ=255)
              gCAGCCATGAGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                                                                                                       <  1:1719078/150‑1 (MQ=255)
               cAGCGATGCGATTGAGCTGGTAGGTAAGAAACTGCGCGGCTATAGGACAGATATGAAACGTATTGCTGTGGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGTGTGGGAGAAAACACTCATGGCCTTCTCCTg                                                                                                                                      <  2:1734930/150‑1 (MQ=38)
               cAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTCAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCGTCTCCTg                                                                                                                                      >  1:2661286/1‑150 (MQ=255)
               cAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                                                                                                                      <  1:1113308/150‑1 (MQ=255)
                  ccATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCGACCCACACCCACACGGAGAGGGCGAAAACACACAAGGCCTTCGCCTGgag                                                                                                                                   >  1:378378/1‑150 (MQ=255)
                    aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                                 <  1:1298104/150‑1 (MQ=255)
                     tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAg                                                                                                                                <  2:119372/150‑1 (MQ=255)
                     tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCACAACCCCAGGGAGAGGGAGAAAAAAAACAAGGCCATCTACGGGAGAAg                                                                                                                                >  1:1750734/1‑150 (MQ=255)
                      gCGATTGAGCAGGTAGGTCAGAAACTGCGCGGCTATATGACAGATATGAAACGTAGTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTGCTCCTGGAGAAgg                                                                                                                               <  2:1964510/150‑1 (MQ=255)
                       cGATTGAGCTGGTAGGTAAGAAACTGCGCGACTATATGACAGATATGACACTTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGc                                                                                                                              <  1:765437/150‑1 (MQ=38)
                       cGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                                                                              <  1:1279729/150‑1 (MQ=255)
                       cGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGGTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                                                                              <  1:2144170/150‑1 (MQ=255)
                        gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGcc                                                                                                                             <  2:223705/150‑1 (MQ=255)
                         aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATCTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                                                                            <  1:2453050/150‑1 (MQ=255)
                         aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCt                                                                                                                            <  1:102344/150‑1 (MQ=255)
                             aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCCCCCACAACCACAGGGAGAGGGAGCAAAAAAACAAAGCCATCTCCAGGAGTAGGCCGTGCt                                                                                                                        >  2:2683443/1‑150 (MQ=255)
                              gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCACGGCCTTCTCCTGGAGAAGGCCTTGCTa                                                                                                                       <  1:2861935/150‑1 (MQ=255)
                              gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                                                                                                                       <  1:856419/150‑1 (MQ=255)
                              gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                                                                                                                       <  1:2861900/150‑1 (MQ=255)
                                 ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                                                                                                                    <  2:1078840/150‑1 (MQ=255)
                                 ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                                                                                                                    <  1:2518605/150‑1 (MQ=255)
                                    aggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCAATCCCACAGGGGGAGGGGGAAAACAAACAAGGCCCTATCCTGGTGAAGGCCGTGCTATTAttt                                                                                                                 >  1:938576/1‑150 (MQ=37)
                                      gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                                                                                                               <  2:1935422/150‑1 (MQ=255)
                                       taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                                                                                                              <  1:2235767/150‑1 (MQ=255)
                                       taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                                                                                                              <  1:258482/150‑1 (MQ=255)
                                        aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGGTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCCCCCCCACCCCCCCGGCGAGGGAGAAAACACACAAGGCCCTCGCCAGGCGCAGGCCGGGCGATTAGGTGCgg                                                                                                             >  2:2352483/1‑150 (MQ=255)
                                            aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa                                                                                                         <  1:2587126/150‑1 (MQ=255)
                                            aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa                                                                                                         <  1:1514288/150‑1 (MQ=255)
                                            aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa                                                                                                         <  1:2761867/150‑1 (MQ=255)
                                            aaCTGCGCGGCGATATGACAGAGATGAAACGTATTGCTGGTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCGCCCCCAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa                                                                                                         <  2:1489241/150‑1 (MQ=38)
                                             aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg                                                                                                        <  2:900628/150‑1 (MQ=255)
                                             aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCGACCCTCACCCACAGGGAGAGGGACAAAACACTCCAGGCCATATCCAGGAGAAGGCCGTGCTATTAGTTGCCGTACAg                                                                                                        >  2:25737/1‑150 (MQ=255)
                                               tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCa                                                                                                      <  1:492050/150‑1 (MQ=255)
                                                    ggCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCTCTCCCACCGGGAGAGGGAGAAAAAACACAATGCCTTCACCAGGGGGATGCCGTTCTAATAGTTGCGGTAACGCACCATa                                                                                                 >  2:819043/1‑150 (MQ=25)
                                                      cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                                                                               <  1:1950176/150‑1 (MQ=255)
                                                      cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                                                                               <  2:2858227/150‑1 (MQ=255)
                                                         tctGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgat                                                                                            <  1:1293954/148‑1 (MQ=255)
                                                         gatctcAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgat                                                                                            <  2:646668/145‑1 (MQ=255)
                                                            gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgtg                                                                                         <  1:193321/150‑1 (MQ=255)
                                                            gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgtg                                                                                         <  1:543184/150‑1 (MQ=255)
                                                               aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCACCCCCCCACCCACACCCCCCGGGAGTGGGGGAAAAAAATCGACGCCATAAGCGGGGGAAGGCGGTGCTGGTATTTGCGGTACAGCACCGTAatggtgagtaa                                                                                      >  1:646668/1‑145 (MQ=255)
                                                                   atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCGATACCACAGGGAGAGGGAGAAAACACTCAAGGCCGTCTCCGGGGGCAGGCCGTGCTAGTAGGTGCGGTACAGCACCGTAATTCTGTGATTTCCg                                                                                  >  2:2047829/1‑150 (MQ=25)
                                                                   atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCACCCCCACCCACACCCAAACGGAGAGGGAGAAAAAACACAAGGCCATCTCCTGGAGAAGGCCTTGCTAATAGTTTCGGTACAGCACCATAAAGATGTGATATACg                                                                                  >  2:2047823/1‑150 (MQ=255)
                                                                      aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                                                                               <  1:1896799/150‑1 (MQ=255)
                                                                      aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCACAGGGAGACGGAGAAAACACTCAAGGCCATCTCCTGGTGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                                                                               >  2:1861197/1‑150 (MQ=255)
                                                                      aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCACCCCGCCCCACAACCACACGGGGAGGGAGCAAACACACAATGCCCTGTACGGGGGCAGGCCGTGGAATTAGTTTCGGGGCACCACCGTAATGATGTGgtatacccac                                                                               >  1:2819499/1‑144 (MQ=255)
                                                                        aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACACAAGGCCTTCTCCGGGAGAAGGCCGTGCTATTAGTTGCGGTACAGCACCATAATGATGTGATATACGAACTg                                                                             >  2:1225334/1‑150 (MQ=38)
                                                                         cGTAGTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                            <  1:2072040/150‑1 (MQ=255)
                                                                            aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCGACCCACACCCACAGGGAGAGGGAGAAAAAACACAAGGCCATATCCAGGAGAAGGCCATGCTATTATTTTCGGTAAAGCACCATAATTATTTGATATACGAAATGGgtg                                                                         >  1:2648315/1‑150 (MQ=255)
                                                                              tGCTGTTGCGGGTTCAGTGCGCGCTGATGCCCACACCCCCACCCTCAACCACAGGGAGAGGGAGAAAACACACAAGGCCGTGGCCTGGGGAACGCCTTGCTCTTAGTTTCGGTACAGCCACATAATGTTTTGATATCCCAACCGGGtgtg                                                                       >  2:2760349/1‑150 (MQ=19)
                                                                                      cGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc                                                               <  2:2012073/150‑1 (MQ=255)
                                                                                      cGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACACCCACCCGGCGAGGGAGAAAACAAACAAGGCCATCTCCGGGACAAGGCCATGCGATTAGTTGCGGTACACCCCCAAAATGATGTGAGATCCGAACAGGGTGTGCCGCCGGc                                                               >  2:1465580/1‑150 (MQ=255)
                                                                                      cGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACACCCACCCGGCGAGGGAGAAAAAAAACAAGGCCATAAACAGGAGAAGGCCATGCTATTAGTTTAGGTACAACACCATAAAGATGTTATATACGAAAAGGGTGTGCACCGGGc                                                               >  2:1047866/1‑150 (MQ=255)
                                                                                        ggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGTGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTAGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCg                                                             <  1:1454877/150‑1 (MQ=255)
                                                                                           tAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCGTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTc                                                          <  2:853891/150‑1 (MQ=255)
                                                                                            aaGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCg                                                         <  2:1780566/150‑1 (MQ=255)
                                                                                                  gcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa                                                   <  2:1091828/150‑1 (MQ=255)
                                                                                                  gcgcTGAGGCCCTCACCCCGACCCTCTCCCACCGGGAGAGGGAGAAAACACGCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa                                                   <  1:2303969/150‑1 (MQ=255)
                                                                                                          gCCCTCCCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCGCCTGGAGCAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGa                                           <  1:314334/150‑1 (MQ=255)
                                                                                                          gCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGa                                           <  1:2234351/150‑1 (MQ=255)
                                                                                                              tCACCCCGACCCTCTCCCACAGGGCGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAag                                       >  2:1586152/1‑150 (MQ=255)
                                                                                                                aCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAagag                                     <  1:2086135/150‑1 (MQ=255)
                                                                                                                 ccccGACCCTCTCCCCCAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCGGGCGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAGGGGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgaga                                    >  2:1366835/1‑150 (MQ=255)
                                                                                                                   ccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaa                                  <  1:1693279/150‑1 (MQ=255)
                                                                                                                    cGACCCTCTCCCACAGGGAGAGGGAGAAACCACCCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGCTATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGGACCGGACAAGAGaaaa                                 <  2:1946369/150‑1 (MQ=255)
                                                                                                                        ccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTTGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc                             >  1:1424661/1‑150 (MQ=255)
                                                                                                                        ccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc                             >  1:844839/1‑150 (MQ=255)
                                                                                                                        ccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCGTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc                             >  1:1424561/1‑150 (MQ=255)
                                                                                                                         ctctCCCACCGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACt                            <  1:358529/150‑1 (MQ=255)
                                                                                                                          tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACtt                           >  1:2756747/1‑150 (MQ=255)
                                                                                                                             cccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTAt                        <  2:1491760/150‑1 (MQ=255)
                                                                                                                                acaGGGAGAGGGAGAAACCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACATTATCGa                     >  2:1648287/1‑150 (MQ=255)
                                                                                                                                     gagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGcc                <  2:1435724/150‑1 (MQ=255)
                                                                                                                                     gagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGcc                <  2:1435637/150‑1 (MQ=255)
                                                                                                                                      agagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCg               <  2:2483332/150‑1 (MQ=255)
                                                                                                                                       gtggtagaAACCACACCAGGCGTTCTCCTGGATATGGCATTGCTATTAGGTGCGGAACACCAAATTAATTATGTGATATCCGAAGTGGGTGTGCAGCGGGCCGGTGGGCACCTTTACCGTAAAATAGAAAACCACTTTATCGAACGCCgg              <  1:1954872/145‑1 (MQ=255)
                                                                                                                                       gtgggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCCGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCgg              <  1:1212996/148‑1 (MQ=255)
                                                                                                                                       gagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCATTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTCCCGGACAAGAGAAAACCACTTTATCGAACGCCgg              <  1:2880092/150‑1 (MQ=255)
                                                                                                                                         gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGaa            >  2:17506/1‑150 (MQ=255)
                                                                                                                                         gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGaa            <  2:2603487/150‑1 (MQ=255)
                                                                                                                                         gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGaa            <  2:284556/150‑1 (MQ=255)
                                                                                                                                          ggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAc           >  1:406631/1‑150 (MQ=255)
                                                                                                                                           gagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAcc          >  1:1077277/1‑150 (MQ=255)
                                                                                                                                            agaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCa         >  2:597315/1‑150 (MQ=255)
                                                                                                                                              aaaaCACTCACGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGTTACAGCACCTTAATGATGTGATATCCGAACTGGGGGTGCAGCGGGCCGTTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATc       <  1:259822/150‑1 (MQ=255)
                                                                                                                                              aaaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATc       >  1:1862713/1‑150 (MQ=255)
                                                                                                                                               aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCt      <  1:187774/150‑1 (MQ=255)
                                                                                                                                               aaaCACACAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGTTGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGAAAAGAGAAAACCACTTTATCGAACGCCGGAACCATCt      <  1:1811357/150‑1 (MQ=255)
                                                                                                                                                   acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACGAGTGCAAACCACTTTATCGGACGCCGGAACCATCTGAcc  >  1:2788102/1‑150 (MQ=255)
                                                                                                                                                   acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAcc  <  2:2661286/150‑1 (MQ=255)
                                                                                                                                                   acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAcc  >  1:806967/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACC  >  NC_000913/3959352‑3959648

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