Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,907,998 A→T 15.6% intergenic (‑57/+31) eno ← / ← pyrG enolase/CTP synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,907,9980AT15.6% 244.1 / 34.7 86intergenic (‑57/+31)eno/pyrGenolase/CTP synthetase
Reads supporting (aligned to +/- strand):  ref base A (42/30);  new base T (7/6);  total (50/36)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

ACCACCCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCAC  >  NC_000913/2907852‑2908146
                                                                                                                                                  |                                                                                                                                                    
accaccCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAgg                                                                                                                                                    <  1:1149147/149‑1 (MQ=255)
accaccCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAgg                                                                                                                                                    <  3:624188/149‑1 (MQ=255)
 ccaccCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAggg                                                                                                                                                   <  1:866564/149‑1 (MQ=255)
  caccCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGt                                                                                                                                                  >  3:712837/1‑149 (MQ=255)
  caccCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCAGACAAACAACGCGTACCCAGGGTa                                                                                                                                                 <  2:22274/149‑1 (MQ=255)
      cTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTAcgc                                                                                                                                              >  4:412594/1‑149 (MQ=255)
      cTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTAcgc                                                                                                                                              <  3:949739/149‑1 (MQ=255)
       tCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTAcgcg                                                                                                                                             >  2:651397/1‑149 (MQ=255)
        ccAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTACGCGt                                                                                                                                            <  3:412594/149‑1 (MQ=255)
        ccAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTACGCGt                                                                                                                                            <  1:782594/149‑1 (MQ=255)
        ccAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGt                                                                                                                                            <  3:524942/149‑1 (MQ=255)
                 cTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGctct                                                                                                                                   <  1:65732/149‑1 (MQ=255)
                  ttCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGtt                                                                                                                                           >  3:210653/1‑140 (MQ=255)
                  ttCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGtt                                                                                                                                           <  4:210653/140‑1 (MQ=255)
                  ttCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGt                                                                                                                                            >  3:393510/1‑139 (MQ=255)
                  ttCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGt                                                                                                                                            <  4:393510/139‑1 (MQ=255)
                    cGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAAc                                                                                                                                <  1:651397/149‑1 (MQ=255)
                         tCAACAGTGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACtttttt                                                                                                                          <  2:752510/149‑1 (MQ=255)
                              aGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttactt                                                                                                                      >  1:346339/1‑149 (MQ=255)
                               gTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttactta                                                                                                                     >  2:556796/1‑149 (MQ=255)
                               gTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttactta                                                                                                                     >  2:436362/1‑149 (MQ=255)
                                tCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttacttac                                                                                                                    <  2:269678/149‑1 (MQ=255)
                                  gggTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTACGCGTTgtt                                                                                                                                        <  1:804059/127‑3 (MQ=255)
                                  gggTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTACGCGTTgtt                                                                                                                                        >  2:804059/1‑125 (MQ=255)
                                   ggTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAgg                                                                                                                                                    >  1:573475/1‑114 (MQ=255)
                                   ggTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAgg                                                                                                                                                    <  2:573475/114‑1 (MQ=255)
                                       aCCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCt                                                                                                             >  2:2468/1‑149 (MQ=255)
                                          aCGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGAc                                                                                                          >  1:183429/1‑149 (MQ=255)
                                          aCGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGAc                                                                                                          >  1:324276/1‑149 (MQ=255)
                                          aCGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGAc                                                                                                          >  1:1129264/1‑149 (MQ=255)
                                              gAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGttt                                                                                                      <  1:436362/149‑1 (MQ=255)
                                               aGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTc                                                                                                     <  1:2468/149‑1 (MQ=255)
                                                   gatgatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTACGCGTTgttt                                                                                                                                       >  3:382389/1‑108 (MQ=255)
                                                   gatgatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTACGCGTTgttt                                                                                                                                       <  4:382389/111‑4 (MQ=255)
                                                   gatgatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTACGCGTTgt                                                                                                                                         >  2:717723/1‑108 (MQ=255)
                                                   gatgatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCTGGGTACGCGTTgt                                                                                                                                         <  1:717723/109‑2 (MQ=255)
                                                    atgatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttactt                                                                                                                      <  4:457495/127‑1 (MQ=255)
                                                    atgatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttactt                                                                                                                      >  3:457495/1‑127 (MQ=255)
                                                    atgatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGaa                                                                                                >  3:68655/1‑149 (MQ=255)
                                                    atgatTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGaa                                                                                                >  4:622628/1‑149 (MQ=255)
                                                         ttCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGc                                                                                           >  2:400944/1‑149 (MQ=255)
                                                                   atgatTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCttt                                                                                 <  1:400944/149‑1 (MQ=255)
                                                                           tACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGc                                                                         <  2:1129264/149‑1 (MQ=255)
                                                                           tACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGc                                                                         >  4:904103/1‑149 (MQ=255)
                                                                               aTTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGc                                                                     <  1:556796/149‑1 (MQ=255)
                                                                                          tAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGt                                                          >  2:1035662/1‑149 (MQ=255)
                                                                                             gTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGa                                                                                                           <  1:31108/97‑1 (MQ=255)
                                                                                             gTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGa                                                                                                           >  2:31108/1‑97 (MQ=255)
                                                                                               tttCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCCGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCa                                                     >  2:982716/1‑149 (MQ=255)
                                                                                                  ccTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCa                                                  >  4:773530/1‑149 (MQ=255)
                                                                                                   cTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCAc                                                 >  1:968185/1‑149 (MQ=255)
                                                                                                    tCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACaaa                                                                           >  4:471204/1‑122 (MQ=255)
                                                                                                    tCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACaaa                                                                           <  3:471204/122‑1 (MQ=255)
                                                                                                    tCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACg                                                >  3:126384/1‑149 (MQ=255)
                                                                                                    tCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACg                                                >  4:893813/1‑149 (MQ=255)
                                                                                                           aCTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTa                                         >  1:193041/1‑149 (MQ=255)
                                                                                                            cTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCttt                                                                                 >  1:153907/1‑108 (MQ=255)
                                                                                                            cTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCttt                                                                                 <  2:153907/108‑1 (MQ=255)
                                                                                                            cTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAg                                        >  2:482368/1‑149 (MQ=255)
                                                                                                             tAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGa                                       >  2:556766/1‑149 (MQ=255)
                                                                                                               gTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACaaa                                                                           >  3:997741/1‑111 (MQ=255)
                                                                                                               gTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACaaa                                                                           <  4:997741/111‑1 (MQ=255)
                                                                                                                ttAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttac                                                                                                                        >  4:1116720/1‑65 (MQ=255)
                                                                                                                ttAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTttac                                                                                                                        <  3:1116720/65‑1 (MQ=255)
                                                                                                                ttAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGt                                    >  3:303833/1‑149 (MQ=255)
                                                                                                                ttAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGt                                    >  4:1028288/1‑149 (MQ=255)
                                                                                                                    aCTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAAc                                >  2:622260/1‑149 (MQ=255)
                                                                                                                      tGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTc                              >  3:2107/1‑149 (MQ=255)
                                                                                                                         aaCTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAGGTAAACTCCgg                           <  3:904103/149‑1 (MQ=255)
                                                                                                                            tCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATg                        >  3:950819/1‑149 (MQ=255)
                                                                                                                                gACAAACAACGTGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAAc                    >  3:1004442/1‑149 (MQ=255)
                                                                                                                                 acaaacaaCGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACt                   >  1:858448/1‑149 (MQ=255)
                                                                                                                                   aaacaaCGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTgg                 <  4:712837/149‑1 (MQ=255)
                                                                                                                                     gcaacGCGTACCCTGGGTACGCGTTGTCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCa               >  3:13089/2‑149 (MQ=255)
                                                                                                                                     acaacGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTg                  >  2:351770/1‑146 (MQ=255)
                                                                                                                                     acaacGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTg                  <  1:351770/146‑1 (MQ=255)
                                                                                                                                      caacGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCaa              <  3:1028288/149‑1 (MQ=255)
                                                                                                                                       aacGCGTACCCTGGGTACACGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCAGGTGACCATTACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAg             >  4:87702/1‑149 (MQ=255)
                                                                                                                                        acGCGTACCCTGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGc            >  1:355790/1‑149 (MQ=255)
                                                                                                                                        acGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGc            >  4:1037771/1‑149 (MQ=255)
                                                                                                                                        acGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGc            <  3:773530/149‑1 (MQ=255)
                                                                                                                                         cgcgTACCCTGGGTACGCGTTGTCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGcc           <  4:13089/149‑1 (MQ=255)
                                                                                                                                          gcgTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCa          >  3:408233/1‑149 (MQ=255)
                                                                                                                                              aCCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGcc           <  1:153248/144‑1 (MQ=255)
                                                                                                                                              aCCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGcc           >  2:153248/1‑144 (MQ=255)
                                                                                                                                               cccAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGTCGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAGACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAAc     >  4:954335/1‑149 (MQ=255)
                                                                                                                                                  aGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCAc  >  2:328537/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ACCACCCTCCAGATGTACTTCGGCTTCAACAGTCGGGTTACCACGGGAGTCGATGATTTCACGACCGATGATTTTTACGATTTTGGACATTAGGTTTTCCTCAAGTCACTAGTTAAACTGAAACTCCAGACAAACAACGCGTACCCAGGGTACGCGTTGCCGCTCTAACTTTTTTACTTACTTCGCCTGACGTTTCTGGAACTCGCTGGCGGCTTTCACAAAGCCTGCAAACAGCGGGTGACCATCACGTGGAGTAGAAGTAAACTCCGGATGGAACTGGCAAGCCACGAACCAC  >  NC_000913/2907852‑2908146

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: