Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,629,382 T→C 23.4% S53G (AGC→GGC)  rbbA ← ribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,629,3820TC23.4% 91.5 / 2.6 46S53G (AGC→GGC) rbbAribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein
Reads supporting (aligned to +/- strand):  ref base T (16/16);  new base C (5/7);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACACCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCA‑GGCCACGGAAGTTCTTCATTCACCCGCACCCACGCT  >  NC_000913/3629158‑3629601
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                             
tCGCACTTCCCGCCCCGCTTTGGCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACCCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCg                                                                                                                                                                                                                             <  1:71932/227‑1 (MQ=255)
 cGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACACCGACAAGGTGTGGTAGAGGGTTTTGCCCAGCGCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGAGGCTTCGGGTCGCGCATATCGCCGCCC‑AGCCCCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAAAGACAACAACCCCGa                                                                                                                                                                                                                            >  1:74782/1‑227 (MQ=255)
              ccgctTTGTCGTGACCGAACAGGCGATCGAAAAAATCGACGTTTTCAAACACCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGt                                                                                                                                                                                                              <  1:327993/228‑1 (MQ=255)
              ccgctTTGTCGTGACCGAACAGGCGAGACAAAATATCGGCGTTTTCATACACCGACAAGGTGTGGTAGAGGTTTGTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCGACGACAACAAGCTCGACTTCCCGACGCCGt                                                                                                                                                                                                              <  1:63074/228‑1 (MQ=255)
                    tGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACACCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAAGGACAACAAGCCGGACTTCCCGACGCCGTCCggg                                                                                                                                                                                                         >  1:209306/1‑227 (MQ=255)
                           aCCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACACCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACCTTGCCCTGTTCAATGACCCGGGCACCGGAAATCAACGACAACAAGCGCGACTTCCCCACGCCGTCCGGGCCACTCa                                                                                                                                                                                                  >  1:273045/1‑227 (MQ=255)
                           aCCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACACCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACCCGGGCACCGCAAATCAACGACAACAAGCCCGACTTCCCGACGCCGTCCGGGCCAATCAg                                                                                                                                                                                                 >  1:312340/1‑228 (MQ=255)
                              gAACGGGCGAGCGAAAAACTCGACGTTTTCATACCCCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGccc                                                                                                                                                                                              <  1:459002/228‑1 (MQ=255)
                                         cgaAAAAATCGACGTTTTCATACACCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCATGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCTCCC‑AGCACCATCACATTGCCCTGTTCACTGACGCGGGCACCCGAAATCAACGACAACAAGCCCCACTTCCCGCCGCCGTCCGGGCCAATCAGCCCGACCAAACacc                                                                                                                                                                                   >  1:623102/1‑226 (MQ=255)
                                                        tttctaacacCGAAAAGGTGTGGTGGAGGTTTTTGCCCAGCCCCTGCGGGATACAGGCGATGCGCGGGCAGACGTCGGGGCGGTGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCga                                                                                                                                                                    <  2:137736/222‑1 (MQ=255)
                                                           catcccccggaaaGGTGGGGTAGGGGTTTTTGCCACGCCCCGGCGGCCTCCAGGGCATGCGCGGGGAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGagag                                                                                                                                                                 <  2:22108/217‑1 (MQ=255)
                                                           cagacacCGACAAGGTGTGGTGGGGGTTGTGGCCCAGCCCCTGGGGCAGCCAGGCGAGGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGAGCATATCGCCGGCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCTCCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGagag                                                                                                                                                                 <  2:199278/225‑1 (MQ=255)
                                                           cATCCACCGACAAGGTGTGGTAGAGTTTTTTCCCCAGCCCCTGCGGCGTCCGGGCGATGCGCGGGCATACGTCGCGGCGATGCGTCGGGTCGCGCATATCGCCGCAC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGagag                                                                                                                                                                 <  2:209306/228‑1 (MQ=255)
                                                           cATACCCCGGAAGGGTGGGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGAGGCGCGGGCAGACGTTGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGagag                                                                                                                                                                 <  2:565615/228‑1 (MQ=255)
                                                                  cGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGCCAACAAGCTCACCTTCCCGCCGCCGTCCGGGCCAATCAGCCCGCCCATAACACGGGCCGGAATATCGACAGGGatat                                                                                                                                                           >  1:587951/1‑227 (MQ=255)
                                                                     cAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGtt                                                                                                                                                       <  1:447418/228‑1 (MQ=255)
                                                                          tgtgGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAgc                                                                                                                                                   >  1:391578/1‑227 (MQ=255)
                                                                          tgtgGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAgc                                                                                                                                                   >  1:217333/1‑227 (MQ=255)
                                                                            ggggagagGTTTTTCCCCAGCCCCGGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑Ca                                                                                                                                               <  1:441336/224‑1 (MQ=255)
                                                                              gTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAAc                                                                                                                                             <  1:218721/228‑1 (MQ=255)
                                                                                    gTTTTTGCCCCGCCCCTTCGGCATCCAGGCGATGCCGGGGCAGACGTCGCGGCGATTCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGttt                                                                                                                                        <  2:317038/227‑1 (MQ=255)
                                                                                    gTTTTTGCCCAGCCCCTGCGGCATCCAGGCGCTGCGCGGGCAGACGTCGCGGCGCGGCTTCGGGTCGCGCATATCGCTGCCC‑AGCACCATCCCATTGCCCTGCTCAATGGCGCGGGCACCGGCAATCAACGCCAACAAGCCCGCCTTCCCGCCGCCGTCCGGGCCAATAAACCCGACACTACACCCGGCCGGAACACCGTCAGGGATATT‑TTCAACGCCAACGGtttt                                                                                                                                       >  2:162183/1‑228 (MQ=255)
                                                                                    gTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACACCAAGCTCGACTTCCCGACCCCGTCCGGGCCAATCAGCCCGACCACACGACGGGCCGGAATATCGGCAGGGACATTGTTCACCC‑CACCGGtttt                                                                                                                                       >  2:441336/1‑228 (MQ=255)
                                                                                     tttttGCCCAGCCCCGGAGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGttt                                                                                                                                        <  1:28776/226‑1 (MQ=255)
                                                                                     tttgtGCCCAGCCCCTGCGGCAGCCCGGCGATGCCCGGGGCGGAGTCGCGGCGATTCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTTGTCAATGACGCGGGCCCCGGAAATCAACGACAACAAGCCCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGttt                                                                                                                                        <  2:54860/226‑1 (MQ=255)
                                                                                       tttGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGACCGACCATACATCGGGCCGGAATATCGAGAGTGATATTGGTCACCG‑CAACGGTTTTTc                                                                                                                                     >  1:294141/1‑227 (MQ=255)
                                                                                               gCCCCTGCGGCATCCAGGCGATGCGCGGGCAGCCGTCCCGGAGATGCTTCGGGGCGCGCATAGCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGc                                                                                                                             <  1:401858/227‑1 (MQ=255)
                                                                                                     gCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACACGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAACGGGCCGGGATATCGAGAGGGGTATTGTTCAGGG‑CCACGGCTTTTCCCTaaaccgggcaca                                                                                                                      >  2:222417/1‑218 (MQ=255)
                                                                                                          aTCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAACCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCCTAATGCTGGCTCAcgccc                                                                                                                 >  1:336635/1‑228 (MQ=255)
                                                                                                                     gcgcGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGCAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGCATATCGGGAGGGATATTGTTCAAGC‑CAACGGTTTTTCCATAAAGCTGGCTCCCGCCCGCCAGtttgcg                                                                                                      >  2:517407/1‑224 (MQ=255)
                                                                                                                            ccgccGTCGCGGCGGTGTTGCGGTGCGGCCAAATCGCCGCCC‑AGCCCCATCCCATGGCCCTTTTCAATGCCGCGGGCACCGGAAATCAACGACAACAAGCGCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAgg                                                                                                <  2:312340/223‑1 (MQ=255)
                                                                                                                              gcctgcgcgGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCCCCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcg                                                                                              <  1:329262/222‑1 (MQ=255)
                                                                                                                              gACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCAAAATGCTGGCCCACGCCCCCCAGGTCCGCGACAAGcgg                                                                                             >  2:269422/1‑228 (MQ=255)
                                                                                                                                        ggttGCTTCGGGGCGCGCCATCCGCCGCCCAACCACCA‑CCCTTTGCCCTGTTCAAGGACGCGGGCACCGGAAATCAACGACAACAAGCCCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAAc                                                                                    <  2:587951/224‑1 (MQ=255)
                                                                                                                                                                cccccc‑acccccaccACATTGCCCGTTTCAAGCAGCGGGGCCCCGGAAATCAACGACACCAAGCCGGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGGGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCa                                                           <  2:336635/225‑1 (MQ=255)
                                                                                                                                                                              tCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGCCAGGCGGGACGGGAACCAGTTGCCGATGCGGCCTTGCGGCAACCTCACCACCAg                                              >  1:500330/1‑228 (MQ=255)
                                                                                                                                                                                   ttGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGggcgt                                         >  1:208528/1‑228 (MQ=255)
                                                                                                                                                                                          gTTCAATGACGCGGGCCCGGGAAAACAACGAAAACAAGCCCGTCTTCCCGACGCCGTCCGGGCCAAACAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCA‑GGCCa                                 <  1:222417/228‑1 (MQ=255)
                                                                                                                                                                                           ttCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCCGAGGGGCCATTGCGGCAACCTCACCACGAGGGCGGCA‑ggga                                  >  1:338161/1‑224 (MQ=255)
                                                                                                                                                                                             ccaggACCCGGGCACCGGAAACCAACGGCAACAAGCCCGACTTCCCGACGCGTGCCGGTCCAAGCAGCCCGACCAGACAGCGGGCCGGAATATCGAGAGGGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCA‑GGCCACgg                              <  2:217333/224‑1 (MQ=255)
                                                                                                                                                                                                tGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCCGTTCCCGATGCGTCATTGCGCGAACCCCACCCCGAGGGCGTCA‑GGcccacgca                            >  1:285443/1‑221 (MQ=255)
                                                                                                                                                                                                tGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACCGGGACCAGTTCCAGATGCGTCCATGCGGCAACCCCCCCCCGAGGGCGTCA‑GGCCACGGaa                            >  1:565404/1‑227 (MQ=255)
                                                                                                                                                                                                     cgGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCCGTTCCCGATGCGGCCATGCGGCGACCCCCCCACCAGGGGGGCAGGGCCCCGAAAGttc                        >  2:228349/1‑227 (MQ=255)
                                                                                                                                                                                                     cgGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATCGGCCCTTGCGGACACGCCACCCCGCGGCCGCCC‑GGCCACCGAAGttctt                      >  2:175227/1‑228 (MQ=255)
                                                                                                                                                                                                               aaaTCAACCACAAAAAGCCCGACTTCCCGCGGCCGGCCGGGCCAATCCGCCCGACCAAACAGCGGGCCGGAATATCGAGAGGGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCA‑GGCCACGGAAGTTCTTCATTCACCCg            <  2:208528/228‑1 (MQ=255)
                                                                                                                                                                                                                          acacaGCCCGACTTCCCGACGCCGTCCTGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCA‑GGCCACGGAAGTTCTTCATTCACCCGCACCCACGCt  <  1:228349/223‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                             
TCGCACTTCCCGCTCCGCTTTGTCGTGACCGAACAGGCGAGCGAAAAAATCGACGTTTTCATACACCGACAAGGTGTGGTAGAGGTTTTTGCCCAGCCCCTGCGGCATCCAGGCGATGCGCGGGCAGACGTCGCGGCGATGCTTCGGGTCGCGCATATCGCCGCCC‑AGCACCATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCG‑CAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCA‑GGCCACGGAAGTTCTTCATTCACCCGCACCCACGCT  >  NC_000913/3629158‑3629601

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: ‑