Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,728,120 T→G 5.3% T418P (ACG→CCG)  xylB ← xylulokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,728,1200TG 5.3% 1166.6 / 7.1 378T418P (ACG→CCG) xylBxylulokinase
Reads supporting (aligned to +/- strand):  ref base T (119/228);  new base G (28/2);  total (148/230)
Fisher's exact test for biased strand distribution p-value = 1.23e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 388 total aligned reads displayed.

TGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGC  >  NC_000913/3727971‑3728269
                                                                                                                                                      |                                                                                                                                                     
  gCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTAc                                                                                                                                                     <  1:1629618/150‑1 (MQ=255)
    tGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGTGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACgg                                                                                                                                                   <  2:1503275/150‑1 (MQ=255)
    tGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACTGCGCCCGCCGTGCgg                                                                                                                                                   >  1:1238206/1‑150 (MQ=255)
    tGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACgg                                                                                                                                                   >  2:890216/1‑150 (MQ=255)
    tGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACgg                                                                                                                                                   >  2:3557570/1‑150 (MQ=255)
      aTAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTa                                                                                                                                                 <  2:5427794/150‑1 (MQ=255)
       tAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTaa                                                                                                                                                <  2:106564/150‑1 (MQ=255)
        aaGCGGCATACCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAAt                                                                                                                                               <  2:6339598/150‑1 (MQ=255)
         aGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATc                                                                                                                                              <  1:4342571/150‑1 (MQ=255)
          gCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCGCCCGCCGGAGGGtcggcg                                                                                                                                             >  2:6293044/1‑145 (MQ=255)
           cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGa                                                                                                                                            <  1:3652580/150‑1 (MQ=255)
           cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGa                                                                                                                                            <  1:5808157/150‑1 (MQ=255)
           cGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGa                                                                                                                                            <  2:3165152/150‑1 (MQ=255)
             gCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgc                                                                                                                                          <  2:6044203/150‑1 (MQ=255)
             gCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgc                                                                                                                                          <  1:1512119/150‑1 (MQ=255)
              cATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgct                                                                                                                                         <  1:2255327/150‑1 (MQ=255)
                 aaCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACCGCCCCCCCCGGGCGGTCATCGGgctgcg                                                                                                                                      >  1:6139745/1‑149 (MQ=255)
                 aaCGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAgctgct                                                                                                                                      <  2:1945890/150‑1 (MQ=255)
                   cGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGa                                                                                                                                    <  1:3435468/150‑1 (MQ=255)
                   cGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCAAATCCCCCCCCGTACGGTATTCGAGCTGCTGa                                                                                                                                    <  2:6229580/150‑1 (MQ=255)
                    gCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGAc                                                                                                                                   <  2:4440120/150‑1 (MQ=255)
                     cTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGGCGTCGAgctgcggggc                                                                                                                                  >  2:2755965/1‑145 (MQ=255)
                        cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCt                                                                                                                               <  1:2665632/150‑1 (MQ=255)
                        cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGCCGCGAGCTGCGGGCCGcg                                                                                                                               >  1:4692530/1‑149 (MQ=255)
                         gcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTg                                                                                                                              <  2:5924576/150‑1 (MQ=255)
                         gcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTg                                                                                                                              <  2:6538883/150‑1 (MQ=255)
                          cATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGa                                                                                                                             <  1:681866/150‑1 (MQ=255)
                          cATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAAATCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGa                                                                                                                             <  2:5555334/150‑1 (MQ=255)
                                  aGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGcc                                                                                                                     <  1:5579227/150‑1 (MQ=255)
                                       cGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCAt                                                                                                                <  2:5225088/150‑1 (MQ=255)
                                        gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATc                                                                                                               >  2:1839485/1‑150 (MQ=255)
                                        gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGTGCTCGGGCTGCTGACCGCTGATGTCCGCCGGccgc                                                                                                               >  2:4227386/1‑147 (MQ=255)
                                           tGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGa                                                                                                            <  2:6121496/150‑1 (MQ=255)
                                             ttCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACg                                                                                                          <  1:2613361/150‑1 (MQ=255)
                                              tCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                                                         <  2:4585601/150‑1 (MQ=255)
                                              tCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                                                         <  1:3237335/150‑1 (MQ=255)
                                               cTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                                                         <  1:5943840/150‑1 (MQ=255)
                                                tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                                       <  1:5439029/150‑1 (MQ=255)
                                                tAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                                       <  1:6048365/150‑1 (MQ=255)
                                                 aaCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAg                                                                                                      <  2:2649786/150‑1 (MQ=255)
                                                 aaCGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGGCATCGCGGTGCTGGCCGCTGATCTCCGCCCGCCTCGGGCGccgg                                                                                                      >  2:204198/1‑148 (MQ=255)
                                                       aGTTGCGGCAACAATTCAATGAGCGCTTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGATGATATCCGCCAGCATCTGACGCCAGTACTCa                                                                                                <  2:6589593/150‑1 (MQ=255)
                                                       aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCa                                                                                                <  1:777416/150‑1 (MQ=255)
                                                           gCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTAc                                                                                            <  1:3787472/150‑1 (MQ=255)
                                                             ggCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTAcgc                                                                                          <  1:4251124/150‑1 (MQ=255)
                                                                   aaTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGGGCGGTGATCGGGCTGCGGTCCGGTGATATCCGCCTGCGTCTGTCGCGGGTAGTGACTGCGcgccccg                                                                                    >  1:2107563/1‑150 (MQ=255)
                                                                    aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGTATCTGACGCCAGTACTCACTACGcgccccgc                                                                                   <  1:1067484/150‑1 (MQ=255)
                                                                    aTTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgc                                                                                   <  2:5157512/150‑1 (MQ=255)
                                                                     ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACg                                                                                                                                                    >  1:5229101/1‑84 (MQ=255)
                                                                     ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgcc                                                                                  <  1:2929670/150‑1 (MQ=255)
                                                                     ttCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTCGCTGCGCCCAGTGCTGGCCCCACCTCCCCCCCCGTTCGGGAATGGGGCTGCTGTGCGGTGATTTCCGCCAGCAGCTGACGCCGGGACACACTACGcgccccgcc                                                                                  >  1:2307386/1‑150 (MQ=255)
                                                                      tCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgcc                                                                                  <  2:692246/150‑1 (MQ=255)
                                                                      tCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGcgccccgccc                                                                                 <  1:4607414/150‑1 (MQ=255)
                                                                       cAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCgccccgcccc                                                                                <  2:1431944/150‑1 (MQ=255)
                                                                        aaTGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCCCCGTGCGGTCGTCGAGCTGGTGGCCGCTGATATCCGCCGGCGGCTGGCGGCGGTGCGGACTGCGCGCGCCGccccc                                                                               >  2:1317390/1‑150 (MQ=255)
                                                                          tGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCaa                                                                             <  2:3760658/150‑1 (MQ=255)
                                                                           gAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAGCCCCCGCCGTGCGGTAGGCGCGCTGGTGGCCGCTGGTATCCGGCAGCGTCTGGCGGGGGTGCTCGCTGGGCGCCCCGCCCCCacg                                                                            >  1:3517602/1‑148 (MQ=255)
                                                                            aGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATc                                                                           <  1:3898812/150‑1 (MQ=255)
                                                                             gCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTGCGGGCCGCGCGCTGCGGACCGGTGCGCTCCGCCCGCAGCGGGCGCCGGTGCTCGCTACGCGCCCCGccccccctcg                                                                          >  2:4885670/1‑145 (MQ=255)
                                                                               gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCTAc                                                                        <  2:6635911/150‑1 (MQ=255)
                                                                                     tctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAaca                                                                  <  2:1683495/150‑1 (MQ=255)
                                                                                       tGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCt                                                                                                                               <  2:2614897/87‑1 (MQ=255)
                                                                                        ggATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACtt                                                               <  2:2696448/150‑1 (MQ=255)
                                                                                            tCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCg                                                           <  1:4115985/150‑1 (MQ=255)
                                                                                                cgcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGGCGGTGCTCGAGCTGGTGACCGCTGTTATCCGCCAGCCTCTGTCGCCGGTACTCGCGACGCGCCCCGCCCCCGATCAGCGGAGCGCTTTGCGGtgg                                                       >  2:3319533/1‑148 (MQ=255)
                                                                                                 gcgATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGttt                                                       <  1:2304024/150‑1 (MQ=255)
                                                                                                   gATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAAt                                                    <  1:6606364/150‑1 (MQ=255)
                                                                                                    cgcTGCGCCAGCCTTGCTGCGCCCAGTGCGGGCCCCACCTCCCCCCCCGTCCGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACCCTTTGCGGTTTAATa                                                   <  2:3684101/148‑1 (MQ=255)
                                                                                                    aTCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATa                                                   <  2:3459626/150‑1 (MQ=255)
                                                                                                       tGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGGACGGTAATCGAGCTGCGGGCCGCTGATCTCCGCCAGCGTCTGGCGCCCGTACGCACTGCGCGCCCCGCCCCCCCGCAACGTAACCCTTTGCGGTGTAATACCg                                                >  2:1744367/1‑150 (MQ=255)
                                                                                                        gcgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGc                                               <  2:2553218/150‑1 (MQ=255)
                                                                                                         cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCa                                              <  2:1778179/150‑1 (MQ=255)
                                                                                                         cgcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTATTAAACCGCa                                              <  1:2473745/150‑1 (MQ=255)
                                                                                                          gcCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAg                                             <  1:859062/150‑1 (MQ=255)
                                                                                                            cAGCCTTGCTGCGCCCAGTGCTGGCCCCACCTCCCCGCCCGTTGGGGTGTCGGGCTGCTGGCCGCGGGTATCCGCCGGCGTCGGGGGCCGGTACTCACTACGCGCCCCGCGCCCAATCGACGGAACCCGTTGCGGGTTAATACCGCGGGc                                           >  2:3972539/1‑150 (MQ=255)
                                                                                                              gCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCAt                                         <  1:6293044/150‑1 (MQ=255)
                                                                                                                cTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGc                                       <  1:5961339/150‑1 (MQ=255)
                                                                                                                 ttGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCa                                      <  1:769066/150‑1 (MQ=255)
                                                                                                                     tgcGCCCAGTGCTGGCCCCACATCCCCCCCCGTCCGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCCTCTGGCGCCAGTACTCACTCCGCGCCCCGCCCCCCATCAACGTCACACGTTGCGGTTTAATGCCGCAGGCGGGCACGAc                                  >  2:5796836/1‑150 (MQ=255)
                                                                                                                      gcgcCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTCCTCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGTCGCCAGTAGTCACTCCGCGCCCCGCCCCCCATCAACGTCCCCCGTTGCGGTTGAATACCGCTGGCAGGCACGACa                                 >  2:2888770/1‑150 (MQ=255)
                                                                                                                         cccAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCCGCATCTGACGCCAGTACTCCCTACGCGCCCCGCCCCCAATCCACGTAACACTTTGCGGTGTGATACCGCAGGCATGCACGTCATcc                              >  2:759761/1‑150 (MQ=255)
                                                                                                                          ccAGTGCTGGCCCCACCTCCCCCCCCGTCCGGTACTCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCCGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACGTTGCGGTTTAATACCGCAGGCATGCACGACATCCa                             >  1:4995459/1‑150 (MQ=255)
                                                                                                                            aGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATg                           <  2:4153861/150‑1 (MQ=255)
                                                                                                                             gTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGc                          <  1:3853425/150‑1 (MQ=255)
                                                                                                                                cTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCAt                       <  2:5479983/150‑1 (MQ=255)
                                                                                                                                 tGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCGGATATCCGCCAGCATCTGACGCCAGTACGCACTACGCGCCCCGCCCCCAATCAACGGACCACTTTGCGGTGTAATACCGCAGGCATGCACGACATCCCTGCCATc                      >  2:1180568/1‑150 (MQ=255)
                                                                                                                                     cccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGAAGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGAAATCCAAGCCATCTGcc                  >  2:4694693/1‑150 (MQ=255)
                                                                                                                                      ccACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCa                 >  1:1023992/1‑150 (MQ=255)
                                                                                                                                             ccccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAACAACGTAACACTTTGCGGTTTAAAACCGCAGGCATGCACGACACCCATGCCATCTGCCAGCGCaa           >  1:1795621/1‑149 (MQ=255)
                                                                                                                                               ccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATAc                                                  >  2:5857818/1‑109 (MQ=255)
                                                                                                                                               ccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATaa         <  2:503925/150‑1 (MQ=255)
                                                                                                                                               ccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATaa         >  2:3844980/1‑150 (MQ=255)
                                                                                                                                               ccccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATaa         <  1:5040618/150‑1 (MQ=255)
                                                                                                                                                 ccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCa                                      >  2:3092930/1‑119 (MQ=255)
                                                                                                                                                 ccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCa                                      <  1:3092930/119‑1 (MQ=255)
                                                                                                                                                 ccccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAAcc       >  2:4037480/1‑150 (MQ=255)
                                                                                                                                                   ccGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAgc               >  2:2390278/1‑140 (MQ=255)
                                                                                                                                                    cGTACGGTAAGCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCAc    <  2:2298899/150‑1 (MQ=255)
                                                                                                                                                      tACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGc  <  1:6576133/150‑1 (MQ=255)
                                                                                                                                                      tACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGc  <  2:5720989/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATCCATGCCATCTGCCAGCGCATAACCCACGC  >  NC_000913/3727971‑3728269

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: