Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,515 T→A 6.9% 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 6.9% 259.2 / 7.4 96intergenic (+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 (34/54);  new base A (3/4);  total (38/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

AGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACC  >  NC_000913/3959368‑3959663
                                                                                                                                                   |                                                                                                                                                    
agggaTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCCATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                                                                                                                                     <  2:806907/145‑1 (MQ=255)
 gCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTgg                                                                                                                                                    <  1:1069993/149‑1 (MQ=255)
  ccATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                                                                                                                   <  1:272922/149‑1 (MQ=255)
  ccATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                                                                                                                   <  2:732126/149‑1 (MQ=255)
     tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                                                                <  1:1563951/149‑1 (MQ=255)
      gCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                                                                                                                               <  2:1421386/149‑1 (MQ=255)
      gCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                                                                                                                               <  2:1415217/149‑1 (MQ=255)
        gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                                                                                             <  1:1048913/149‑1 (MQ=255)
        gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGCCTCCTCCAGGAGAAGGc                                                                                                                                             <  1:587198/149‑1 (MQ=255)
            gAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                                                                                                                                         <  2:1567211/149‑1 (MQ=255)
            gAGCAGGTAGGTAAGAAACTGCGCGGATATATGACAGATATGAAACTTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                                                                                                                                         <  2:427782/149‑1 (MQ=255)
               cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                                                                                                                                      <  2:1568868/149‑1 (MQ=255)
                aggtaggtATGAAACTGCGCGGATATATGACATATATGAACCGGATTGCTGTTGCGGGTTAAGTGCGCGCGTCGGCCCTCACCCCGACACTCTCCAACAGGGAGAGGGAGAAAACACTCAAGGCGTTCTCCTGGAGAAGGCCTTGCTAt                                                                                                                                     <  2:90153/149‑1 (MQ=255)
                aggtaggtAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                                                                                                                                     <  1:647224/149‑1 (MQ=255)
                 ggtaggtcAGATCCTGCGCGGCTATATGACAGATATGACACGTATTGCTGTTGCGGGTGAAGTTCGCGCTGAGGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGACTTCTCCTGGAGAAGGCCTTgagtgt                                                                                                                                    <  1:234787/149‑6 (MQ=37)
                 ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTAACCCCGACCCACTCCCACACGGAGAGGGGGAAAACACACAAAGCCTTCTCCAGGAGAAGGCCGTGCTAtt                                                                                                                                    >  2:1321425/1‑149 (MQ=255)
                  gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGCCTTCTCCAGGAGAAGGCCTTgagtgtt                                                                                                                                   <  1:1562761/149‑7 (MQ=255)
                   taggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCTTgagtgttt                                                                                                                                  <  2:328907/149‑8 (MQ=255)
                     ggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                                                                                                                                <  1:35236/149‑1 (MQ=255)
                      gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                                                                                                                               <  1:1321425/149‑1 (MQ=255)
                        aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                                                                                                                             >  1:1114089/1‑149 (MQ=255)
                         aGAAACTGCGAGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTATCCAACAGGGAGAGGGAGAAAACACTCACGGCCATCTACAAGAGAAGGCCTTGCTATTAGTTgggg                                                                                                                            >  1:668090/1‑149 (MQ=37)
                          gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt                                                                                                                           >  1:932669/1‑149 (MQ=255)
                            aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTAc                                                                                                                         <  2:1574773/149‑1 (MQ=255)
                            aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTAc                                                                                                                         <  1:1163609/149‑1 (MQ=255)
                             aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa                                                                                                                        >  2:48574/1‑149 (MQ=255)
                               tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc                                                                                                                      >  2:1469045/1‑149 (MQ=255)
                                gcgcGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGGTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCa                                                                                                                     <  1:182097/149‑1 (MQ=255)
                                   cgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt                                                                                                                  <  2:1342195/149‑1 (MQ=255)
                                    ggCTATATGACATATATGAAACGTATTGCTGTTGGGGGTTTAGCGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAAGACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt                                                                                                                 <  2:684958/149‑1 (MQ=38)
                                      cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGACCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa                                                                                                               <  1:1000674/149‑1 (MQ=255)
                                      cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa                                                                                                               <  2:641799/149‑1 (MQ=255)
                                         tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatga                                                                                                            >  2:257831/1‑149 (MQ=255)
                                          atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgat                                                                                                           <  2:1114089/149‑1 (MQ=255)
                                           tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                                                                                          <  2:1291376/149‑1 (MQ=255)
                                            gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGGGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACTGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                                                                                         <  1:315593/149‑1 (MQ=255)
                                            gACAGATATGAAACGGATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                                                                                         <  2:660390/149‑1 (MQ=255)
                                               aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTGGCTATTAGTTGCGGTACAGCACCGTAATGATGTGat                                                                                                      >  2:992099/1‑149 (MQ=255)
                                               aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTACCACAGGGAGAGGGAGCAAACACTCAAGGCCTTCTCCTGGTGAAGGCCGTGCGATTAGTTGCGGTACAGCACCTTAATGATGTGat                                                                                                      >  1:1511865/1‑149 (MQ=38)
                                                 atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACACCGACCCTCTCCCCCCGGGAGAGGGAGAAAACACACAAGGCCTACACCTGGAGAAGGCCTAGCTAATAGGTGCGGTACAGCACCTTAATGATGTGataa                                                                                                    >  2:74745/1‑148 (MQ=255)
                                                  tatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCCCCTTAATGATGTGctata                                                                                                   >  2:146846/1‑148 (MQ=255)
                                                  tatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATc                                                                                                   <  1:1133950/149‑1 (MQ=255)
                                                  tatGAAACGTATTGCTGTTGCGGGGTAAGGGCGCGCGGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATc                                                                                                   <  2:111681/149‑1 (MQ=255)
                                                     gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGa                                                                                                <  2:1110526/149‑1 (MQ=255)
                                                     gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCGTCGCCGGTCGCAGGCCGTGCTATTAGTGGCGGGACAGCACCTTAATTCCGTGATATACGa                                                                                                >  2:784605/1‑149 (MQ=25)
                                                      aaaCGTATTGCTGTGGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGaa                                                                                               <  2:301666/149‑1 (MQ=255)
                                                       aaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGACCTCACCCCGACCCTCTCCCACAGTGAGAGGGAGATCACACTCAAGGCCTTCGCATGGAGAAGGCCTTTCTATTAGGTGCGGTACATCACCTTAATGATGTGATATGCGAAc                                                                                              >  1:266822/1‑149 (MQ=255)
                                                          gTATTGCTGTTGCGGGTTATGAGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                                           <  2:558434/149‑1 (MQ=255)
                                                          gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                                           <  1:668488/149‑1 (MQ=255)
                                                          gTATTGCTGTTGCGGGATAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTgg                                                                                           <  1:344354/149‑1 (MQ=255)
                                                           tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTggg                                                                                          >  2:19849/1‑149 (MQ=255)
                                                            aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTATTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGgt                                                                                         >  2:1138283/1‑149 (MQ=255)
                                                                cTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGc                                                                                     >  2:276164/1‑149 (MQ=255)
                                                                      cGGGTTAAGTGCGCGCTGATGCCCTCCCCCCGCCCCTCTCCCACAGGGAGAGGGAGAAAACACTCACGGCCTTCTCCTGGAGAAGGCCTGGCTATTAGTTTCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCggg                                                                               <  2:235941/149‑1 (MQ=38)
                                                                       gggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc                                                                              <  1:950605/149‑1 (MQ=255)
                                                                        ggTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCACAGGGAGAGGGAGAAAACACACAAGGCCTTATCCTGGAGAAGGCCTTGCTAATAGTTTCGGTACAGCACCTGAATGATGTGATATACGACCTGGGTGTGCAGCGGGcc                                                                             >  2:327950/1‑149 (MQ=37)
                                                                        ggTTAAGTGCGCGCTGATGCCCACACCCCGACCCTCTACCACACGGAGAGGGAGAAAACACTCAAAGCCTTCTCCTGGAGAAGGCCTTGGTATTAGTTTCGGTACAGCACCTTAATTATGATAAATCCGAACTGGGTGTGACGCGGGcc                                                                             >  2:733191/1‑149 (MQ=255)
                                                                         gTTAAGTGCGCGCTGATGCCCTCCCCCCGACCCAATCCCACAGGGAGAGGGAGAAAACCCTCAAGGCCGTCTCCTGGAGAAGGCCGTGCTATTAGGTTCGGTACAGCACCGGAATGATGTGATATCCTCACTGGGTGTGGAGCGGGCCg                                                                            >  1:301322/1‑149 (MQ=25)
                                                                         gTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCg                                                                            <  1:733191/149‑1 (MQ=255)
                                                                           tAAGTGCGCGCTGATGCCCTGACCCCGTCCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGt                                                                          <  1:1216694/149‑1 (MQ=255)
                                                                              gTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCgg                                                                       <  2:1151682/149‑1 (MQ=255)
                                                                                     cTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGt                                                                <  2:1531167/149‑1 (MQ=255)
                                                                                        aTGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTTCAGCGGGCCGGTCGGCTCCAGTAcc                                                             <  2:295075/149‑1 (MQ=255)
                                                                                         tGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCg                                                            <  1:257831/149‑1 (MQ=255)
                                                                                             cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACa                                                        >  2:1279069/1‑149 (MQ=255)
                                                                                                  cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgaga                                                   >  1:1430616/1‑149 (MQ=255)
                                                                                                  cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgaga                                                   >  2:840889/1‑149 (MQ=255)
                                                                                                  cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgaga                                                   >  2:809094/1‑149 (MQ=255)
                                                                                                     gACCCTCTCCCCCAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTTATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaaa                                                <  2:668090/149‑1 (MQ=255)
                                                                                                      tcccTCTCCCTGTGGGAGAGGGTGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc                                               >  2:702389/2‑149 (MQ=255)
                                                                                                      tcccTCTCCCTGTGGGAGAGGGTGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCGTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc                                               >  1:263884/2‑149 (MQ=255)
                                                                                                      aCCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAc                                               <  2:830219/149‑1 (MQ=255)
                                                                                                        ccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAAccc                                             >  1:770561/1‑148 (MQ=255)
                                                                                                          tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACt                                           >  2:412368/1‑149 (MQ=255)
                                                                                                          tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACt                                           >  1:1139525/1‑149 (MQ=255)
                                                                                                              ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCTCCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGCCAAGAGAAAACCACTTTAt                                       >  1:651277/1‑149 (MQ=255)
                                                                                                               cacctggagagggagaAAACACTCAAGGCCATCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGAAATCCGAACTGGGTGTTCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc                                      <  2:1038990/144‑1 (MQ=255)
                                                                                                               cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc                                      <  2:982272/149‑1 (MQ=255)
                                                                                                               cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc                                      <  1:856866/149‑1 (MQ=255)
                                                                                                                 caGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGa                                    >  1:762199/1‑149 (MQ=255)
                                                                                                                         gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCTGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGa                            >  1:25608/1‑149 (MQ=255)
                                                                                                                           gagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGCGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAc                          <  1:48574/149‑1 (MQ=255)
                                                                                                                              aaaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCAt                       <  2:1430616/149‑1 (MQ=255)
                                                                                                                                  cacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGa                   <  1:1279069/149‑1 (MQ=255)
                                                                                                                                   acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGCACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAGCCCCTTTATCGAGCGCCGGAACCATCTGAc                  >  1:774615/1‑149 (MQ=255)
                                                                                                                                   acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAc                  <  1:412368/149‑1 (MQ=255)
                                                                                                                                   acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAc                  >  1:546733/1‑149 (MQ=255)
                                                                                                                                    cTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAcc                 <  1:146846/149‑1 (MQ=255)
                                                                                                                                       aaGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTg              >  2:772380/1‑149 (MQ=255)
                                                                                                                                         ggCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGc            >  2:1264139/1‑149 (MQ=255)
                                                                                                                                             ttCTCCTGGAGAAGGCCTTGCTCTTAGTTGCGGTACAGCACCTAAATGACGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGaa        <  1:784605/149‑1 (MQ=255)
                                                                                                                                             ttCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGGGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGaa        <  2:301322/149‑1 (MQ=255)
                                                                                                                                                  cTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGACCCATCTGACCCTGGCGGAATTCAc   >  1:108773/1‑149 (MQ=255)
                                                                                                                                                  cTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCAc   >  1:733125/1‑149 (MQ=255)
                                                                                                                                                  cTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCAc   >  2:955350/1‑149 (MQ=255)
                                                                                                                                                  cTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCAc   <  2:546733/149‑1 (MQ=255)
                                                                                                                                                   tGGAGAAGGCCTTGCTATTAGTTGCGTTACAGCTCCTTAATGATGTGATATCCGAACAGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACGTCTTTATCGAACGCCGGAAGCATCTGACCCTGGCGGAATTCAcc  <  1:74745/149‑1 (MQ=255)
                                                                                                                                                   tGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCAcc  >  2:740441/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAATTCACC  >  NC_000913/3959368‑3959663

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: