Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,956,973:1 +A 8.7% coding (848/894 nt) ilvY ← DNA‑binding transcriptional dual regulator IlvY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,956,9731.A 8.7% 595.2 / 56.8 161coding (848/894 nt)ilvYDNA‑binding transcriptional dual regulator IlvY
Reads supporting (aligned to +/- strand):  ref base . (80/67);  new base A (4/10);  total (84/77)
Fisher's exact test for biased strand distribution p-value = 9.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.16e-01

Only 100 of 165 total aligned reads displayed.

TTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGCCGAGTGCC  >  NC_000913/3956824‑3957122
                                                                                                                                                      |                                                                                                                                                     
ttGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCt                                                                                                                                                        >  2:502404/1‑150 (MQ=255)
ttGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCt                                                                                                                                                        <  1:287031/150‑1 (MQ=255)
  gCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCa                                                                                                                                                     <  2:1154527/150‑1 (MQ=255)
    cACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATg                                                                                                                                                   <  1:678371/150‑1 (MQ=255)
        acgaAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGc                                                                                                                                               <  1:929636/150‑1 (MQ=255)
         cgaAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGcc                                                                                                                                              >  2:484149/1‑150 (MQ=255)
            aaCCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCt                                                                                                                                           <  2:267806/150‑1 (MQ=255)
             aCCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCt                                                                                                                                           >  1:444583/1‑149 (MQ=255)
             aCCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCt                                                                                                                                           <  2:444583/149‑1 (MQ=255)
               cAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaat                                                                                                                                        <  1:588111/150‑4 (MQ=255)
               cAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaat                                                                                                                                        <  2:401907/150‑4 (MQ=255)
               cAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaat                                                                                                                                        <  2:49068/150‑4 (MQ=255)
               cAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaat                                                                                                                                        <  2:646532/150‑4 (MQ=255)
                  tAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCt                                                                                                                                     >  1:916713/1‑150 (MQ=255)
                    aCCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTgc                                                                                                                                   >  1:1156144/1‑150 (MQ=255)
                              cAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCACTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGc                                                                                                                         >  1:1322581/1‑150 (MQ=255)
                                    cTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCt                                                                                                                   >  2:357611/1‑150 (MQ=255)
                                          ggCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACg                                                                                                             <  1:1083013/150‑1 (MQ=255)
                                             gggTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCg                                                                                                          <  2:1206233/150‑1 (MQ=255)
                                                  aGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtc                                                                                                     <  2:614883/150‑1 (MQ=255)
                                                  aGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtc                                                                                                     <  1:798759/150‑1 (MQ=255)
                                                   gggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtct                                                                                                    <  1:1216360/150‑1 (MQ=255)
                                                   gggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtct                                                                                                    >  2:1060568/1‑150 (MQ=255)
                                                   gggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtct                                                                                                    >  2:818495/1‑150 (MQ=255)
                                                   gggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtct                                                                                                    >  2:144675/1‑150 (MQ=255)
                                                   gggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGCCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtct                                                                                                    >  1:13010/1‑150 (MQ=255)
                                                   gggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGCCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtct                                                                                                    >  1:818813/1‑150 (MQ=255)
                                                    ggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtctc                                                                                                   >  1:892433/1‑150 (MQ=255)
                                                       aaaaTGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATc                                                                                                <  2:537828/150‑1 (MQ=255)
                                                          aTGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACCCGCCCAGCTAAAACGGCGGTTGCTCAACgag                                                                                             >  1:771001/1‑148 (MQ=255)
                                                           tGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTg                                                                                            <  2:292702/150‑1 (MQ=255)
                                                            gCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTgc                                                                                           >  1:610290/1‑150 (MQ=255)
                                                            gCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTgc                                                                                           >  2:675059/1‑150 (MQ=255)
                                                               tGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGct                                                                                        >  2:472816/1‑150 (MQ=255)
                                                                gATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGGCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGctc                                                                                       <  1:973663/150‑1 (MQ=255)
                                                                  tAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTa                                                                                     >  1:91592/1‑150 (MQ=255)
                                                                   aGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGGTTTCTCATCGCTGCGCTCTaa                                                                                    <  1:158335/150‑1 (MQ=255)
                                                                    gcgcTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTg                                                                                            <  2:1260677/141‑1 (MQ=255)
                                                                    gcgcTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTg                                                                                            >  1:1260677/1‑141 (MQ=255)
                                                                     cgcttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTaaaa                                                                                  <  2:282737/150‑1 (MQ=255)
                                                                       cttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAatc                                                                                >  2:683351/1‑150 (MQ=255)
                                                                       cttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAatc                                                                                >  2:688304/1‑150 (MQ=255)
                                                                       cttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAatc                                                                                <  1:725057/150‑1 (MQ=255)
                                                                          cgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAatcatc                                                                             <  2:610290/150‑1 (MQ=255)
                                                                          cgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAatcatc                                                                             >  1:251566/1‑150 (MQ=255)
                                                                          cgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAatcatc                                                                             >  1:501378/1‑150 (MQ=255)
                                                                           gcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAACCGGCGTTTTCTCATCGCTGCGCTCTAAAAtcatca                                                                            >  2:346189/1‑150 (MQ=255)
                                                                            cttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCAc                                                                           >  1:192268/1‑150 (MQ=255)
                                                                            cttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCAc                                                                           >  1:473230/1‑150 (MQ=255)
                                                                            cttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCAc                                                                           >  2:432514/1‑150 (MQ=255)
                                                                              tATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCAcgc                                                                         <  1:472816/150‑1 (MQ=255)
                                                                               aTCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCTGCAAAATCTTCCATAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTgcg                                                                                          <  1:955624/132‑1 (MQ=255)
                                                                               aTCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCAcgcg                                                                        <  1:695989/150‑1 (MQ=255)
                                                                               aTCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCAGCAAAATCTTCCATAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTgcg                                                                                          >  2:955624/1‑132 (MQ=255)
                                                                                   ggCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTa                                                                    >  2:1280958/1‑150 (MQ=255)
                                                                                   ggCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTa                                                                    <  2:1322581/150‑1 (MQ=255)
                                                                                    gCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTAc                                                                   >  1:95545/1‑150 (MQ=255)
                                                                                    gCCTACCCGCGCGACAACGCCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTAc                                                                   <  2:818813/150‑1 (MQ=255)
                                                                                     ccTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACg                                                                  <  1:1060568/150‑1 (MQ=255)
                                                                                         cccGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCAcc                                                              >  2:770361/1‑150 (MQ=255)
                                                                                          ccGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGTGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCg                                                             >  1:602139/1‑150 (MQ=255)
                                                                                           cgcgcgACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCgg                                                            >  2:13488/1‑150 (MQ=255)
                                                                                           cgcgcgACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCgg                                                            <  1:688304/150‑1 (MQ=255)
                                                                                           cgcgcgACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCgg                                                            <  1:484149/150‑1 (MQ=255)
                                                                                                  cAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCgggg                                                     <  2:566919/150‑1 (MQ=255)
                                                                                                   aaCGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGc                                                    >  2:860476/1‑150 (MQ=255)
                                                                                                    aCGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCt                                                   >  2:858542/1‑150 (MQ=255)
                                                                                                      gTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGt                                                 <  1:358242/150‑1 (MQ=255)
                                                                                                        cATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGttt                                               >  1:1397742/1‑150 (MQ=255)
                                                                                                          tttGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTc                                             >  2:1112476/1‑150 (MQ=255)
                                                                                                          tttGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTc                                             >  1:894077/1‑150 (MQ=255)
                                                                                                           ttGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTGCGCACCGGTTCGGGGCTGTTTTcc                                            <  1:1031093/150‑1 (MQ=255)
                                                                                                           ttGTGGTGCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTcc                                            <  1:1023934/150‑1 (MQ=255)
                                                                                                              tgGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGc                                         <  2:1223939/150‑1 (MQ=255)
                                                                                                              tgGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGc                                         >  2:1216726/1‑150 (MQ=255)
                                                                                                                 ttCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGGTGTTTTCCAGcacc                                      <  2:882429/150‑1 (MQ=255)
                                                                                                                 ttCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGcacc                                      <  2:1156144/150‑1 (MQ=255)
                                                                                                                   cGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGcaccac                                    <  2:192268/150‑1 (MQ=255)
                                                                                                                    ggCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACt                                   >  2:781340/1‑150 (MQ=255)
                                                                                                                     gCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACtt                                  <  1:860476/150‑1 (MQ=255)
                                                                                                                       aaaaTCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTcc                                >  2:833787/1‑150 (MQ=255)
                                                                                                                       aaaaTCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTcc                                >  2:518936/1‑150 (MQ=255)
                                                                                                                          aTCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGc                             <  2:916713/150‑1 (MQ=255)
                                                                                                                           tCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGca                            <  1:826175/150‑1 (MQ=255)
                                                                                                                               ccAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGcaacaa                        >  2:423327/1‑150 (MQ=255)
                                                                                                                                  gAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgc                     <  2:1397742/150‑1 (MQ=255)
                                                                                                                                   aaTGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgcc                    <  2:1270670/150‑1 (MQ=255)
                                                                                                                                       ccTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCTCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgccacgc                >  1:842878/1‑150 (MQ=255)
                                                                                                                                       ccTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgccacgc                >  2:1052761/1‑150 (MQ=255)
                                                                                                                                        cagaaaaaGCGGCTACATGAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgccacgcc               >  2:570569/8‑150 (MQ=255)
                                                                                                                                        cTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCCCCACGTCCGGCAACAATgccacgcc               >  1:1315312/1‑150 (MQ=255)
                                                                                                                                        cTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgccacgcc               >  1:638587/1‑150 (MQ=255)
                                                                                                                                            aaaaGCGGCTACATGAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAg           >  1:484360/4‑150 (MQ=255)
                                                                                                                                             ttAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGc          <  1:432514/150‑1 (MQ=255)
                                                                                                                                              tAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGcc         >  1:1071423/1‑150 (MQ=255)
                                                                                                                                                gCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGCCGa       >  1:684370/1‑150 (MQ=255)
                                                                                                                                                 cGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGCCGAg      <  1:409127/150‑1 (MQ=255)
                                                                                                                                                 cGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGCCGAg      <  2:1220299/150‑1 (MQ=255)
                                                                                                                                                  ggCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAACTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCCCCACTTCCGGCAACAATGCCACGCCACCGCCGAGt     >  2:740436/1‑150 (MQ=255)
                                                                                                                                                     tCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGCCGAGTGcc  <  2:91830/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGCCGAGTGCC  >  NC_000913/3956824‑3957122

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