Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 657,377 T→C V29A (GTG→GCG)  cspE → constitutive cold shock family transcription antitermination protein; negative regulator of cspA transcription; RNA melting protein; ssDNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913657,3770TC98.5% 241.7 / NA 67V29A (GTG→GCG) cspEconstitutive cold shock family transcription antitermination protein; negative regulator of cspA transcription; RNA melting protein; ssDNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (28/38);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAA  >  NC_000913/657231‑657523
                                                                                                                                                  |                                                                                                                                                  
aaTTTGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGtt                                                                                                                                                 <  1:684396/150‑1 (MQ=255)
aaTTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTc                                                                                                                                                >  1:894184/1‑151 (MQ=255)
 aTTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCg                                                                                                                                               <  1:522195/151‑1 (MQ=255)
   tGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGt                                                                                                                                              <  2:1045228/150‑1 (MQ=255)
    ggTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTa                                                                                                                                             <  1:493086/150‑1 (MQ=255)
         aCTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTtc                                                                                                                                       >  1:330389/1‑151 (MQ=255)
            ggTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTtctc                                                                                                                                     <  2:365463/150‑1 (MQ=255)
             gTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTg                                                                                                                                   >  1:116246/1‑151 (MQ=255)
                     cacaGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAg                                                                                                                           <  2:262486/151‑1 (MQ=255)
                      acaGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATcc                                                                                                                             >  2:193772/1‑148 (MQ=255)
                        aGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTATGCAATCCAGACt                                                                                                                        <  2:1191603/151‑1 (MQ=255)
                           gtttggggcaaTTTTTCATGAAGAGAAAATTTGGATTGCCCAGTTTAAAGGTACCGGTAAGCGGTTTAATGAATCAAAAGATTACGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTTGAACCCTTCGCTGCAATCCAGACTAAt                                                                                                                     >  1:1026783/11‑151 (MQ=255)
                              tgtgTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTTCGGTAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATgg                                                                                                                   >  2:711725/1‑150 (MQ=255)
                               gtgtCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGt                                                                                                                  <  1:410641/150‑1 (MQ=255)
                                      tttttCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaaa                                                                                                           <  2:848761/150‑1 (MQ=255)
                                       ttttCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAAc                                                                                                          <  1:998719/150‑1 (MQ=255)
                                       ttttCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAAc                                                                                                          <  2:832085/150‑1 (MQ=255)
                                           cATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTg                                                                                                     >  2:756450/1‑151 (MQ=255)
                                            aTGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGc                                                                                                    >  1:1001335/1‑151 (MQ=255)
                                            aTGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGc                                                                                                    >  2:776289/1‑151 (MQ=255)
                                               tAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTg                                                                                                  >  1:274940/1‑150 (MQ=255)
                                                aaaGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGaa                                                                                                <  2:1109227/151‑1 (MQ=255)
                                                    gTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTc                                                                                            <  1:927898/151‑1 (MQ=255)
                                                     tAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTc                                                                                            >  1:27882/1‑150 (MQ=255)
                                                      aaTTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTc                                                                                            >  2:991165/1‑149 (MQ=255)
                                                       aTTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgc                                                                                         >  1:263106/1‑151 (MQ=255)
                                                          ttGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcgc                                                                                       <  1:1082418/150‑1 (MQ=255)
                                                          ttGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCC‑AAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGt                                                                                     >  2:635780/1‑151 (MQ=255)
                                                            gATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGt                                                                                     >  2:1070089/1‑150 (MQ=255)
                                                             aTGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTag                                                                                   <  2:479818/151‑1 (MQ=255)
                                                             aTGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTa                                                                                    <  1:415056/150‑1 (MQ=255)
                                                               gTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTagag                                                                                 <  2:198228/151‑1 (MQ=255)
                                                                   aaGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCTATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTc                                                                              <  2:511302/150‑1 (MQ=255)
                                                                   aaGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGtt                                                                               <  1:821033/149‑1 (MQ=255)
                                                                   aaGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGt                                                                                <  1:979412/148‑1 (MQ=255)
                                                                    aGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCg                                                                             >  2:858519/1‑150 (MQ=255)
                                                                         aaaGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATc                                                                        <  1:172448/150‑1 (MQ=255)
                                                                               aaCGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACg                                                                 >  1:571996/1‑151 (MQ=255)
                                                                                aCGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACgg                                                                >  1:2601/1‑151 (MQ=255)
                                                                                   ttAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGc                                                             >  2:176037/1‑151 (MQ=255)
                                                                                      aGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCa                                                           >  2:929634/1‑150 (MQ=255)
                                                                                         ggTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAgg                                                       >  2:575808/1‑151 (MQ=255)
                                                                                             tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCt                                                   >  1:541889/1‑151 (MQ=255)
                                                                                              aaTGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCt                                                   <  1:809833/150‑1 (MQ=255)
                                                                                                tGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTct                                                >  1:621638/1‑151 (MQ=255)
                                                                                                 gAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTctg                                               <  2:924717/151‑1 (MQ=255)
                                                                                                  aGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTctgc                                              <  2:948586/151‑1 (MQ=255)
                                                                                                     ccAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTctgctg                                            <  2:526516/150‑1 (MQ=255)
                                                                                                     ccAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCtgctgc                                           <  1:130356/151‑1 (MQ=255)
                                                                                                     ccAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCtgctgc                                           >  1:264313/1‑151 (MQ=255)
                                                                                                       aaaGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCa                                          >  1:843316/1‑150 (MQ=255)
                                                                                                         aGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAAc                                       <  2:572137/151‑1 (MQ=255)
                                                                                                                   ttCATTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGctc                             >  2:1080188/1‑151 (MQ=255)
                                                                                                                      aTTACTCCGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGt                          <  1:1022912/151‑1 (MQ=255)
                                                                                                                            ccGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGAt                     <  2:257764/150‑1 (MQ=255)
                                                                                                                            ccGGAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGAt                     <  1:756278/150‑1 (MQ=255)
                                                                                                                              ggAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGccc                                                    <  1:1001065/117‑1 (MQ=255)
                                                                                                                              ggAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGccc                                                    >  2:1001298/1‑117 (MQ=255)
                                                                                                                               gAAGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGt                 >  1:57751/1‑151 (MQ=255)
                                                                                                                                 aGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCa               <  1:281273/151‑1 (MQ=255)
                                                                                                                                 aGACGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCa               <  2:434680/151‑1 (MQ=255)
                                                                                                                                   aCGGCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGc             <  2:356127/151‑1 (MQ=255)
                                                                                                                                     ggCAGCAAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCaa           <  2:862602/151‑1 (MQ=255)
                                                                                                                                      gcagcaAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAg          <  1:193729/151‑1 (MQ=255)
                                                                                                                                         gcaAAGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAAt       <  2:717649/151‑1 (MQ=255)
                                                                                                                                           aaaGACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTc     >  1:570920/1‑151 (MQ=255)
                                                                                                                                              gACGCGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCaaa  <  2:149838/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AATTGGTGTACTGGTAACCGACACAGCATTTGTGTCTATTTTTCATGTAAAGGTAATTTTGATGTCTAAGATTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAAGATACGTCAGCAAGAATTCAAA  >  NC_000913/657231‑657523

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑