Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,273,229 A→T 5.1% intergenic (+27/+17) chaC → / ← ychN cation transport regulator/putative sulfur relay protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,273,2290AT 5.1% 272.6 / 4.3 79intergenic (+27/+17)chaC/ychNcation transport regulator/putative sulfur relay protein
Reads supporting (aligned to +/- strand):  ref base A (33/42);  new base T (2/2);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

TGGAACAGGAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTC  >  NC_000913/1273081‑1273375
                                                                                                                                                    |                                                                                                                                                  
tGGAACAGGAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCa                                                                                                                                                    >  1:54527/1‑149 (MQ=255)
     cAGGAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGc                                                                                                                                               <  2:108032/149‑1 (MQ=255)
        gAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCCCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGcccc                                                                                                                                            >  1:412956/1‑149 (MQ=255)
            ttATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCa                                                                                                                                        >  1:121667/1‑149 (MQ=255)
             tATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGccccactat                                                                                                                                       <  2:508496/149‑6 (MQ=255)
               tCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAAc                                                                                                                                     <  3:295019/149‑1 (MQ=255)
                    cTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGcccca                                                                                                                                           >  1:710526/1‑137 (MQ=255)
                    cTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCTGACGcccca                                                                                                                                           <  2:710526/138‑2 (MQ=255)
                     tGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTaaa                                                                                                                               >  1:659950/1‑149 (MQ=255)
                       ggAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAAt                                                                                                                             <  2:163573/149‑1 (MQ=255)
                       ggAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAAt                                                                                                                             >  4:481352/1‑149 (MQ=255)
                         aaTGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAAtgt                                                                                                                           <  3:786246/149‑1 (MQ=255)
                         aaTGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAAtgt                                                                                                                           <  3:785824/149‑1 (MQ=255)
                         aaTGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAAtgt                                                                                                                           <  2:315952/149‑1 (MQ=255)
                         aaTGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAAtgt                                                                                                                           <  1:506228/149‑1 (MQ=255)
                         aaTGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGCATGTTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAAagt                                                                                                                           <  2:897259/149‑1 (MQ=255)
                           tGCAGGCCGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAACCTACTGGCGAAGAATTTTCCGGATGGTGTGTCACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTAAGACCCCCCCTCAAACATTAAAATGTGa                                                                                                                         <  2:169492/149‑1 (MQ=255)
                            gCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAg                                                                                                                        >  2:890652/1‑149 (MQ=255)
                                    aTGGGCTGAATGATTTGCTGGTATCGGTAAAAACACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTAt                                                                                                                <  4:413938/149‑1 (MQ=255)
                                     cgatcTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATc                                                                                                               <  1:258429/145‑1 (MQ=255)
                                       ggCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCgg                                                                                                             <  2:509041/149‑1 (MQ=255)
                                        gCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGc                                                                                                            <  1:511700/149‑1 (MQ=255)
                                         cTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCt                                                                                                           <  4:448394/149‑1 (MQ=255)
                                         cTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCt                                                                                                           <  1:890652/149‑1 (MQ=255)
                                          tGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTg                                                                                                          >  2:190769/1‑149 (MQ=255)
                                            aaTGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGAc                                                                                                        >  3:555233/1‑149 (MQ=255)
                                            aaTGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGAc                                                                                                        >  2:162718/1‑149 (MQ=255)
                                                aTTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCg                                                                                                    <  1:770697/149‑1 (MQ=255)
                                                 tttGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGt                                                                                                   >  2:192218/1‑149 (MQ=255)
                                                  ttGCTGGTATCGGTAAAAAAACTGCTGGCGGAGTATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTc                                                                                                  <  4:583797/149‑1 (MQ=255)
                                                          aTCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTgcgc                                                                                          <  4:506335/149‑1 (MQ=255)
                                                           tCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGcc                                                                                         >  1:582376/1‑149 (MQ=255)
                                                             ggTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAg                                                                                       >  1:179750/1‑149 (MQ=255)
                                                               tAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATGGGCTGACAGCGTCCATTGCGCCAGtt                                                                                     >  2:379085/1‑149 (MQ=255)
                                                                aaaaaaaCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTc                                                                                    >  2:542832/1‑149 (MQ=255)
                                                                   aaaaCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCcac                                                                                 <  3:747185/149‑1 (MQ=255)
                                                                     aaCTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACaa                                                                               <  1:379085/149‑1 (MQ=255)
                                                                     aaCTGCTGGCGGAGAATTTGCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCACCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACaa                                                                               <  4:604318/149‑1 (MQ=255)
                                                                       ctgctgGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCTCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAaga                                                                             <  2:26772/149‑1 (MQ=255)
                                                                       ctgctgGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAaga                                                                             <  3:58675/149‑1 (MQ=255)
                                                                        tgctgGCGGAGAATTTTCCTGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAagag                                                                            >  3:517802/1‑149 (MQ=255)
                                                                        tgctgGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGt                                                                                      >  1:646697/1‑139 (MQ=255)
                                                                        tgctgGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGt                                                                                      <  2:646697/139‑1 (MQ=255)
                                                                        tgctgGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAagag                                                                            <  4:278458/149‑1 (MQ=255)
                                                                          ctgGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTa                                                                          <  3:481352/149‑1 (MQ=255)
                                                                               ggAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGtt                                                                                     <  3:716636/133‑1 (MQ=255)
                                                                               ggAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGtt                                                                                     >  4:716636/1‑133 (MQ=255)
                                                                                gagaATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGAtt                                                                    >  2:382808/1‑149 (MQ=255)
                                                                                gagaATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGAtt                                                                    >  3:402826/1‑149 (MQ=255)
                                                                                  gaATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTc                                                                  <  2:589199/149‑1 (MQ=255)
                                                                                      tttCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCAcc                                                              <  4:593776/149‑1 (MQ=255)
                                                                                          cGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCat                                                          <  1:537717/149‑1 (MQ=255)
                                                                                            gATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCt                                                                                                           >  4:581408/1‑98 (MQ=255)
                                                                                            gATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCt                                                                                                           <  3:581408/98‑1 (MQ=255)
                                                                                            gATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCatca                                                        >  1:177708/1‑149 (MQ=255)
                                                                                               ggTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCatcaatc                                                     >  1:75085/1‑149 (MQ=255)
                                                                                                  gtgtTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCaga                                                  >  4:229103/1‑149 (MQ=255)
                                                                                                   tgtTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCagag                                                 <  3:300423/149‑1 (MQ=255)
                                                                                                      tACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGaa                                              <  3:217052/149‑1 (MQ=255)
                                                                                                         gTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGa                                                                                                         <  1:361158/87‑1 (MQ=255)
                                                                                                         gTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGa                                                                                                         >  2:361158/1‑87 (MQ=255)
                                                                                                             gggATTCGCCTTGGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACt                                       >  4:123178/1‑149 (MQ=255)
                                                                                                                      cTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAgtg                                           <  2:724793/136‑1 (MQ=255)
                                                                                                                      cTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAgtg                                           >  1:724793/1‑136 (MQ=255)
                                                                                                                      cTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCcgcg                              >  2:374439/1‑149 (MQ=255)
                                                                                                                        gAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGcc                            >  1:51740/1‑149 (MQ=255)
                                                                                                                          gTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTc                                                                                    <  3:118751/91‑1 (MQ=255)
                                                                                                                          gTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTc                                                                                    >  4:118751/1‑91 (MQ=255)
                                                                                                                           tAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGt                         >  2:280526/1‑149 (MQ=255)
                                                                                                                            aaaCTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTc                        >  2:284213/1‑149 (MQ=255)
                                                                                                                            aaaCTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTc                        >  2:69834/1‑149 (MQ=255)
                                                                                                                             aaCTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCg                       >  4:428157/1‑149 (MQ=255)
                                                                                                                             aaCTTCCCGCATAGTGGGGCGACAGACGCGCCCTCAAACATCAAAATGTGAGCACTTTATCGCCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAAACAGAGGAAGTGTACTAATCCCGCGCCCGTCg                       <  3:543954/149‑1 (MQ=255)
                                                                                                                                  cccGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACa                  <  2:177708/149‑1 (MQ=255)
                                                                                                                                   ccGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAg                 <  2:65520/149‑1 (MQ=255)
                                                                                                                                    cGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAgg                >  2:181656/1‑149 (MQ=255)
                                                                                                                                        tAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGtttt            <  4:517069/149‑1 (MQ=255)
                                                                                                                                            ggggCGTCTGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGcac        >  2:781281/1‑149 (MQ=255)
                                                                                                                                               gCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCAGACAGCATCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGGACAGGTTTTGCACAAt     <  4:539842/149‑1 (MQ=255)
                                                                                                                                                  tCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTc  <  2:121667/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGGAACAGGAGCTTATCAAACTGGGAATGCAGGACGATGGGCTGAATGATTTGCTGGTATCGGTAAAAAAACTGCTGGCGGAGAATTTTCCGGATGGTGTGTTACGTCCGGGATTCGCCTGAGTAAACTTCCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTGTACTAATCCCGCGCCCGTCGGTACAGGTTTTGCACAATTTC  >  NC_000913/1273081‑1273375

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