Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,620,224 C→G 10.2% intergenic (+17/+14) marB → / ← eamA periplasmic mar operon regulator/cysteine and O‑acetyl‑L‑serine efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,620,2240CG10.2% 292.2 / 21.4 98intergenic (+17/+14)marB/eamAperiplasmic mar operon regulator/cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base C (33/54);  new base G (3/7);  total (36/62)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

TACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCC  >  NC_000913/1620075‑1620370
                                                                                                                                                     |                                                                                                                                                   
tACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAc                                                                                                                                                     <  1:1496208/150‑1 (MQ=255)
 aCTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGACCCCCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGTTAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGCCa                                                                                                                                                    >  1:1157492/1‑150 (MQ=255)
 aCTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGGTGGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACa                                                                                                                                                    <  2:1689182/150‑1 (MQ=255)
   ttcttGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtc                                                                                                                                                  <  2:439729/150‑1 (MQ=255)
   ttcttGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGAAGTTTGTTCTGGCCCCGACAtc                                                                                                                                                  >  2:993182/1‑150 (MQ=255)
     cttGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTACATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtctc                                                                                                                                                <  1:625060/150‑1 (MQ=255)
     cttGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtctc                                                                                                                                                <  1:383547/150‑1 (MQ=255)
        gtgCCAATGTCGTGGTTGTTCCGCCATCGCTGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCggg                                                                                                                                             <  1:119857/150‑1 (MQ=255)
        gtgCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACACTATGTAGTTTGTTCTGGCCCCGACATCTCggg                                                                                                                                             <  2:734142/150‑1 (MQ=255)
         tgCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCgggg                                                                                                                                            <  1:1433241/150‑1 (MQ=255)
          gCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGc                                                                                                                                           <  1:1145240/150‑1 (MQ=255)
           ccAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGcc                                                                                                                                          <  2:1405262/150‑2 (MQ=255)
            cAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCtt                                                                                                                                         >  1:1014017/1‑150 (MQ=255)
             aaTGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccag                                                                                                                                        <  2:1357102/150‑4 (MQ=255)
             aaTGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccag                                                                                                                                        <  2:1363761/150‑4 (MQ=255)
             aaTGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccag                                                                                                                                        <  1:1431830/150‑4 (MQ=255)
             aaTGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGccag                                                                                                                                        <  1:705283/150‑4 (MQ=255)
             aaTGGCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTAAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccag                                                                                                                                        <  1:1287964/150‑4 (MQ=255)
              aTGTCGTGGTTGTTGCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGTGATGTCGGGGccaga                                                                                                                                       <  2:281812/150‑5 (MQ=255)
               tGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccagaa                                                                                                                                      <  1:1350842/150‑6 (MQ=255)
                  cGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGGCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAAc                                                                                                                                   <  2:321929/150‑1 (MQ=255)
                   gTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACt                                                                                                                                  <  1:525665/150‑1 (MQ=255)
                    tGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACtt                                                                                                                                 <  2:438221/150‑1 (MQ=255)
                    tGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCGCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACtt                                                                                                                                 <  2:372199/150‑1 (MQ=255)
                     ggTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTc                                                                                                                                >  1:248174/1‑150 (MQ=255)
                      gttgttCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTcc                                                                                                                               <  1:891507/150‑1 (MQ=255)
                       ttgttCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTccc                                                                                                                              <  1:202023/150‑1 (MQ=255)
                         gttCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCACCACGCTATGTAGTGTGTTCTGGCCCGGACATCGCGGGGCTTATTAACTTCCCAc                                                                                                                            >  2:1129984/1‑150 (MQ=255)
                            cccccATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCtt                                                                                                                         >  1:740520/1‑150 (MQ=255)
                             ccccATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCttt                                                                                                                        >  1:625718/1‑150 (MQ=255)
                              cccATCGCAGGAACACCCACCGTTTGATTTAAATCACATGAGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTa                                                                                                                       <  2:220053/150‑1 (MQ=255)
                               ccATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                                                      <  1:1588286/150‑1 (MQ=255)
                               ccATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAc                                                                                                                      <  2:1558885/150‑1 (MQ=255)
                                cATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAcc                                                                                                                     >  2:646946/1‑150 (MQ=255)
                                 aTCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCTGCGTGCCCTATTATAATCAACACGATATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCg                                                                                                                    >  1:77990/1‑150 (MQ=255)
                                  tCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGc                                                                                                                   <  1:539936/150‑1 (MQ=255)
                                   cGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTAGGTAGTGTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCt                                                                                                                  >  2:825386/1‑150 (MQ=255)
                                    gCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCtt                                                                                                                 >  1:740508/1‑150 (MQ=255)
                                       ggAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTAc                                                                                                              <  1:1410110/150‑1 (MQ=255)
                                           cacccaccGTTTGATTTTAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGccc                                                                                                          >  1:315648/1‑149 (MQ=255)
                                           cacccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCa                                                                                                          >  1:1398444/1‑150 (MQ=255)
                                              ccaccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCg                                                                                                       <  1:1278648/150‑1 (MQ=255)
                                               caccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGc                                                                                                      <  1:893440/150‑1 (MQ=255)
                                               caccGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGc                                                                                                      >  2:818295/1‑150 (MQ=255)
                                                   gTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGc                                                                                                  <  1:818295/150‑1 (MQ=255)
                                                      tGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCGTAATATAATCCACACGCTATGGAGTTTGTTCTGGCCCCGACAGCTCGGGTAGTATTACCTTGACACCTTTAACGCTTAACGCCACCGCAAGCCaa                                                                                               >  1:1165290/1‑150 (MQ=255)
                                                       gATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTACCTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCaaa                                                                                              >  2:193959/1‑150 (MQ=255)
                                                        aTTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAAt                                                                                             <  1:1457114/150‑1 (MQ=255)
                                                              aTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTg                                                                                       <  2:730613/150‑1 (MQ=255)
                                                              aTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTg                                                                                       <  2:35679/150‑1 (MQ=255)
                                                               tCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGa                                                                                      <  2:908761/150‑1 (MQ=255)
                                                                  caTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGatat                                                                                   <  2:1365714/150‑1 (MQ=255)
                                                                      ggTACTGGCAGTGTTAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑Ag                                                                              <  2:1223421/150‑1 (MQ=255)
                                                                       gTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTTTTATAATCAACACGCTATGTAGATTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCGACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGTTATAC‑AGc                                                                             <  2:821908/150‑1 (MQ=255)
                                                                         acaggCAGTGATAAGTCGGATGCGCTCGGTGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGccc                                                                           <  2:702525/147‑1 (MQ=255)
                                                                               aGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCAt                                                                     <  2:1635184/150‑1 (MQ=255)
                                                                                gTGATAAGTCGGATGCGCTCGGCGTGGCCTATTATAATCAACACGCTATGTTGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACAGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATa                                                                    <  2:1076295/150‑1 (MQ=255)
                                                                                gTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCCAGCCAAATTCATTGATATAC‑TGCCAGGTCATa                                                                    >  1:259820/1‑150 (MQ=255)
                                                                                  gATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAAt                                                                  >  1:504591/1‑150 (MQ=255)
                                                                                   cttaggcGGATGCGCTCTGCGAGCCCTATTATAATAAACACTCTATTTAGTTTGTTCTGGCCACGAAATCTCGGTGATTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATTCATTGATATAC‑AGCCCGGTCATAATg                                                                 <  1:288725/144‑1 (MQ=255)
                                                                                   aTAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATg                                                                 <  2:1235971/150‑1 (MQ=255)
                                                                                   aTAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCCACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATg                                                                 >  2:1569845/1‑150 (MQ=255)
                                                                                   aTAAGTCGGATGAGCTCGGCGTGCAATATTATAATCAACAAGCTATGTAGTTTGTTATGGCAACGACATCTATGTGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATg                                                                 <  1:1157746/150‑1 (MQ=255)
                                                                                    tAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCGCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGa                                                                >  1:1638306/1‑150 (MQ=255)
                                                                                     aaGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAg                                                               <  2:1014017/150‑1 (MQ=255)
                                                                                      aGTCGGATGCGCTCGGCGTGCCCTATTATATTCAACAGGATATGTAGTTTGTTATGCCCCCGACTTCTCGGGGCTTATTAACTTCCCCCCTTTCCCGCTTTACGCCACCGCAAGCCAAATAC‑TTTATTTACAACCCCGGGCATAATGAgc                                                              >  2:648541/1‑150 (MQ=255)
                                                                                      aGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAgc                                                              <  1:1686338/150‑1 (MQ=255)
                                                                                         cGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑ATCCCGGTCATAATGAgcacc                                                           >  2:1363575/1‑150 (MQ=255)
                                                                                         cGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAgcacc                                                           <  1:927869/150‑1 (MQ=255)
                                                                                               cgcTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCt                                                     <  2:1398069/150‑1 (MQ=255)
                                                                                               cgcTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCt                                                     >  1:911598/1‑150 (MQ=255)
                                                                                                     gCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAAt                                               >  1:637318/1‑150 (MQ=255)
                                                                                                         gCCCTATTATAATCACCACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCTCCGCACCTAAAAATTGCa                                           >  1:1721712/1‑150 (MQ=255)
                                                                                                           ccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGa                                         <  1:1367351/150‑1 (MQ=255)
                                                                                                              gttatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTAAGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGGCATAATGAGCACCGCACCTAAAAATTGCAGAccc                                      >  1:1198395/2‑150 (MQ=255)
                                                                                                                 ataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGtt                                   <  2:1194377/150‑1 (MQ=255)
                                                                                                                  taATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGTCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTa                                  >  2:1363597/1‑150 (MQ=255)
                                                                                                                   aaTCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGGTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGCCCCGTTaa                                 >  2:1246747/1‑150 (MQ=255)
                                                                                                                      cAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATAAAC‑AGCCCGGTCATAATGCGCACCGCACCTAATAATTGCAGACCCGTTAAGCg                              >  2:227275/1‑150 (MQ=255)
                                                                                                                        acacGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGtt                            <  1:33800/150‑1 (MQ=255)
                                                                                                                        acacGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGtt                            <  1:917242/150‑1 (MQ=255)
                                                                                                                         cacGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTc                           <  2:740520/150‑1 (MQ=255)
                                                                                                                          acGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTGAAGCGTTCa                          >  2:881227/1‑150 (MQ=255)
                                                                                                                            gCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATc                        <  2:911598/150‑1 (MQ=255)
                                                                                                                             cTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATcc                       <  2:1398444/150‑1 (MQ=255)
                                                                                                                             cTATGTAGTTTGTTATGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCTCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATcc                       <  2:315648/150‑1 (MQ=255)
                                                                                                                              tATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCca                      <  2:1157492/150‑1 (MQ=255)
                                                                                                                                tGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTCTACGCCACCGCAATCGAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGCCCCGGTTAGCGGTCATCcaac                    >  2:295349/1‑150 (MQ=255)
                                                                                                                                 gTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCcaaca                   <  1:1122343/150‑1 (MQ=255)
                                                                                                                                 gTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCcaaca                   >  1:588882/1‑150 (MQ=255)
                                                                                                                                     tttGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAg               <  2:637318/150‑1 (MQ=255)
                                                                                                                                        taataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCCCCGCCAGCCAAATACATTGAGATAG‑AGCCCGGGCATAAGGAGCACCGCACCTAAAAATTGCAGACCCGTAAAGCGTTCATCCACCAATAGTGc            >  2:842219/7‑150 (MQ=255)
                                                                                                                                        taataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAATTTGCAGACCCGTTAAGCGTTCATCCAACAATAGTga            >  2:315984/7‑149 (MQ=255)
                                                                                                                                        gTTCTGGCCCCGACATCTCTGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGAACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGc            <  1:324027/150‑1 (MQ=255)
                                                                                                                                         aataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGGCATAATGAGCACCGCACCAAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGcc           >  2:665545/6‑150 (MQ=255)
                                                                                                                                              gCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACt      <  1:1037455/150‑1 (MQ=255)
                                                                                                                                                 ccGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGTCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGc   <  1:1363597/150‑1 (MQ=255)
                                                                                                                                                  cGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGcc  <  1:1523515/150‑1 (MQ=255)
                                                                                                                                                  cGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGcc  >  2:308329/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TACTTCTTGTGCCAATGTCGTGGTTGTTCCCCCATCGCAGGAACACCCACCGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAC‑AGCCCGGTCATAATGAGCACCGCACCTAAAAATTGCAGACCCGTTAAGCGTTCATCCAACAATAGTGCCGCACTTGCC  >  NC_000913/1620075‑1620370

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