Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,806,897 C→A R296R (CGC→CGA)  proW → glycine betaine ABC transporter membrane subunit ProW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,806,8970CA91.1% 119.5 / ‑3.0 49R296R (CGC→CGA) proWglycine betaine ABC transporter membrane subunit ProW
Reads supporting (aligned to +/- strand):  ref base C (2/1);  major base A (4/37);  minor base G (3/1);  total (9/40)
Fisher's exact test for biased strand distribution p-value = 9.16e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCTCTGCCGCCGATTATCCGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATC‑TGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGCGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCG  >  NC_000913/2806659‑2807136
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
ttGCTCTGCCGCCGATTATCCGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATC‑TGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGgtc                                                                                                                                                                                                                                       <  1:277953/249‑1 (MQ=255)
  gCTCTGCCGCCGATTATCCGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATC‑TGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGgtcgtc                                                                                                                                                                                                                                    <  1:625459/250‑1 (MQ=255)
        ccggggaggctCCGTATGCCCATTCTGGGGATTAACCAGGGTCCGGTGGATCGTAATT‑AAGTCTCGGGCTCATTCGGGGCTAGCCCGCGCCAGATGCTGTTAAAAGTTCAGTTACCTCTTGCGATACCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATTGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATTGGCCGTCTGGata                                                                                                                                                                                                                              <  2:157963/240‑1 (MQ=255)
         cgccgATTATCCGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATC‑TGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATg                                                                                                                                                                                                                            <  1:541624/251‑1 (MQ=255)
               gtATCCGGCTGCCCATCCGGGGGATTAACCAGGTTCCGGCGGCTC‑TGATTGAAGCCGCGCGCTCATTCGGGGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGTTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACCT‑CGCGGTATCGGTCGTCTGGATATGGGGcc                                                                                                                                                                                                                       <  2:83936/249‑2 (MQ=255)
                tATCCGCCTGACCATTCTGGGGATTAACCGGGTTCCGGCGGATC‑TGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCTCGCCAGATGCTGGTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCGGGCCCTTTCTATGGTGGT‑CATCGCTTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCtt                                                                                                                                                                                                                      <  2:428650/250‑1 (MQ=255)
                    cgtcggcccATTCTGGGGATTAACCAGGTTCCGGCGGATC‑TGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTCCCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGctac                                                                                                                                                                                                                 <  2:486004/244‑4 (MQ=255)
                                  gtggtttacccAGGTTCCGGCTGTTC‑TGATTGTAGCCTCGCGCTCATTAGGTTCCAGCCCGCGCCAGATGCTGTTCAAAGGTCAGTTACCTCTGGCGATGCCGTCCATTATGGCTGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGTTGGT‑CGTCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGTGGCTTCCCCCCGTTGGCTGCGTc                                                                                                                                                                                                    <  2:138578/241‑1 (MQ=255)
                                        tAACCAGGTTCCGGCGGATC‑TGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGAtt                                                                                                                                                                                              <  1:52865/250‑1 (MQ=255)
                                         aaCCAGGTTCCGGCGGATC‑TGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGTGATGCCGTCCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTG‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATtg                                                                                                                                                                                             <  1:503329/250‑1 (MQ=255)
                                         aaCCAGGTTCCGGCGGATC‑TGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATtg                                                                                                                                                                                             <  1:649413/250‑1 (MQ=255)
                                         aaCCAGGTTCCGGCGGATC‑TGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATtg                                                                                                                                                                                             <  1:532319/250‑1 (MQ=255)
                                                  ccGGCGGATC‑TGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCg                                                                                                                                                                                    <  1:677166/250‑1 (MQ=255)
                                                              gATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTAt                                                                                                                                                                             <  1:54332/246‑1 (MQ=255)
                                                                           cgcacagttggTGCCAGCCCGCGCAAGATGCTGTTCAAAGGTCAGTTACCGCTTGCGATGCCTACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTAGCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCa                                                                                                                                                           <  1:381630/242‑1 (MQ=255)
                                                                            ccccagtCGGGGCCAGACCGCGCCAGAGGCTTTTCAAAGTTCAGTGCCTGCTGCCGTTTCCGACCATTATGGCGGGCGTTA‑ACCGTACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGGCTGACGCAg                                                                                                                                                          <  2:184131/245‑1 (MQ=255)
                                                                               cATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAgg                                                                                                                                                         <  2:53805/249‑1 (MQ=255)
                                                                                    ggTGCCAGCCCGCGCCCGATGCTGTTCAAAGTTCCGTTACCGCTGGGGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCGGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGAACGTCTGACGCAGGCCGTTgg                                                                                                                                                  <  2:289151/251‑1 (MQ=255)
                                                                                     gtgccgtccctcgcCATATGCTTTTCAAAGTTCATTTACCGCTGGCGCTGCCGACCATTATGGCTGGATTTA‑ACCAGACGCGGATGCTGGCCCTTTCTATGGTGGT‑CCTCGCCTCGATGATTGCCTTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCCGTCGTCTGGATATGGGGCTTGCCCCCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTggg                                                                                                                                                 <  2:555877/240‑1 (MQ=255)
                                                                                       gCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGgcg                                                                                                                                               <  1:398046/251‑1 (MQ=255)
                                                                                       gCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGgcg                                                                                                                                               <  1:472244/251‑1 (MQ=255)
                                                                                              cgcgCCGGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCGCGATGATTGCCGGCGGCGGGTTGGGTCAGATGGTACTT‑CGAGTTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGAGCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCa                                                                                                                                        <  1:286719/251‑1 (MQ=255)
                                                                                                  ccataggcggTTCAAAGTTCAGTTTCCTGTGGCTATGCCGACCATTAGGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCATTTCTATGGTGGT‑CATCGCATCGATGATTGCCATCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACgta                                                                                                                                    <  2:206248/242‑3 (MQ=255)
                                                                                                      aTGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGGTGGGTCCGATGGGACTTTCAGGGT‑TCGGGCGTCTGGGGATGGGGGTTTCCACCCGTGGGGGCGTCGGGGGTTTGGTCCTCGCCCTTTTCCTCCATCCTCTGACGCAGGCCGGTggggggggccccccccggcg                                                                                                                                >  1:661375/1‑234 (MQ=255)
                                                                                                              cAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAAcc                                                                                                                        <  1:574020/251‑1 (MQ=255)
                                                                                                                  gTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCg                                                                                                                    <  1:450132/251‑1 (MQ=255)
                                                                                                                     catttctcGCTGGCGGTGCCGACCATTATGGCGGGAGTTA‑ACCAGACGCTGAGGCTGGCCCTGTCTTGGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGAGGGTACTT‑CGAGGTATCGGTCGTCGGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCTCCATTATCCTCGATCGTCGGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGTCAATCGTCGCt                                                                                                                  <  2:79455/243‑1 (MQ=255)
                                                                                                                             gCTGGCGAGGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTAcacca                                                                                                         <  1:318734/251‑1 (MQ=255)
                                                                                                                                ggCGATGCCGACCATTATGGCGGGCGTTA‑ACCCGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGACTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGTTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGTGGTTGGATTTTGATCCTCTCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTg                                                                                                      <  1:631777/251‑1 (MQ=255)
                                                                                                                                ggCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTTATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCTGCGGGTTGGGTCAGATGGTACTT‑CGAGGTAGCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGACGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACt                                                                                                       <  1:50992/250‑1 (MQ=255)
                                                                                                                                               gttgtcgggcgTGT‑ACCAGACTCTGATGGCGGCCCTTGCGATGGGGGG‑AATCGCCTCGATTATGGCCGTCGGCGGGGTGGGTCAGAGGGTACTC‑CGAGGTATCGGTCGTCTGGATAGGGGGCTTGCCACCGTTGGCGTCGTCGGGATGGTGATCTTCGCCAGTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGGCGCTGGTACGCCACTGGCCCTGTTGGCctg                                                                                        <  2:34596/245‑1 (MQ=255)
                                                                                                                                                ttGGCGGG‑GTTACACCAGACTCGGGTGCTGGCCCTTTCTATGGTGGT‑CAGCGCCTCGATGATTGCCGTCGGCTGGGTGGTTCAGATGGTACTT‑CGAGGTTTCGGTCGTCTGGATAGGGGGCTTGCCACGGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGGCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTctgc                                                                                       <  2:499324/249‑1 (MQ=255)
                                                                                                                                                        cgTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGGGGGATCGGTCGTCTGGATCTGGGGCTTGCCACCGTTGGCGGCGTCGGGCTTGTGATCCTCGCCATTATCCTCGATCGGCTGACGCAGGCCGTTGGGCGCGACTCACACAGTCGCGGCAACCGACGCTGGTCCACCACTGGCCCTGTTGGTCTGCTGCCCCGcc                                                                              >  1:617763/1‑251 (MQ=255)
                                                                                                                                                              acccgACGCTGATGCTGGCCCTTTCTATGGTGTT‑CATCGCCTCGATGATTGCCGGCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGTGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCcattc                                                                         <  1:397794/247‑1 (MQ=255)
                                                                                                                                                                        gATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACCCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAAc                                                               >  1:20729/1‑251 (MQ=255)
                                                                                                                                                                           gCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGCGGTATCGGTCGTCTGGGTATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCCGGCCGTTGGGCGCGGCTCCCCCCGTCGCGGCCACCGTCGCTGGGACACCACTGGCCCCGTTGGTCTGCCGACACGCCCCTTCATTAAGTAActct                                                            >  2:578144/1‑251 (MQ=255)
                                                                                                                                                                             tGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTTGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCCCCGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAActcttc                                                          >  1:316248/1‑249 (MQ=255)
                                                                                                                                                                              ggCCCTTTCTATG‑TGGTGCAGCGCTTCGTTGCTTGCCTTCGGCGGGGTGGGTCCGACGGTTCTT‑TGTGGTATCGTGCGTCTGGATAGGGGTCTTGCCACCGTTGGCGGGGGCGGGATTGTGAGCCTCGCCATTATCCTCGAGCGTCTGACGCAGTCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTATACCACTGGTCCTGTTGGTCTGCTTTCCCGCCCATTCATTAAGTAACGCTGc                                                          <  2:437377/250‑1 (MQ=255)
                                                                                                                                                                                    ttCTATGGTGGT‑CAGCGCCGAGATGATGGCCGTCGGCGGGTGGGGTCAGAGGGTGCTT‑CGAGGGATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGGTTTCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTg                                                     <  2:523883/249‑1 (MQ=255)
                                                                                                                                                                                          ggtggt‑CATCGCCTCGATGATTGCCGTCGGCTGGTTGGGTCAGATGGTACTT‑CGCGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCca                                                                            >  1:26734/1‑220 (MQ=255)
                                                                                                                                                                                               t‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCCGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGaacggctg                                                                                                                                                                >  1:497449/1‑124 (MQ=255)
                                                                                                                                                                                               t‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTg                                                                                                                                                                <  2:497449/131‑1 (MQ=255)
                                                                                                                                                                                                 tctCGCCTCGGTGATTGCCGTCGTCGGGTTGGGTCAGCTGGTACTT‑CGGGGTATCGGTCGGTTGGAGATGGGGCTGGCCACCGTTGGCCGCGTCGTGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGGCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCg                                        <  2:58150/248‑1 (MQ=255)
                                                                                                                                                                                                         cGATGATTGCCGGCGGCGGGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATc                               <  1:646283/251‑1 (MQ=255)
                                                                                                                                                                                                                 gCCGTCGGCGGGTTGGGTCAGATGGTACTA‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGACCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGc                        <  1:316307/250‑1 (MQ=255)
                                                                                                                                                                                                                      cggcggGTTGGGTCCAATGGTACTT‑CGGGGGATTGGGCGTCTGGGTATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCCTTTTCCTCGATCGTATGGCGCCGGCCGTTTGGGGCGGCTCCCTCCGTCGCCGCAACCGTCGCTGGGTCACCCCTGGCCCCTTTTGTCTCTTCACCCGCCCCTTTATTTAGGTACTTTCCACTTGCCCCGTGACGCCCGGCCTTATCcccccctccaaa                   >  2:160872/1‑242 (MQ=255)
                                                                                                                                                                                                                        gcggGTTGGGTCAGATGGTACTT‑CGAGGTCTCGGTCGTCTGGTTATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCaaaaaaa                <  1:123722/251‑1 (MQ=255)
                                                                                                                                                                                                                         cggGTTGGGTCAGATGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAg               <  1:578144/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aTGGTACTT‑CGAGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGCCTGCTGACCAGCCCATTCATTAAGTAAATCTGCACCTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAACAAGGGATAACAATGCg  >  2:158926/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
TTGCTCTGCCGCCGATTATCCGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATC‑TGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTA‑ACCAGACGCTGATGCTGGCCCTTTCTATGGTGGT‑CATCGCCTCGATGATTGCCGTCGGCGGGTTGGGTCAGATGGTACTT‑CGCGGTATCGGTCGTCTGGATATGGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCG  >  NC_000913/2806659‑2807136

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