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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,959,5010TA10.0% 1247.1 / 29.7 416intergenic (+56/+31)ilvC/ppiCketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Reads supporting (aligned to +/- strand):  ref base T (134/221);  new base A (47/2);  total (190/226)
Fisher's exact test for biased strand distribution p-value = 4.39e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 421 total aligned reads displayed.

AACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCT  >  NC_000913/3959353‑3959650
                                                                                                                                                    |                                                                                                                                                     
aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTc                                                                                                                                                      >  2:4617495/1‑150 (MQ=255)
aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCAATCCCACCGGGAGAGGGAGAAAACACTc                                                                                                                                                      >  1:10955025/1‑150 (MQ=255)
  cGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                                                                    <  1:350707/150‑1 (MQ=255)
  cGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                                                                    <  1:1444847/150‑1 (MQ=255)
   gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAGACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                                                                                                                   >  2:7337611/1‑150 (MQ=255)
   gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACAACCACACGGAGAGGGAGAAAAAACACaaa                                                                                                                                                   >  1:4147323/1‑149 (MQ=255)
   gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCACCCCGACCCTCACCCACAGGGGGTGGGGGTAAACACTCAAg                                                                                                                                                   >  1:1815031/1‑150 (MQ=255)
    aaGCGATTCGCAGCCATGCTTTTGAGCAGGTAGGTAAGAATCTGCGCGGGTATATGACAGATATGAAACGTATTGCTGTTTCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                                                                                                                  <  2:982913/150‑1 (MQ=38)
    aaGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACTCCCACCGGGAGAGGGAGAAAAAACGCaacg                                                                                                                                                  >  1:296555/1‑148 (MQ=255)
    aaGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCACCCCGCCCCGCGCCCCCCGGGCGCGGGAGGAAACACACaacg                                                                                                                                                  >  1:5759229/1‑148 (MQ=255)
     aGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGc                                                                                                                                                 <  1:3296807/150‑1 (MQ=255)
     aGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGc                                                                                                                                                 <  1:9862665/150‑1 (MQ=255)
       cGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCACAGGGAGAGGGAGAAAACACACAAGGCCt                                                                                                                                               >  1:2721578/1‑150 (MQ=255)
         aTTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGAGATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCACGGCCTtc                                                                                                                                             <  2:3962385/150‑1 (MQ=255)
           tCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATTCCCACACCCCCACCCACACCCCCAGGGAGAGGGAGAAAAAACAAAAGGccatatc                                                                                                                                           >  2:4977971/1‑145 (MQ=255)
            cGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGGGAGGGAGAAAACACTCAAGGCCTTCTcc                                                                                                                                          >  1:2224649/1‑150 (MQ=255)
             gCAGCCCTGCGAGTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATAGGAAACGTATTGCTGTTGCGGGTGAAGTGCGCGCGGCGGCCCTCACCCCGACCCTCTCCCACAGGGAGCGGGAGAAAACACGCACGGCCTTCGccg                                                                                                                                         >  2:6284166/1‑149 (MQ=255)
             gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                                                                                                         <  2:7937526/150‑1 (MQ=255)
               aGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCACACCCACAGGGAGAGGGAGAAAAAACACAAAGCCTTCTCCTgg                                                                                                                                       >  2:7649310/1‑150 (MQ=39)
                gCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                                                                                                                      <  1:9664625/150‑1 (MQ=255)
                     gCGATTGAGCATGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGGGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                                                                                                                 <  1:6840184/150‑1 (MQ=255)
                        aTTGAGCAGGTAGGTAAGAATCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                                                                              <  1:5255425/150‑1 (MQ=255)
                        aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                                                                              <  1:10719899/150‑1 (MQ=255)
                          tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACCCTCCAGGCCTTCTCCTGGAGAAGGCCTTg                                                                                                                            <  2:10591578/150‑1 (MQ=255)
                           gAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                                                                                                                           <  2:2748378/150‑1 (MQ=255)
                           gAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATACGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACGACCACAGGGCGAGGGAGAAAAAACAAAAAGCCATCTCCTGGAGAAGGCCATGc                                                                                                                           >  2:307135/1‑150 (MQ=255)
                             gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                                                                                                                         <  2:10263993/150‑1 (MQ=255)
                                    ggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                                                                                                                  <  2:1358092/150‑1 (MQ=255)
                                    ggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                                                                                                                  >  1:9075489/1‑150 (MQ=255)
                                     gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCACCCACGCCCACACGGGGTGGGCGAAAACACACAAGGCCATCTACGGGGGAAGGCCGTGCGATTAGTTGc                                                                                                                 >  2:2330812/1‑150 (MQ=255)
                                     gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACTCCCACAGGGAGTGGGAGAAAACACACAAGGCCATAACCAGGAGGAGGCCGTGCTATTAGtttc                                                                                                                 >  1:10882101/1‑148 (MQ=255)
                                       taGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTGAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCACACAGGGAGAGTGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                                                                                                               <  2:3758978/149‑1 (MQ=255)
                                       aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTGACGTGCGCGCTGAGGCCCGCACCCCGTCCCGCTCCCCCCGGGCGGGGGCGAAACCACGCAAGGCCCTCGCCGGGGGGCGGGCGTGGGATTAGGTGCgg                                                                                                               >  2:3443992/1‑150 (MQ=255)
                                              tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCa                                                                                                        <  2:9216335/150‑1 (MQ=255)
                                              tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCa                                                                                                        <  1:7671123/150‑1 (MQ=255)
                                              tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACACCCCCAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCATGCTATTAGTTGCGGTACAGCa                                                                                                        >  2:5307019/1‑150 (MQ=39)
                                              tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACTACCAAAGGGGGAGGGGGAAAACAAACAAGGCCAAGTCCGGGAGAAGGCCTGGCTATTAGTTGGGGGACAGCa                                                                                                        >  2:9577177/1‑150 (MQ=255)
                                                cgcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc                                                                                                      <  1:2499413/150‑1 (MQ=255)
                                                cgcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc                                                                                                      <  2:10706610/150‑1 (MQ=255)
                                                cgcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc                                                                                                      <  2:9385348/150‑1 (MQ=255)
                                                     cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                                                                                 <  1:10018185/150‑1 (MQ=255)
                                                       atatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatga                                                                                               <  1:3326912/150‑1 (MQ=255)
                                                         atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCACCCAAACCCCCAGGGCGAGGGAGAAAACAAACAAGGCCATCAACAGGAGAAGGCCGTGCTATTATTTGCGGTACACCACCAAAatgatg                                                                                             >  2:8646170/1‑150 (MQ=255)
                                                          tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                                                                            <  1:5094421/150‑1 (MQ=255)
                                                          tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                                                                            <  1:3184245/150‑1 (MQ=255)
                                                           gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgtg                                                                                           <  1:6347263/150‑1 (MQ=255)
                                                            aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGa                                                                                          <  1:2386625/150‑1 (MQ=255)
                                                             cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGat                                                                                         <  1:398226/150‑1 (MQ=255)
                                                              aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCACCCGCAACCACACGGGGCGGGGGAAAAAACAAAAGGCCATAGCCGGGGGGAGGCCGTGCGATTAGTTGCGGTACAGCCCCATAATGGTGTGtta                                                                                        >  2:2268250/1‑150 (MQ=255)
                                                                atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATc                                                                                      <  2:4106597/150‑1 (MQ=255)
                                                                           aTTGCTGTTGCGGGTTAAGTGCGCGCTGATTCCCTCACCCCGACCCACTCCCACAGGGAGAGGGAGAAAAAACTCAAGGCCTTATCCTGGAGAAGGCCTTGGTATTAGTTGAGGTACAGCACCATAAAGAATTGATATCCGAACTGGgtg                                                                           >  2:1638317/1‑150 (MQ=37)
                                                                                  ttGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCg                                                                    <  1:4467332/150‑1 (MQ=255)
                                                                                   tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCgg                                                                   <  1:8783102/150‑1 (MQ=255)
                                                                                    gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCggg                                                                  <  2:6771778/150‑1 (MQ=255)
                                                                                     cGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc                                                                 <  1:10497069/150‑1 (MQ=255)
                                                                                     cGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc                                                                 <  1:6018145/150‑1 (MQ=255)
                                                                                     cGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGc                                                                 <  1:6090253/150‑1 (MQ=255)
                                                                                       ggTTAAGTGCGCGCTGATGCCCTCACCCCCAACCAAAACCACACGGAGAGGGAGAAAAAACAAAAAGCCGTGTACTGGTGGAGGCCCTGCGATTATTTGCGGGTCAGCACCATAAATATGTGATATCCCAACAGGGTGTTAAGCGGGCCg                                                               >  2:4043718/1‑150 (MQ=255)
                                                                                         ttAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGGGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGt                                                             <  2:6129759/150‑1 (MQ=255)
                                                                                         ttAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGt                                                             <  2:6435276/150‑1 (MQ=255)
                                                                                         ttAAGTGCGCGCTGATGCCCTCACCCCCACCCACAACCACAGGGAGAGGGAGAAAAAACAAAAGGCCTTAAACAGGAGAAGGCCATTCTATTAGTTTCGGTACAGCACCATAATGATGTGATAAACGAACAGGGTGTGCAGCGGGCCGGt                                                             >  2:697630/1‑150 (MQ=25)
                                                                                         ttAAGTGCGCGCTGATGCCCTAACCCCGACCCACTACCACCGGGAGAGGGAGAAAACACACAATGCCTTCTCCTGGAGAAGGCCTTGCGATTAGTTGCGGTACAGCACCTTAATGATGTTATATACGAAATGGGTGTGCAGCGGGCCggg                                                             >  2:4357634/1‑149 (MQ=255)
                                                                                           aaGTGCGCGCTGATGCCCTCACCCCGACCCTCTACCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTATTTGCGGTACAGCACCTTAATGATGTGATATCCGAACAGGGTGTGCAGCGGGCCGGTCg                                                           >  2:10056293/1‑150 (MQ=255)
                                                                                                 gcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa                                                     <  2:10955025/150‑1 (MQ=255)
                                                                                                 gcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCGTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa                                                     <  1:7031875/150‑1 (MQ=255)
                                                                                                      gATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAcc                                                <  1:1626694/150‑1 (MQ=255)
                                                                                                      gATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAcc                                                <  1:7338286/150‑1 (MQ=255)
                                                                                                           ccTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACa                                           >  2:2550786/1‑150 (MQ=255)
                                                                                                            cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACac                                          >  2:7116350/1‑149 (MQ=255)
                                                                                                            cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                                          <  1:8145310/150‑1 (MQ=255)
                                                                                                            cTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACaa                                          >  2:5972110/1‑150 (MQ=255)
                                                                                                                 cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaa                                     <  1:6600177/150‑1 (MQ=255)
                                                                                                                  ccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaa                                    >  2:4573538/1‑150 (MQ=255)
                                                                                                                   cGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGaaaa                                   <  1:3892053/150‑1 (MQ=255)
                                                                                                                     tcccTCTCCCTGTGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                                                                                                                   <  2:2095106/67‑1 (MQ=38)
                                                                                                                      cccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCa                                <  2:5638917/150‑1 (MQ=255)
                                                                                                                        ctctCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGGACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACACGAGAAAACCACt                              >  2:8856791/1‑150 (MQ=255)
                                                                                                                         tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGACAACCACtt                             >  2:8160380/1‑150 (MQ=255)
                                                                                                                         tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACtt                             >  1:7961347/1‑150 (MQ=255)
                                                                                                                         tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACtt                             >  2:7461606/1‑150 (MQ=255)
                                                                                                                           tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTa                           >  2:7772800/1‑150 (MQ=255)
                                                                                                                           tcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTa                           >  1:4681425/1‑150 (MQ=255)
                                                                                                                             ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc                         >  1:9895080/1‑150 (MQ=255)
                                                                                                                             ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc                         >  1:3842610/1‑150 (MQ=255)
                                                                                                                             ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTACAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc                         >  1:8611989/1‑150 (MQ=255)
                                                                                                                             ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTAGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATc                         >  1:8329832/1‑150 (MQ=255)
                                                                                                                              cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCg                        >  2:2013091/1‑150 (MQ=255)
                                                                                                                              cacaGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCg                        >  2:2966205/1‑150 (MQ=255)
                                                                                                                                caGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGaa                      <  2:6893675/150‑1 (MQ=255)
                                                                                                                                caGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGaa                      >  1:9388280/1‑150 (MQ=255)
                                                                                                                                     agagggataACACACTCACGGCCTGCTCCTGGAGACGGCCTTGCTATTAGTTGCGGTACAGCACCTTATTTATGTGAAATCCGAACTGGGTGTGCAGCGGGCGGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCg                 <  1:2868904/150‑1 (MQ=255)
                                                                                                                                        gggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGaa              <  2:1477643/150‑1 (MQ=255)
                                                                                                                                         ggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAc             >  1:1939860/1‑150 (MQ=255)
                                                                                                                                         ggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAc             >  2:1159031/1‑150 (MQ=255)
                                                                                                                                         ggagaAAACACTCAAGCCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGTTGTGATATCCGACCTGGGTGTGCAGCGGGCCGGTCGGCTCCTGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAc             <  1:4302277/150‑1 (MQ=255)
                                                                                                                                          gagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAAcc            >  1:1504794/1‑150 (MQ=255)
                                                                                                                                            gaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTCTTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGCACCAt          >  2:9269786/1‑150 (MQ=255)
                                                                                                                                               aaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCa                                <  2:5179858/125‑1 (MQ=255)
                                                                                                                                                 cacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAc     >  1:9958122/1‑150 (MQ=255)
                                                                                                                                                 cacTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAc     <  2:9914765/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCT  >  NC_000913/3959353‑3959650

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