Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 304,085 A→C 6.1% I33S (ATT→AGT)  ykgJ ← UPF0153 cysteine cluster protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913304,0850AC 6.1% 362.3 / 13.0 132I33S (ATT→AGT) ykgJUPF0153 cysteine cluster protein
Reads supporting (aligned to +/- strand):  ref base A (66/58);  new base C (3/5);  total (69/63)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-01

Only 100 of 137 total aligned reads displayed.

CTCTGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTATGT  >  NC_000913/303936‑304234
                                                                                                                                                     |                                                                                                                                                     
ctctGCATGTGGACGCTCGATTTTTATAAATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGaa                                                                                                                                                       <  1:44630/150‑1 (MQ=255)
ctctGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGaa                                                                                                                                                       >  2:790349/1‑150 (MQ=255)
ctctGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGCGGCGTTTTCTGATTGGTACCGCTCATGCACCGTTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGaa                                                                                                                                                       >  1:413581/1‑150 (MQ=255)
 tctGCATGTGGACGATAGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGGTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGCTTGGTACCGCGCATGCATCGGTGAAAAGGGGATATTTGTCCAGTGAGCCTGGCGGGAAt                                                                                                                                                      <  2:81215/150‑1 (MQ=255)
      aTGTGGACGATCGATTTTTATATATCGTACAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTAc                                                                                                                                                 <  2:942113/150‑1 (MQ=255)
       tgtgGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGTGTTTTTATTATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTAcc                                                                                                                                                >  2:335969/1‑150 (MQ=255)
       tgtgGACGATCGATTTGTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTAcc                                                                                                                                                <  2:1532204/150‑1 (MQ=255)
        gtgGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCg                                                                                                                                               >  1:1370692/1‑150 (MQ=255)
         tgGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGc                                                                                                                                              >  1:773853/1‑150 (MQ=255)
          ggACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGcc                                                                                                                                             >  1:790798/1‑150 (MQ=255)
          ggACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGcc                                                                                                                                             >  2:218636/1‑150 (MQ=255)
          ggACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGcc                                                                                                                                             <  2:538345/150‑1 (MQ=255)
           gACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCa                                                                                                                                            >  1:567211/1‑150 (MQ=255)
            aagatcgttTTTTATATTTCGTGCAGCAGGCATTTTTCCCCGGGGTTCCTGCAAGGGCAATACTTCGTGTGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGTAATAGTACCGCCAg                                                                                                                                           <  2:78772/148‑1 (MQ=255)
             cgatcgatTTTTATATATCGTGCAGCTGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGTACTATTCCCGCCAgg                                                                                                                                          <  1:1002259/150‑2 (MQ=255)
             cgatcgatTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGTACTATTCCCGCCAgg                                                                                                                                          <  1:604289/150‑2 (MQ=255)
             cgatcgatTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGTACTATTCCCGCCAgg                                                                                                                                          <  2:1146678/150‑2 (MQ=255)
             cgatcgatTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTCATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTAccgcccgc                                                                                                                                          >  2:426473/1‑150 (MQ=255)
               atcgatTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCAt                                                                                                                                        >  1:401995/1‑150 (MQ=255)
                tcgatTTTTATATATCGTGCAGCAGGTATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGTACTATTCCCGCCAggctc                                                                                                                                       <  1:1226461/150‑5 (MQ=255)
                    tttttATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCAtcgtcg                                                                                                                                   >  1:94963/1‑150 (MQ=255)
                      tttATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGTGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCAtcgtctga                                                                                                                                 >  2:521988/1‑147 (MQ=255)
                        tatataTCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCtt                                                                                                                               <  1:1029272/150‑1 (MQ=255)
                         atatatCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTc                                                                                                                              <  2:767672/150‑1 (MQ=255)
                          tatatCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCg                                                                                                                             >  1:1269601/1‑150 (MQ=255)
                          tatatCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCg                                                                                                                             >  2:638428/1‑150 (MQ=255)
                            tatCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGc                                                                                                                           >  2:366521/1‑150 (MQ=255)
                             atCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGGTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGcc                                                                                                                          <  1:174972/150‑1 (MQ=255)
                               cGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCa                                                                                                                        >  1:1343309/1‑150 (MQ=255)
                                  gcagcaGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTTCGCCTGGCGTGAATAGTACCGACAGCATCGTCGGCTTCGGCCCAGTa                                                                                                                     >  2:829931/1‑150 (MQ=255)
                                  gcagcaGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTa                                                                                                                     >  1:1711308/1‑150 (MQ=255)
                                    agcagGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTaaa                                                                                                                   >  1:1482565/1‑150 (MQ=255)
                                       agGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAag                                                                                                                <  2:1338236/150‑1 (MQ=255)
                                       agGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAag                                                                                                                <  2:83402/150‑1 (MQ=255)
                                              tttGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGAcgcg                                                                                                         <  1:1028231/150‑1 (MQ=255)
                                               ttGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCgg                                                                                                        <  1:151107/150‑1 (MQ=255)
                                               ttGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCgg                                                                                                        <  2:1269601/150‑1 (MQ=255)
                                                tGCCCGGGGTGCCTGCAAGGGCAATACATCGGGGGTGTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGa                                                                                                       <  2:1606602/150‑1 (MQ=255)
                                                  cccGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGaaa                                                                                                     <  1:633801/150‑1 (MQ=255)
                                                   ccGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGaaaa                                                                                                    >  1:1574635/1‑150 (MQ=255)
                                                       ggTTCCTGCAAGGGCAATACATCGGGGGGTTTTCTGATGGGTACCGCTCATGCATCGGTGAAAAGGGGAGATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGc                                                                                                <  2:436732/150‑1 (MQ=255)
                                                        gTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGca                                                                                               >  1:848995/1‑150 (MQ=255)
                                                        gTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCAGCGTCGGCTGCGGCCCAGTAAAAAGAGCCGCGGAAAAATgac                                                                                               >  1:431591/1‑148 (MQ=255)
                                                           ccTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCacaa                                                                                            >  2:310802/1‑150 (MQ=255)
                                                           ccTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAAAAGTACCGCCAGCATCGTCGGCTTCGGCCCAGAAAAAAGAGACGCGGAAAAATGGacaa                                                                                            >  1:536631/1‑150 (MQ=255)
                                                             tGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCaccaca                                                                                          >  1:278871/1‑150 (MQ=255)
                                                              gCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAg                                                                                         >  1:1122999/1‑150 (MQ=255)
                                                                 aGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGCACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGca                                                                                      <  1:681066/150‑1 (MQ=255)
                                                                    gCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGCGCCTGGCGGGAATAGTCCCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATTCACAACAGGcaccc                                                                                   >  2:67600/1‑149 (MQ=255)
                                                                        tACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACg                                                                               <  2:1128440/150‑1 (MQ=255)
                                                                              gggggTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTcatgc                                                                         <  1:566655/150‑1 (MQ=255)
                                                                              gggggTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCAAACGTcatgc                                                                         >  2:125534/1‑150 (MQ=255)
                                                                                gggTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACCCGAcatgcat                                                                       >  1:1354244/1‑150 (MQ=255)
                                                                                   tttttCTGATTGGTACCGCTCATGCATCTGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCATCGGCCCAGTAAAAAGAGACCCGGAAAAAAGCACAAAATGCAACACACGTCATGCATTGa                                                                    >  2:107550/1‑150 (MQ=255)
                                                                                   tttttCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGa                                                                    >  2:389552/1‑150 (MQ=255)
                                                                                   tttttCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGa                                                                    >  2:1360751/1‑150 (MQ=255)
                                                                                      ttCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGATAAATGCACAACAGGCACCACACGTCATGCATGGATTc                                                                 >  2:1198918/1‑150 (MQ=255)
                                                                                        cTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAg                                                               >  2:526570/1‑150 (MQ=255)
                                                                                        cTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAg                                                               <  2:1370692/150‑1 (MQ=255)
                                                                                          gATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGAt                                                             <  2:36134/150‑1 (MQ=255)
                                                                                            ttGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTg                                                           <  2:970721/150‑1 (MQ=255)
                                                                                              ggTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCt                                                         <  2:220502/150‑1 (MQ=255)
                                                                                              ggTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCt                                                         <  2:617088/150‑1 (MQ=255)
                                                                                               gTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTc                                                        >  1:1155024/1‑150 (MQ=255)
                                                                                                   cGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATTGATTCAGATTGCTCATaa                                                    >  1:217100/1‑150 (MQ=255)
                                                                                                   cGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATaa                                                    >  1:1491499/1‑150 (MQ=255)
                                                                                                       caggCATCGGTGAAAAGGGGATATTTGGTCAGTGAGCCTGGCGGGAATATTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGACCAACCGTCACCACACGTCCGGCAGGGCTTCAGATTGCTCATTATTcc                                                >  1:1060749/4‑149 (MQ=255)
                                                                                                         tGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCAcc                                              <  2:1711308/150‑1 (MQ=255)
                                                                                                         tGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCAcc                                              <  1:790349/150‑1 (MQ=255)
                                                                                                         tGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTAGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCAcc                                              >  1:582761/1‑150 (MQ=255)
                                                                                                          gCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCa                                             >  2:756237/1‑150 (MQ=255)
                                                                                                          gCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCa                                             <  2:1314605/150‑1 (MQ=255)
                                                                                                           cATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAg                                            >  2:875839/1‑150 (MQ=255)
                                                                                                            aTCGGTGAAAAGGGTATATTTGTTAAGTGATCCTGGCGGGAATAGTACCTCAAGCATCGTCGGCTTCGGCCAAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGt                                           <  1:77235/150‑1 (MQ=255)
                                                                                                            aTCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTa                                          >  1:509641/1‑150 (MQ=255)
                                                                                                               ggTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACa                                        >  1:1311679/1‑150 (MQ=255)
                                                                                                                 ggAACAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGGAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGc                                      <  2:832056/149‑1 (MQ=255)
                                                                                                                  gAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTAGGCCCAGTAAAAAGAGACGCGGAAAAATGCACACAAGGCACCACACCTCATGAATGGATTCAGATTGCTCATTATTCACCAGGACAGca                                     >  2:880827/1‑149 (MQ=255)
                                                                                                                       ggggATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGGCAAGCATGGATTCAGATTGCTCATAATTCACCACTACAGCTaaaaa                                >  2:165640/1‑150 (MQ=255)
                                                                                                                          gATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCg                             >  1:1545360/1‑150 (MQ=255)
                                                                                                                          gATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCg                             <  1:521988/150‑1 (MQ=255)
                                                                                                                          gATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCg                             <  1:376012/150‑1 (MQ=255)
                                                                                                                          gATATTTGTTCAGTGAGACTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAATAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCTTGGATTCAGATTGCTCATACTTCACCAGTACAGCTAAAAATCg                             <  1:107550/150‑1 (MQ=255)
                                                                                                                               ttGTTCAGTGTGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCGTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCTCCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAag                        <  2:310852/150‑1 (MQ=255)
                                                                                                                               ttGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAag                        <  1:310802/150‑1 (MQ=255)
                                                                                                                               ttGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAag                        >  1:1616779/1‑150 (MQ=255)
                                                                                                                                tGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTTAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAaga                       >  1:642389/1‑150 (MQ=255)
                                                                                                                                tGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAaga                       <  2:78590/150‑1 (MQ=255)
                                                                                                                                 gTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAagag                      <  2:773853/150‑1 (MQ=255)
                                                                                                                                 gTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCACACTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATAAAGATTGCTCAAAATTCACCAGAACAGCTAAAAATCGTAAagag                      >  2:1503995/1‑150 (MQ=255)
                                                                                                                                    cAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTAACCAGTACAGCTAAAAATCGTAAAGAGAAc                   >  2:1551116/1‑150 (MQ=255)
                                                                                                                                      gTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATCGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAg                 >  1:896964/1‑150 (MQ=255)
                                                                                                                                        gatgcTGGCGGTACTATTCCCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCg               >  2:1520076/5‑150 (MQ=255)
                                                                                                                                        gatgcTGGCGGTACTATTCCCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCg               >  2:1389509/5‑150 (MQ=255)
                                                                                                                                        gatgcTGGCGGTACTATTCCCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGGCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCg               >  2:1667349/5‑150 (MQ=255)
                                                                                                                                        gAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCg               >  2:9760/1‑150 (MQ=255)
                                                                                                                                         aGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCgg              <  2:1338344/150‑1 (MQ=255)
                                                                                                                                           ccTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTa            <  2:1636472/150‑1 (MQ=255)
                                                                                                                                           ccTCGCGGTAATAGTACCGCCAGCATCGTCGGCTTCGGACGAGTAAAAAGAGACTCGGAAAAATGCACAACAGTCACCACACGTCATGAATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTa            <  1:829931/150‑1 (MQ=255)
                                                                                                                                                    aaTAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCAAACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTATg   <  1:125534/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCTGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTATGT  >  NC_000913/303936‑304234

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