Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,515 T→A 8.6% intergenic (+70/+17) 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,5150TA 8.6% 1159.0 / 53.6 409intergenic (+70/+17)ilvC/ppiCketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Reads supporting (aligned to +/- strand):  ref base T (143/197);  new base A (29/10);  total (199/210)
Fisher's exact test for biased strand distribution p-value = 1.44e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 416 total aligned reads displayed.

GCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCT  >  NC_000913/3959366‑3959664
                                                                                                                                                     |                                                                                                                                                     
   gCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                                                                                                                    <  2:2781499/150‑1 (MQ=255)
      aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAATCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCACAGGGAGAGGGGGAAAACACACAAGGCCGTCTCCgggggga                                                                                                                                                 >  1:454879/1‑146 (MQ=38)
      aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                                                 >  2:1426683/1‑150 (MQ=255)
      aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                                                 <  2:7604096/150‑1 (MQ=255)
      aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                                                 <  2:39707/150‑1 (MQ=255)
       tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                                                                                                                                <  1:3540747/150‑1 (MQ=255)
       tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCCACCCACTACCACAGGGCGAGGGAGAAAAAAAAAAAAGCCATCACCAGGAGAAg                                                                                                                                                >  1:4232358/1‑150 (MQ=255)
         cGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                                                                                              <  1:8195135/150‑1 (MQ=255)
         cGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                                                                                              <  1:3932519/150‑1 (MQ=255)
          gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCACCCACAACCACAGGGAGAGGGCGAAAAAACGAAAGGCCATATCCAGGAGAAGGcc                                                                                                                                             >  1:5577341/1‑150 (MQ=255)
            ttGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                                                                                                                                           <  2:3614659/150‑1 (MQ=255)
             tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                                                                                                                                          <  2:1737806/150‑1 (MQ=255)
              gAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTAACCCCCACCCGCAACCCCCGGGAGAGGGAGAAAACAAAAAAGGCCTTATCCTGGTGAAGGCCTTgg                                                                                                                                         >  1:4692317/1‑149 (MQ=255)
               aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCACTCCCACAGGGAGAGGGGGAAAAAACTCAATGCCATATCCTGGAGAAGGCCGTGCt                                                                                                                                        >  1:5061130/1‑150 (MQ=38)
                         taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGAACCCCGACCCTCTCCCACCCGGCGAGGGAGAAAACACACAATGCCGTATCCGGGAGAAGGCCGGGCGATTATTTGCg                                                                                                                              >  2:7065611/1‑150 (MQ=255)
                          aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                                                                                                                             <  2:3376926/150‑1 (MQ=255)
                             aaaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTACCACAGGGGGAGGGTGAAAAAACACAAGGCCGTCGCCAGGGGCAGGGCTGTGTATTAGTTGGGGtgc                                                                                                                          >  1:838152/1‑148 (MQ=25)
                                 tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCACCCACCGGGAGAGGGAGAAAACACTCAAGGCCGTCAACTGGAGAACGCCTTGCGAGTAGTTGCGGTACCGcc                                                                                                                      >  2:4538730/1‑149 (MQ=37)
                                 tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGGTACGTGCGCGCTGATGCCCGCACCCCGCCCCGCTCCCCCCGGGGGAGGGCGAAAACACGCAAGGCCGTCTCCCGGCGGAGGCCGGTCTATTTTTTGCGGGACAGcc                                                                                                                      >  2:110342/1‑149 (MQ=255)
                                   cgcgGCTATATGACAGATATGAACCGTATTGCTGTTGCGGGTTAAGTGCGCGCGGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTCGTTGCGGTACAGCAcc                                                                                                                    <  2:1000638/150‑1 (MQ=255)
                                   cgcgGCTATATGACAGATATGAAACGGATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCCCAGGGAGAGGGAGAACACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc                                                                                                                    <  2:4464812/150‑1 (MQ=255)
                                    gcgGCTATATGACGGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCGACCCTATACCACAGGGAGAGGGAGAAAAAAAAAAACGCCGTATCCTGGAGAAGGCCGTGCTATTATTTGCGGTAAAGCAcca                                                                                                                   >  1:6521507/1‑149 (MQ=25)
                                    gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCACCCAAAACCACAGGGCGAGGGGGAAAAAACACAACGCCTACTCCTGGGGGAGGGCGGGCTATTAGTTTGGGTACAGCaacg                                                                                                                   >  1:6835395/1‑147 (MQ=255)
                                      ggCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGGGAGGGCGAAAACACTCAAGGCCTTCTCCTGGAGTAGGCCTTGATATTAGTTGCGGGACCGCACCTTa                                                                                                                 >  2:7470607/1‑150 (MQ=38)
                                       gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa                                                                                                                <  1:3317679/150‑1 (MQ=255)
                                           tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgat                                                                                                            <  1:1618686/150‑1 (MQ=255)
                                            atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                                                                                           <  1:5525713/150‑1 (MQ=255)
                                             tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCACAACCCCACCCACTCCCCCAGGGAGAGGGAGAAAACACACAACGCCGAATCCAGGAGAAGGGCTTGCTAATAATTTCGGTAAAGCACCATaaatatgg                                                                                                          >  2:6731432/1‑144 (MQ=255)
                                              gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt                                                                                                                  <  1:5780597/141‑1 (MQ=255)
                                               aCAGAGATGAAACGTATTGCTGTTGCGGGTTTAGTGCGCGCTGAGGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGa                                                                                                        <  1:2130488/150‑1 (MQ=255)
                                                cAGATATGAACCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCACCCCGGCCCGCACCCACCGGGGGGGGGAGGAAAAACGCAACGCCGTGACCGGGCGCACGCCGTGCGAATAGGTGCGGGTCAGCCCCGTAatgatggggt                                                                                                       >  2:4024219/1‑146 (MQ=255)
                                                cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGat                                                                                                       <  2:6179444/150‑1 (MQ=255)
                                                cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCCCCCACAACCACCGGGGGGGGGAGAAAACACACAAGGCCGTCTCCGGGTGAAGGCCGGGCGATTAGTTTCGGTACAGCACCATAATGTtgtggt                                                                                                       >  2:8105113/1‑148 (MQ=255)
                                                 aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCAACCCGACCAACAACCACACGGAGAGGGAGAAAAAAAAAAAGGCCGTATCCAGGAGAAGCCCATGCTAATATTTGAGGTAAATCACCTTAATGATGTtata                                                                                                      >  2:2454262/1‑150 (MQ=255)
                                                 aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCACACCCCCAACCAAAACCACAGGGCGAGGGGGAAAACAAACAAGGCCATCACCTGGAGAAGGCCATTCTAATAGTTGCGGTAAAGCACCATAATGATTTGata                                                                                                      >  2:6017038/1‑150 (MQ=255)
                                                    tatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATcc                                                                                                   <  1:4715390/150‑1 (MQ=255)
                                                      tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGa                                                                                                 <  1:7864353/150‑1 (MQ=255)
                                                        aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATTCCCTACCCCCCACCCTCTCCAACAGGGAGAGCGAGAAAACACTCAAGGCCTTATCCTGGAGAAGGCCATGCTATTAGTTGAGGGACAGCACCATAATGAAGTGATATCCGAAc                                                                                               >  2:8596023/1‑150 (MQ=255)
                                                        aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                                                                                                                                      >  2:4901493/1‑95 (MQ=255)
                                                           cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                                            <  1:7067419/150‑1 (MQ=255)
                                                           cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                                            <  1:7713732/150‑1 (MQ=255)
                                                            gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCTCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTggg                                                                                           >  2:6462464/1‑150 (MQ=255)
                                                             tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCACGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgt                                                                                          <  1:7509026/150‑1 (MQ=255)
                                                             tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgt                                                                                          <  1:6218988/150‑1 (MQ=255)
                                                               gatcTGTGGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGGTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgtgt                                                                                        <  1:4179190/147‑1 (MQ=255)
                                                                    gTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGc                                                                                   <  2:5412758/150‑1 (MQ=255)
                                                                     tttcGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGTGGGAGAAAACACTCTAGGCCTTCTCCTGGAGAAGGCCTTGATATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCg                                                                                  <  1:4417415/147‑1 (MQ=255)
                                                                        cGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc                                                                               <  2:5090213/150‑1 (MQ=255)
                                                                          ggTTAAGTGCGCGCTGATGGCCACAACCCCACCCAAAACCACCCGGAGAGGGCGAAAAAACACAAAGCCTTCACCGGGGGAAAGCCGGTCGAATATTTGCGGTACAGAACCAAAAATATGTGTTAAGCGAACAGGGGGTGCAGCGGGccc                                                                             >  2:426407/1‑149 (MQ=255)
                                                                               aGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCgg                                                                        <  1:1011507/150‑1 (MQ=255)
                                                                                  gcgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGGGTGCAGCGGGCCGGTCGGCTc                                                                     <  1:6036165/150‑1 (MQ=255)
                                                                                  gcgcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCGTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTc                                                                     <  1:1730459/150‑1 (MQ=255)
                                                                                      gcTGATGCCCACACCCCCACCCACACCCACCGGGAGAGGGCGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCTTTCTAATAGTTTCGGTACACCACCATAATTATGTGATAAACGAACGGGGTGTGCCGCGGGCCGGTCGGCTCCAGt                                                                 >  2:8398980/1‑150 (MQ=255)
                                                                                          aTGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAACACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCg                                                             <  1:3676128/150‑1 (MQ=255)
                                                                                              ccTCCCCCCGACCCTCTCCCACAGGGAGCGGGAGAAAACACTCAAGGCCTTCGCCGGGAGAAGGCCGGGCGATTAGGTGCGGGACAGCACCTGCAGGATGGGATAGCCGAACGGGGTGTGCCGCGGGCCGGTCGGCTCCAGTACCGGACa                                                         >  2:7169796/1‑150 (MQ=255)
                                                                                                tCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACCCTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAag                                                       <  2:3584459/150‑1 (MQ=255)
                                                                                                  cccccGACCCTCTCCCTCTGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAagag                                                     <  1:520223/149‑1 (MQ=255)
                                                                                                  aCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAagag                                                     <  2:3706849/150‑1 (MQ=255)
                                                                                                    cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaa                                                   <  1:5316810/150‑1 (MQ=255)
                                                                                                    cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaa                                                   <  1:4470246/150‑1 (MQ=255)
                                                                                                         cccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCa                                              >  1:5980969/1‑150 (MQ=255)
                                                                                                          ccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc                                             >  1:3189753/1‑150 (MQ=255)
                                                                                                          ccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCAc                                             >  1:7290551/1‑150 (MQ=255)
                                                                                                                ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc                                       >  2:7875641/1‑150 (MQ=255)
                                                                                                                   caGGGAGAGGGAGAAAACCCTCAAGGCCTGCTCCTGGAGAAGGCCTTGCTATTCGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGTTCGGCTCAAGTACGGTACAAGAGAAAACCACTTTATCGac                                    <  1:283488/150‑2 (MQ=255)
                                                                                                                    agggagagggagAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTCCCGGACAAGAGAAAACCACTTTATCGAAc                                   >  2:4309395/1‑150 (MQ=255)
                                                                                                                     gggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACg                                  >  1:1899929/1‑150 (MQ=255)
                                                                                                                      ggagagggagaAAACACTCTAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATCTCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTATAGCGAACGc                                 >  2:5830422/1‑150 (MQ=255)
                                                                                                                      ggagagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTTGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCGGCGGGCCGGTCGGCTCCTGTACCGGACAAGAGAAAACCACTTTATCGAACGc                                 <  1:5117526/150‑1 (MQ=255)
                                                                                                                        agagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCg                               >  2:2213932/1‑150 (MQ=255)
                                                                                                                         gagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCgg                              >  2:6879076/1‑150 (MQ=255)
                                                                                                                         gagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCgg                              >  1:137822/1‑150 (MQ=255)
                                                                                                                           gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGaa                            <  1:6514858/150‑1 (MQ=255)
                                                                                                                            ggagaAAACACTCAAGGCCTTCTCCTGGAGACGGCCGTGCTATTAGTTGCGGTACAGCACCTTAAGGAGGGGATATCCGAACGGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTGTATCGAACGCCGGacc                           >  2:1943131/1‑148 (MQ=255)
                                                                                                                            ggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                                                                                                                             >  1:1497002/1‑52 (MQ=255)
                                                                                                                            ggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAc                           >  2:6282429/1‑150 (MQ=255)
                                                                                                                              agaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                                                                                               <  2:6577394/80‑1 (MQ=255)
                                                                                                                               gaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCAt                        >  2:2720531/1‑150 (MQ=255)
                                                                                                                               gaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCAt                        >  2:6820378/1‑150 (MQ=255)
                                                                                                                                  aaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACCGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTCCCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTg                     >  2:8329547/1‑150 (MQ=255)
                                                                                                                                  aaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTg                     >  1:2649287/1‑150 (MQ=255)
                                                                                                                                   acacTCACTGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGa                    <  2:7308348/150‑1 (MQ=255)
                                                                                                                                   acacTCACGGCCGTCGCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATTATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGAAAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGa                    <  1:2510394/150‑1 (MQ=255)
                                                                                                                                   acacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGa                    <  2:4215401/150‑1 (MQ=255)
                                                                                                                                   acacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGa                    <  2:2977403/150‑1 (MQ=255)
                                                                                                                                   acacTCAAGGCCGTCTCCTGGAGAAGGCATTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGa                    <  1:3714348/150‑1 (MQ=255)
                                                                                                                                         aaGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTgg              >  2:3807755/1‑150 (MQ=255)
                                                                                                                                         aaGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTgg              <  1:8666841/150‑1 (MQ=255)
                                                                                                                                          cggCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGc             <  1:6730239/149‑1 (MQ=255)
                                                                                                                                           ggCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCg            >  1:4867957/1‑150 (MQ=255)
                                                                                                                                              cTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGaa         <  1:426407/150‑1 (MQ=255)
                                                                                                                                              cTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGaa         <  1:789806/150‑1 (MQ=255)
                                                                                                                                               ttCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAAt        >  2:5695957/1‑150 (MQ=255)
                                                                                                                                               ttCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAAt        >  2:465825/1‑150 (MQ=255)
                                                                                                                                                  tcCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCa     >  1:6964123/1‑150 (MQ=255)
                                                                                                                                                   ccTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCAc    >  2:1551163/1‑150 (MQ=255)
                                                                                                                                                   ccTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCAc    <  2:555933/150‑1 (MQ=255)
                                                                                                                                                   ccTGGAGAAGGCCGAGCTATGAGTTGCGGTCCAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGGACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCAc    <  2:5125483/150‑1 (MQ=255)
                                                                                                                                                    cTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCAcc   >  2:2005965/1‑150 (MQ=255)
                                                                                                                                                    cTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCAcc   >  1:129353/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACCT  >  NC_000913/3959366‑3959664

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: