Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,480 A→T 5.8% intergenic (+35/+52) 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,4800AT 5.8% 232.3 / 4.4 98intergenic (+35/+52)ilvC/ppiCketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Reads supporting (aligned to +/- strand):  ref base A (34/54);  new base T (3/3);  total (41/57)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

AACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTAT  >  NC_000913/3959332‑3959626
                                                                                                                                                    |                                                                                                                                                   
aaCGGGCAACTGCGTGATGTGAACGAAGCGATTCGCATCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGTCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCTCGCTGATGCCCTCACCCCGACCCTCTCCcac                                                                                                                                                    <  1:1174713/150‑1 (MQ=255)
 aCGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCCCCGACCCTCTACcaca                                                                                                                                                   >  2:746245/1‑150 (MQ=255)
  cGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACag                                                                                                                                                  <  1:906827/150‑1 (MQ=255)
  cGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACag                                                                                                                                                  >  1:1105189/1‑150 (MQ=255)
    ggCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACaggg                                                                                                                                                <  1:1204976/150‑1 (MQ=255)
    ggCAACTGCGGGATGTGAACGCAGCGAATCGCAGCCATGCGATTGAGCAGGTAGGGAAGAAACTGCGCGGCTATATGACAGATCTGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACaggg                                                                                                                                                <  1:1129496/150‑1 (MQ=39)
         cTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCTGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTATGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagagg                                                                                                                                           <  2:62141/150‑1 (MQ=255)
          tGCGTGATGTGAACGACGCGATTAGCATCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagaggg                                                                                                                                          <  2:335328/150‑1 (MQ=255)
          tGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGGGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagaggg                                                                                                                                          <  2:676111/150‑1 (MQ=255)
          tGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagaggg                                                                                                                                          <  2:73524/150‑1 (MQ=255)
          tGCGGGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagaggg                                                                                                                                          <  2:1042565/150‑1 (MQ=255)
           gCGGGATGTGAACGAAGCGATTCCCAGCCATGCGATTGAGCAGGTAGGTACGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagaggga                                                                                                                                         <  2:1681617/150‑1 (MQ=255)
              tGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGaa                                                                                                                                      >  1:951573/1‑150 (MQ=255)
               gATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGaaa                                                                                                                                     <  2:1105189/150‑1 (MQ=255)
                aTGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGaaaa                                                                                                                                    >  2:546875/1‑150 (MQ=255)
                  gtgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAaca                                                                                                                                  <  1:1174888/150‑1 (MQ=255)
                  gtgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAaca                                                                                                                                  >  1:1250244/1‑150 (MQ=255)
                     aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTc                                                                                                                               <  1:713732/150‑1 (MQ=255)
                       cGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                                             <  1:911162/150‑1 (MQ=255)
                           gCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGcc                                                                                                                         >  1:563936/1‑150 (MQ=255)
                           gCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGcc                                                                                                                         <  2:733729/150‑1 (MQ=255)
                            cGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCt                                                                                                                        <  2:1113516/150‑1 (MQ=255)
                            cGAGTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCt                                                                                                                        <  2:1569230/150‑1 (MQ=255)
                              aTTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGGTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCCCAGGGAGAGGGAGAAAACACACAAGGCCttg                                                                                                                      >  2:1092709/1‑149 (MQ=255)
                               ttCGCAGCCATGCGATTTAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAATACACTCAAGGCCTtct                                                                                                                     <  2:77727/150‑1 (MQ=255)
                                  gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                                                                                  >  2:1466138/1‑150 (MQ=255)
                                   cAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                                                                                                 >  2:1166538/1‑150 (MQ=255)
                                     gCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                                                                               >  1:1182865/1‑150 (MQ=255)
                                      ccATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgag                                                                                                              <  1:337993/150‑1 (MQ=255)
                                        aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                            >  2:1031707/1‑150 (MQ=255)
                                          gataTTGATCATGTAGGTAAGATTATGCGCGGCTATATGTCAGATATTAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACACTCTCCGACAGGGAGTGGGAGAATTCAGTCAATGACTTTTCGTGGAGATgg                                                                                                          <  2:1644876/147‑1 (MQ=21)
                                           cGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCCCCTGGAGAAGGc                                                                                                         >  2:322596/1‑150 (MQ=255)
                                           agATTGAGCAGGTAGGTAATAAAATGCGCTTCTATATTACAGATATGAAAAGTATTGCTGTTGCTGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGc                                                                                                         <  2:1450199/149‑1 (MQ=37)
                                            taTTGAGCAGGTAGGAAAGAAACTGGGCGGCTATATGACAGATATTAAACGTATTGCTGTTGCGGGTAAAGTGCGCGCGGATGCCCGCACCCCGACCGTCTCCCACAGGGTGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGcc                                                                                                        <  1:877345/149‑1 (MQ=255)
                                            gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATAGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                                                                                        <  1:1023892/150‑1 (MQ=255)
                                              ttGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                                                                                                      <  1:841378/150‑1 (MQ=255)
                                               tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTGAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGATAACACTCAAGGCCTTCTCCAGGAGAAGGCCTTg                                                                                                     <  2:822772/150‑1 (MQ=255)
                                                gAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                                                                                                    <  1:1425792/150‑1 (MQ=255)
                                                gAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                                                                                                    <  1:2777/150‑1 (MQ=255)
                                                gAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                                                                                                    <  1:9410/150‑1 (MQ=255)
                                                gAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTAGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGATAAGGACTTGc                                                                                                    >  2:317759/1‑150 (MQ=255)
                                                 aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                                                                                                   <  2:1455640/150‑1 (MQ=255)
                                                 aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                                                                                                   <  1:966954/150‑1 (MQ=255)
                                                  gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                                                                                                  <  2:1673527/150‑1 (MQ=255)
                                                  gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                                                                                                  <  2:206918/150‑1 (MQ=255)
                                                     ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                                                                                               <  1:1466138/150‑1 (MQ=255)
                                                         ggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCCCCCACCCAATACCACACGGAGACGGAGAAAACACACAATGCCTTGAACAGGATAAGGACATTCTATTAGttt                                                                                           >  2:886450/1‑149 (MQ=255)
                                                          gtaagtaaCTGCGCGGGTATATGACAGCTATGAAACGTATTGCTGTTGCGGGTAAAGTGCGCGCTGATGCCCTCACCCCGCCCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                                                                                          <  1:239745/150‑1 (MQ=38)
                                                          gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACACAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                                                                                          >  2:465891/1‑150 (MQ=255)
                                                           taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                                                                                         <  2:957135/150‑1 (MQ=255)
                                                            aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGCCTTCTCCTGGAGAAGGCCTTGGTATTAGTTGCgg                                                                                        >  2:1716075/1‑150 (MQ=255)
                                                             aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt                                                                                       <  1:1632470/150‑1 (MQ=255)
                                                                  cTGCGCGGCTATATGCCAGATATGAAACGAATAGCTGTTGCGGGTTCAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAACGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc                                                                                  <  1:1088547/150‑1 (MQ=38)
                                                                  cTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc                                                                                  <  1:123719/150‑1 (MQ=255)
                                                                  cTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCGACCCACTCCCACAGGGAGAGGGAGAAAAAACACAACGCCTATTCCTGGAGAAGGCCATGGTATTAGTTTCGGGAAAGc                                                                                  >  2:184320/1‑150 (MQ=37)
                                                                     cgcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc                                                                               <  1:1111103/150‑1 (MQ=255)
                                                                              atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATTCCCTCACCCCCACCCTCACCCACAGGGAGAGGGAGAAAACACTCAAGTCCCTGAACTTGGGAAAGCCGTGCTATTAATTGCGGGACACAACC‑ATAattaag                                                                     >  2:1265300/1‑146 (MQ=21)
                                                                              atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCACCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAatgatg                                                                     <  1:1069926/150‑1 (MQ=255)
                                                                               tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGTGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGAtgt                                                                    >  1:745764/1‑150 (MQ=255)
                                                                               tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGAtgt                                                                    >  1:846781/1‑150 (MQ=255)
                                                                                gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTTTTAGTTGCGGTGCAGCACC‑TTAATGAtgtg                                                                   >  2:163702/1‑150 (MQ=255)
                                                                                  cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACACCGACCCTATACCACAGGGAGAGGGAGAAAACACACAAGGCCTTCAAATGGAGAAGGCCATGCAATTAGTTTCGGTAACGAACC‑TTAATTATGTGat                                                                 >  1:838654/1‑150 (MQ=25)
                                                                                    gATATGAAACGTATTGCTGGTGCGGGTTACGTGCGCGCTGATGCCCTAACACCGACCCTCTCCCTCAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGatat                                                               <  2:878545/150‑1 (MQ=255)
                                                                                      tatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTAACCCCGACCCAATCCCACAGGGAGAGGGGGAAAACACTAAAGGCCGTCTCCTGGAGGAGGCCTTGCTATTTTTTTCGGTACAGCACC‑ATAATTATTTGatatac                                                             >  2:1481097/1‑148 (MQ=255)
                                                                                       atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCCCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCg                                                            <  1:322596/150‑1 (MQ=255)
                                                                                           aaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAAcg                                                        >  2:252536/1‑149 (MQ=255)
                                                                                               tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGgt                                                    >  2:37997/1‑150 (MQ=255)
                                                                                                aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGgtg                                                   <  2:563936/150‑1 (MQ=255)
                                                                                                  tGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGtgtg                                                 <  1:886450/150‑1 (MQ=255)
                                                                                                    cTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGct                                               <  1:546875/150‑2 (MQ=255)
                                                                                                    cTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCACGACCCTATCCCACATGGAGAGGGAGAAAAAACTCAAGGCCATAAACTGGAGAAGGCCTTGCTATTATTTGCGGTACAGCCCC‑TTAATGATGTGATATACGAACTGGGTATGCa                                               >  1:836716/1‑150 (MQ=37)
                                                                                                     tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAg                                              <  1:662738/150‑1 (MQ=255)
                                                                                                       ttGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCg                                            <  2:989889/150‑1 (MQ=255)
                                                                                                          cGGGTTAAGTGCGCGCTGATGCCCTCCCCCCGACCCACACCCCCCCGGGTAGGGTGAAAACACACAACGGCATATCCAGGGGAAGGGCGTTGGATTATTTTGCGTGCACCACCGTTAATGA‑GTGGTATACGCACAGGGGGTTCAGCGGGc                                         >  2:1257795/1‑150 (MQ=255)
                                                                                                             gTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCgg                                      >  1:115228/1‑150 (MQ=255)
                                                                                                             gTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCgg                                      <  2:951573/150‑1 (MQ=255)
                                                                                                                 aGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCgg                                  >  1:545/1‑150 (MQ=255)
                                                                                                                   tGCGCGCTGATGCCCTCACCCCGACCCACTCCCCCAGGGAGAGGGAGAAAACACACAACGCCTTCTCCGGGAGAAGGGCGTGCTTTTAGTTGCGGTACAGCACC‑GTAATGATGGGTTATACGAACTGGGTGTGCCGCGGGCCGGTCGGcg                                >  2:1487180/1‑149 (MQ=25)
                                                                                                                   ggcgcgcTGACGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCCTCGCCGGGAGAAGGCCTTGCTATTAGTCGCGGGACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCt                                <  2:80113/149‑1 (MQ=38)
                                                                                                                     cgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTcc                              <  2:1715074/150‑1 (MQ=255)
                                                                                                                     cgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTcc                              <  1:566243/150‑1 (MQ=255)
                                                                                                                     cgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTcc                              <  2:1530934/150‑1 (MQ=255)
                                                                                                                          tGATGCCCTCACCCCGACCCTCACCCACAGGGAGAGGGAGAAAACACACAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGAGGTGATATCCGAACTGGGTGTGCACGGGGCCGGACGGCTCCAGTAc                         >  2:648657/1‑150 (MQ=39)
                                                                                                                           gATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAGACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAcc                        >  2:112138/1‑150 (MQ=255)
                                                                                                                            aTGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCg                       <  1:1703195/150‑1 (MQ=255)
                                                                                                                                 cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                  >  1:602693/1‑150 (MQ=255)
                                                                                                                                    aCCCCCACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAagag               <  1:676215/150‑1 (MQ=255)
                                                                                                                                     ttttctcccTCTCCCTGTGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgaga              >  1:1030035/7‑150 (MQ=255)
                                                                                                                                     ccccGACCCTCTCCCCCAGGGAGAGGGAGAAAACACTCAAGGCCTTCACCTGCAGAAGGCCGTGCGATTAGTTGCGGGCCAGCACC‑GTAATGAAGAGATATCCGATCGGGGTGTGCAGCGGGCCGGTCGGCTGCAGTACCGGACAAgaga              >  1:240508/1‑150 (MQ=255)
                                                                                                                                     ccccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgaga              <  2:1667140/150‑1 (MQ=255)
                                                                                                                                      tttctcccTCTCCCTGTGGGAGAGGGTGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaa             >  1:641423/6‑150 (MQ=255)
                                                                                                                                       ttctcccTCTCCCTGTGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaa            >  1:1335224/5‑150 (MQ=255)
                                                                                                                                         gACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc          <  2:998524/150‑1 (MQ=255)
                                                                                                                                           cccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGgagaagagaacgcca        >  1:256846/1‑145 (MQ=255)
                                                                                                                                           cccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCa        >  1:1173142/1‑150 (MQ=255)
                                                                                                                                            ccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc       <  1:927993/150‑1 (MQ=255)
                                                                                                                                             ctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTCCCCCCCC‑TTAATGATGTGATTTCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACt      >  2:1153272/1‑150 (MQ=255)
                                                                                                                                             ctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATCTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACt      >  1:1384212/1‑150 (MQ=255)
                                                                                                                                               ctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACttt    <  2:1182865/150‑1 (MQ=255)
                                                                                                                                                 cccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCACCAGTACCGGACAAGAGAAAACCACTTTAt  >  1:938063/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACC‑TTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTAT  >  NC_000913/3959332‑3959626

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑