Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,059,948 T→A 5.2% intergenic (+22/+16) asnU → / ← cbl tRNA‑Asn/ssuEADCB/tauABCD operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,059,9480TA 5.2% 290.2 / 4.5 97intergenic (+22/+16)asnU/cbltRNA‑Asn/ssuEADCB/tauABCD operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (48/44);  new base A (2/3);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-01

AATCGGGTTTGACAAAAGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCT  >  NC_000913/2059799‑2060091
                                                                                                                                                     |                                                                                                                                               
aaTCGGGTTTGACAAAAGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCtttt                                                                                                                                                 >  1:1537796/1‑150 (MQ=255)
 aTCGGGTTTGACAAATGATTTTTCGCCGTTAAGATGTGCCTCAACAACTATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTta                                                                                                                                                <  1:810947/150‑1 (MQ=255)
     ggTTTGACAAAAGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAATGCCCGCTTTTATAGc                                                                                                                                            >  2:60235/1‑150 (MQ=255)
       tttGACAAAAGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCgg                                                                                                                                          <  1:8029/150‑1 (MQ=255)
          gACAAAAGATTTTTCGCCGTTAAGATGGGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCACGCTTTTATAGCGGGAt                                                                                                                                       >  2:1280532/1‑150 (MQ=255)
             aaaaGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGAtttt                                                                                                                                    >  1:1187765/1‑150 (MQ=255)
              aaaGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGAttttt                                                                                                                                   >  2:1487397/1‑150 (MQ=255)
              aaaGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTATAAAAGCGGGCttttc                                                                                                                                   <  1:704217/150‑2 (MQ=255)
               aaGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCTGACGGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAATCCCGCTATAAAAGCGGGCttttca                                                                                                                                  <  2:121382/150‑3 (MQ=21)
               aaGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTg                                                                                                                                  <  2:526094/150‑1 (MQ=255)
               aaGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTGCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTg                                                                                                                                  <  1:229575/150‑1 (MQ=255)
               aaGATTTTTCGCAGTTAAGATGTGCCACCACCACGATTCCTCTGTAGTTCAGTCGGTAGAACGCCGTACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAATCCCGCTATAAAAGCGGGCttttca                                                                                                                                  <  2:1204986/150‑3 (MQ=17)
                cgATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGc                                                                                                                                 <  2:140675/149‑1 (MQ=255)
                aGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGc                                                                                                                                 <  2:1288487/150‑1 (MQ=255)
                        cGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTg                                                                                                                         <  1:634067/150‑1 (MQ=255)
                         gCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATATCTGGATTTTTGCTATATCTGa                                                                                                                        >  1:626519/1‑150 (MQ=255)
                          ccGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGAt                                                                                                                       >  1:1357274/1‑150 (MQ=255)
                          ccGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGAt                                                                                                                       <  1:50786/150‑1 (MQ=255)
                           cGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGCGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGAGTTTTGCTTTATCTGATa                                                                                                                      >  2:271026/1‑150 (MQ=255)
                           cGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATa                                                                                                                      <  1:446158/150‑1 (MQ=255)
                            gTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATaa                                                                                                                     <  2:997013/150‑1 (MQ=255)
                                  aTGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGGTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCGTTTATAGCGGGATTTTTGCTATATCTGATAATCAAtt                                                                                                               <  1:49169/150‑1 (MQ=255)
                                   tgtgCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAAttt                                                                                                              >  1:1470394/1‑150 (MQ=255)
                                     tgCCTCAACAACGATTACTATGTAGTTCAGTCTTTAGAACGGAGGACTGTTAATCCGTATGGCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTcc                                                                                                            <  2:125072/150‑1 (MQ=25)
                                      gCCTCAACTACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGAGTGTTAATCCGTATGTATCTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCct                                                                                                           <  2:1654890/150‑1 (MQ=255)
                                       ccTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCctc                                                                                                          >  2:46516/1‑150 (MQ=255)
                                        cTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCctct                                                                                                         >  1:1694993/1‑150 (MQ=255)
                                        cTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTACTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCctct                                                                                                         >  1:1701849/1‑150 (MQ=255)
                                          caacaaCGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTc                                                                                                       <  1:990631/150‑1 (MQ=255)
                                          caacaaCGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACCGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCTGTCTTTTTGCTATATCTGATAATCAATTTCCTCTTc                                                                                                       >  1:76704/1‑150 (MQ=255)
                                             caacGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACt                                                                                                    >  2:1237598/1‑150 (MQ=255)
                                             caacGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACt                                                                                                    >  1:1798071/1‑150 (MQ=255)
                                             caacGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACt                                                                                                    >  1:1257887/1‑150 (MQ=255)
                                             caacGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACt                                                                                                    >  2:784491/1‑150 (MQ=255)
                                               acGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGc                                                                                                  <  1:451296/150‑1 (MQ=255)
                                                cGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCt                                                                                                 >  1:28372/1‑150 (MQ=255)
                                                     ccTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTAAATCCGTATGTCACTGGTTAGAGTCCAGTCAGAGGAGCCAAATTCCTGAATAGCCCGCTTTTATTGCGGGATGTTTGCTATATCTGATAATCAATTTCCTCGTCACTGCTTTCCa                                                                                            >  2:1601956/1‑150 (MQ=25)
                                                       tctGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATc                                                                                          <  1:1186163/150‑1 (MQ=255)
                                                        ctGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCGAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCa                                                                                         >  1:958158/1‑150 (MQ=255)
                                                          gTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCAcc                                                                                       <  1:1231055/150‑1 (MQ=255)
                                                            aGTTCAGTCGGTAGATCGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGAATATTGCTATAACTGTAAATCAATTTCCTCTTCACTGCGTTCCATCACCTg                                                                                     >  1:889273/1‑150 (MQ=25)
                                                              ttCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGcc                                                                                   <  1:512712/150‑1 (MQ=255)
                                                                cAGTCGGAAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTAGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTACACTGCTTTCCATCACCTGCCGc                                                                                 >  2:1276303/1‑150 (MQ=255)
                                                                  gTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCtt                                                                               <  2:1740661/150‑1 (MQ=255)
                                                                  gTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCtt                                                                               >  1:1525031/1‑150 (MQ=255)
                                                                   gcGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTg                                                                              <  2:314348/149‑1 (MQ=255)
                                                                    cGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGa                                                                             <  2:228203/150‑1 (MQ=255)
                                                                    cGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGa                                                                             <  2:548753/150‑1 (MQ=255)
                                                                     ggTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGat                                                                            <  2:841382/150‑1 (MQ=255)
                                                                     ggTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGat                                                                            <  2:704207/150‑1 (MQ=255)
                                                                        aGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATc                                                                         >  1:603666/1‑150 (MQ=255)
                                                                         taaCGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATcc                                                                        <  1:1201321/149‑1 (MQ=255)
                                                                         gAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATcc                                                                        >  2:893158/1‑150 (MQ=255)
                                                                         gAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATcc                                                                        >  2:670298/1‑150 (MQ=255)
                                                                          aaCGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCt                                                                       >  1:1766804/1‑150 (MQ=255)
                                                                            cggcggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCa                                                                     >  2:505697/1‑150 (MQ=255)
                                                                              gcggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAAc                                                                   <  1:1311725/150‑1 (MQ=255)
                                                                                ggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTg                                                                 >  1:1544964/1‑150 (MQ=255)
                                                                                  aCTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGAc                                                               <  2:603124/150‑1 (MQ=255)
                                                                                     gTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGCCAGt                                                            >  1:393609/1‑150 (MQ=255)
                                                                                       tAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTcc                                                          >  1:743218/1‑150 (MQ=255)
                                                                                        aaTCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCt                                                         >  1:31450/1‑150 (MQ=255)
                                                                                        aaTCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCt                                                         >  1:1448079/1‑150 (MQ=255)
                                                                                         aTCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTg                                                        <  1:721086/150‑1 (MQ=255)
                                                                                         aTCCGTATGTCACTGGATCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTg                                                        <  1:471226/150‑1 (MQ=255)
                                                                                          tCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGc                                                       >  1:1525623/1‑150 (MQ=255)
                                                                                               aTGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCAGTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTAc                                                  <  2:1694993/150‑1 (MQ=255)
                                                                                               aTGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTAc                                                  <  2:1701849/150‑1 (MQ=255)
                                                                                                tGTCGCTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACa                                                 <  1:158963/150‑1 (MQ=255)
                                                                                                   cACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTTTATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAg                                              >  2:1452474/1‑150 (MQ=255)
                                                                                                    aCTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGt                                             >  1:258229/1‑150 (MQ=255)
                                                                                                     cTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGtt                                            <  2:1470394/150‑1 (MQ=255)
                                                                                                      tGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTc                                           <  1:1673523/150‑1 (MQ=255)
                                                                                                      tGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTc                                           <  2:485867/150‑1 (MQ=255)
                                                                                                        gTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCa                                         >  1:1425788/1‑150 (MQ=255)
                                                                                                         ttCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAg                                        <  1:1487397/150‑1 (MQ=255)
                                                                                                             aGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAg                                    >  2:1044202/1‑150 (MQ=255)
                                                                                                               tCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAgcg                                  <  2:1741107/150‑1 (MQ=255)
                                                                                                                 cAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGcc                                >  2:1260378/1‑150 (MQ=255)
                                                                                                                    tCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGa                             >  1:1807395/1‑150 (MQ=255)
                                                                                                                    tCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGa                             <  2:28372/150‑1 (MQ=255)
                                                                                                                     cAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGAc                            >  1:924000/1‑150 (MQ=255)
                                                                                                                      agagGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACa                           >  1:427934/1‑150 (MQ=255)
                                                                                                                          gagCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGt                       >  1:699365/1‑150 (MQ=255)
                                                                                                                            gCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTa                     <  2:529357/150‑1 (MQ=255)
                                                                                                                            gCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTa                     <  1:1128467/150‑1 (MQ=255)
                                                                                                                            gCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTa                     >  2:162929/1‑150 (MQ=255)
                                                                                                                            gCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTa                     >  1:1575004/1‑150 (MQ=255)
                                                                                                                               aaaTTCCTGAAAAGCGCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGc                  <  1:347206/150‑1 (MQ=255)
                                                                                                                               aaaTTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGc                  >  2:414409/1‑150 (MQ=255)
                                                                                                                               aaaTTACTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGc                  <  1:1645515/150‑1 (MQ=255)
                                                                                                                                 aTTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTg                >  2:580843/1‑150 (MQ=255)
                                                                                                                                    ccTGAAAAGCCCGCTTTTATAGCGGGATTTTTGTTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAg             <  1:780313/150‑1 (MQ=255)
                                                                                                                                    ccTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAg             <  2:1357741/150‑1 (MQ=255)
                                                                                                                                     cTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGt            >  1:59779/1‑150 (MQ=255)
                                                                                                                                          aatcccGCTATAAAAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTcc       >  2:1628770/4‑150 (MQ=255)
                                                                                                                                          aaGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTcc       <  2:258229/150‑1 (MQ=255)
                                                                                                                                           atcccGCTATAAAAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCt      >  1:546976/3‑150 (MQ=255)
                                                                                                                                               cGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCt  >  1:568019/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
AATCGGGTTTGACAAAAGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCT  >  NC_000913/2059799‑2060091

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: