Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 357,491 G→A 5.9% intergenic (+37/+300) codA → / ← cynR cytosine/isoguanine deaminase/transcriptional activator of cyn operon; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,4910GA 5.9% 386.1 / 7.2 118intergenic (+37/+300)codA/cynRcytosine/isoguanine deaminase/transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base G (54/57);  new base A (4/3);  total (58/60)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

Only 100 of 121 total aligned reads displayed.

GTCAGGTTCCGGTACGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACG  >  NC_000913/357342‑357634
                                                                                                                                                     |                                                                                                                                               
gTCAGGTTCCGGTACGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATg                                                                                                                                                 <  2:140882/150‑1 (MQ=255)
  cAGGTTCCGGTACGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAACCGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATgcg                                                                                                                                               >  2:666335/1‑150 (MQ=255)
    ggTTCCGGTACGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATgcggc                                                                                                                                             >  2:1641946/1‑150 (MQ=255)
            tACGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAACCCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                                     <  2:786343/150‑1 (MQ=255)
            tACGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGc                                                                                                                                     >  1:707564/1‑150 (MQ=255)
              cGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCt                                                                                                                                   <  2:1064585/150‑1 (MQ=255)
              cGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTGGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCt                                                                                                                                   <  2:1521955/150‑1 (MQ=255)
               gTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                  >  2:356879/1‑150 (MQ=255)
                ttATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                 >  1:33270/1‑150 (MQ=255)
                   ttCGGTACGTGGCGGCAAGGGGATTGCCAGCACACAACCGGCACAAACCACAGTATATCTGGAGCATCCAGAAGCCACCGATTACAAACGTTGAACGTCTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCATATcc                                                                                                                              <  2:147581/150‑1 (MQ=255)
                       gTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                          <  2:179572/150‑1 (MQ=255)
                       gTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                                          <  2:1120981/150‑2 (MQ=255)
                            tGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGCGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGt                                                                                                                     >  1:818289/1‑150 (MQ=255)
                            tGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGCATaaac                                                                                                                     <  1:1101469/150‑4 (MQ=38)
                             ggcggcAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTg                                                                                                                                          >  2:488704/1‑128 (MQ=255)
                             ggcggcAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTa                                                                                                                    >  1:530817/1‑150 (MQ=255)
                                    aGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACt                                                                                                             >  1:741756/1‑150 (MQ=255)
                                    aGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACt                                                                                                             >  1:994051/1‑150 (MQ=255)
                                    aGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACt                                                                                                             <  1:1063711/150‑1 (MQ=255)
                                     ggTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGACATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTg                                                                                                            <  1:398762/150‑1 (MQ=255)
                                       tGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGaa                                                                                                          >  1:795626/1‑150 (MQ=255)
                                        gtgtGCCATCACACAACCGGCACAAACCCCCGTATATCTGGAGCAGCCAGTAGCCATCGATTACACCCGTCGAACCACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGATCACTGAAc                                                                                                         <  1:123752/147‑1 (MQ=37)
                                        gATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCTGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAAc                                                                                                         <  1:991043/150‑1 (MQ=255)
                                        gATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAAc                                                                                                         <  2:818289/150‑1 (MQ=255)
                                        gATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAAc                                                                                                         <  1:356879/150‑1 (MQ=255)
                                         aTTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACt                                                                                                        <  1:1237036/150‑1 (MQ=255)
                                            gCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGt                                                                                                     >  2:1201228/1‑150 (MQ=255)
                                             ccAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTa                                                                                                    >  1:1483590/1‑150 (MQ=255)
                                              cAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAg                                                                                                   <  1:922332/150‑1 (MQ=255)
                                              cAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAg                                                                                                   >  2:363209/1‑150 (MQ=255)
                                                      aaCCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATa                                                                                           <  1:942911/150‑1 (MQ=255)
                                                       aCCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATaa                                                                                          <  1:363209/150‑1 (MQ=255)
                                                       aCCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATaa                                                                                          >  1:1061562/1‑150 (MQ=255)
                                                           gCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGt                                                                                      <  2:1022557/150‑1 (MQ=255)
                                                           gCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGt                                                                                      >  2:47933/1‑150 (MQ=255)
                                                             acaAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAg                                                                                    <  1:613260/150‑1 (MQ=255)
                                                              caAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAgc                                                                                   >  1:1670669/1‑150 (MQ=255)
                                                                aaCCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAgcgc                                                                                 <  2:1076960/150‑1 (MQ=255)
                                                                aaCCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAgcgc                                                                                 >  1:1403652/1‑150 (MQ=255)
                                                                  ccaccGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCAt                                                                               <  2:907466/150‑1 (MQ=255)
                                                                    accGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCa                                                                             <  1:1402819/150‑1 (MQ=255)
                                                                        tataTCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                         >  2:1704618/1‑150 (MQ=255)
                                                                        tataTCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                         <  2:33270/150‑1 (MQ=255)
                                                                          tatCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAAt                                                                       >  1:1576500/1‑150 (MQ=255)
                                                                           atCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAAtt                                                                      >  2:8562/1‑150 (MQ=255)
                                                                           atCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAAtt                                                                      >  1:833935/1‑150 (MQ=255)
                                                                            tCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTc                                                                     <  2:684084/150‑1 (MQ=255)
                                                                             cTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTcc                                                                    <  2:121607/150‑1 (MQ=255)
                                                                              tGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACTCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCa                                                                   >  1:580997/1‑150 (MQ=255)
                                                                              tGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCa                                                                   >  1:1371571/1‑150 (MQ=255)
                                                                              tGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCa                                                                   >  2:504/1‑150 (MQ=255)
                                                                              tGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCa                                                                   >  1:674817/1‑150 (MQ=255)
                                                                                 cgcagcCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGAGGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGcc                                                                <  1:524172/149‑1 (MQ=255)
                                                                                    cgcCAGAAGCCTTCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGCTTATGCCGGCTGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCt                                                             <  1:666335/149‑1 (MQ=255)
                                                                                    agcCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCt                                                             >  1:1302704/1‑150 (MQ=255)
                                                                                     gcCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTg                                                            >  1:7311/1‑150 (MQ=255)
                                                                                     gcCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTg                                                            >  1:1615707/1‑150 (MQ=255)
                                                                                     gcCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCATGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTg                                                            >  1:523917/1‑150 (MQ=255)
                                                                                       cAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGAt                                                          >  2:794780/1‑150 (MQ=255)
                                                                                          aaGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCtg                                                       <  2:693839/150‑1 (MQ=255)
                                                                                           aGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCtgt                                                      >  1:1137050/1‑150 (MQ=255)
                                                                                                tCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAg                                                 <  2:707564/150‑1 (MQ=255)
                                                                                                 cGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGATTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCAATGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCAGCAGCGCTGTGTCAGc                                                >  1:71145/1‑150 (MQ=17)
                                                                                                  gATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCg                                               <  2:795626/150‑1 (MQ=255)
                                                                                                  gATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCg                                               <  1:1595984/150‑1 (MQ=255)
                                                                                                  gATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCg                                               >  1:64477/1‑150 (MQ=255)
                                                                                                    ttACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGGCAGCGGc                                             <  1:641551/150‑1 (MQ=255)
                                                                                                     tACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGATCTGTGTCAGCGGCt                                            >  2:468452/1‑150 (MQ=255)
                                                                                                      aCAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTa                                           <  1:1641946/150‑1 (MQ=255)
                                                                                                       cAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTAc                                          <  2:1576500/150‑1 (MQ=255)
                                                                                                       cAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTAc                                          >  1:1482543/1‑150 (MQ=255)
                                                                                                       cAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTAc                                          >  1:838566/1‑150 (MQ=255)
                                                                                                         aaCGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCg                                        >  1:491307/1‑150 (MQ=255)
                                                                                                          aCGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGt                                       <  1:281964/150‑1 (MQ=255)
                                                                                                             ttGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCATGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGat                                    <  2:523917/150‑1 (MQ=255)
                                                                                                              tGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGctt                                   >  2:533018/1‑150 (MQ=255)
                                                                                                               gAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGattc                                  <  2:574923/150‑1 (MQ=255)
                                                                                                               gAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGattc                                  >  1:1586345/1‑150 (MQ=255)
                                                                                                               gAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGattc                                  <  2:442796/150‑1 (MQ=255)
                                                                                                                 aCGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGattcat                                <  2:1483590/150‑1 (MQ=255)
                                                                                                                  cGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGattcatt                               >  1:800401/1‑150 (MQ=255)
                                                                                                                     cTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTccc                            <  1:1331339/150‑1 (MQ=255)
                                                                                                                     cTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTccc                            >  1:415970/1‑150 (MQ=255)
                                                                                                                      tGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCg                           <  2:610557/150‑1 (MQ=255)
                                                                                                                       gggTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTGATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGc                          <  2:1061562/150‑1 (MQ=255)
                                                                                                                              aGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCcaacacc                   >  1:1442150/1‑148 (MQ=255)
                                                                                                                               gCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGAACTGTGTCAGCGGCTACCGTGATTCATTACCGCCAACAAcc                  >  2:590992/1‑150 (MQ=255)
                                                                                                                                cGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACcg                 >  2:1204096/1‑150 (MQ=255)
                                                                                                                                   gCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCgcgc              <  2:86631/150‑1 (MQ=39)
                                                                                                                                   gCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGGGGCAGCGGCTACCGGGATTCCGTCACGCCAACAACCgcgc              >  1:1206045/1‑150 (MQ=14)
                                                                                                                                         tttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTACAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATtcct        <  1:129399/150‑1 (MQ=255)
                                                                                                                                           tATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATtcctcc      <  2:7311/150‑1 (MQ=40)
                                                                                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCaa    >  2:1617923/2‑150 (MQ=14)
                                                                                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAAc   >  2:395943/1‑150 (MQ=14)
                                                                                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAAc   >  2:402092/1‑150 (MQ=14)
                                                                                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAAc   >  1:967617/1‑150 (MQ=14)
                                                                                                                                               ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACg  >  2:1227500/1‑150 (MQ=255)
                                                                                                                                               ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACg  <  2:1615707/150‑1 (MQ=255)
                                                                                                                                               ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACg  >  2:84719/1‑150 (MQ=255)
                                                                                                                                               ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACg  <  1:1054316/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GTCAGGTTCCGGTACGTTATTCGGTACGTGGCGGCAAGGTGATTGCCAGCACACAACCGGCACAAACCACCGTATATCTGGAGCAGCCAGAAGCCATCGATTACAAACGTTGAACGACTGGGTTACAGCGAGCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACG  >  NC_000913/357342‑357634

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