Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 71,877 A→G 24.9% E176G (GAG→GGG)  yabI → DedA family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091371,8770AG24.9% 59.7 / 3.6 37E176G (GAG→GGG) yabIDedA family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (4/22);  new base G (7/3);  total (12/25)
Fisher's exact test for biased strand distribution p-value = 3.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGAT‑C‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCCG‑CGGTCGTTTGTTGTGGCTG‑ACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAA  >  NC_000913/71660‑72099
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                 
gCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGTATG‑CTGGATCTGCCGGTGGCCCAATTTATTCCGCCGACTA‑TTTACGGCTGCCTGCTGTGG‑CCGCCGTTTTTCTTCTTGCCAGGGATTTTGGCGGGCGCGGCCGT‑CGGATTTTCTTCCCGGAAT‑CAGAGCGGTGGGGTTTaaag                                                                                                                                                                                                                           >  1:209057/1‑226 (MQ=255)
 cGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGGT‑C‑GCTATTCCTGCCGGAATGCCGAGCGGG‑GGGTTTAAATggc                                                                                                                                                                                                                        >  2:155575/1‑227 (MQ=255)
     gCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAAGA‑TTATCGGCTGCCTGCTGTTG‑CCGCCGTTTTACTTCCTGCCCGGGATTCTGGCGGGGGCGGCGGT‑C‑CATATTCCTGCCGGAATGCAGAGGGGT‑GGGGGTAAATggggggg                                                                                                                                                                                                                    >  1:293514/1‑223 (MQ=255)
                gCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCCATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCCAAATTTATTACCCCGAATA‑TTATCCGCTGCCTGCTGTTG‑CCGCCGGTTTACTTCCCTTCCGGGGTTCTGCCGGGCGCGGGGATCC‑GACA‑TCCGTCCGGAATGCCGAGAGGG‑GAGttttacagggtggcgctggcccacg                                                                                                                                                                                                         >  1:432211/1‑206 (MQ=255)
                gCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGGTTTTGGCGGGGGCGGCGGC‑C‑GATATTTCTTCCCGCAAGAGGAGCGGG‑GGGGTTAAAGGGTGGTGGGTGGCAAAAg                                                                                                                                                                                                         >  1:304889/1‑228 (MQ=255)
                  ttGTTCCCCGTTTTGGGGGGTTTTTTTGTTGGCCGGGGGGGTCGGCTGGTGCCATGGGTGGCGGGAATG‑CTGGACCTCCGGGGGGTTGAGTTTATTCCGCCGAATA‑TTACCGGGTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGAGTCTGGCGGGCGCGGCGAT‑C‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCg                                                                                                                                                                                                       <  1:325177/227‑1 (MQ=255)
                              cTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGCTTCTGGCGGGCGCGGGGAT‑C‑GATATTCCTGCCGGAATGCAGGGCGGG‑TGGTTTGAATGGTTGCTGCTGGCAACAGCGGGGGTTTTgggg                                                                                                                                                                                           >  2:104262/1‑228 (MQ=255)
                               tGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTTTGG‑GCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGGGAT‑C‑GCAATTCCTTCCGGCATGCAGAGCCGG‑GCGGTTACATGGTGGCTGCTGGTAACAGCGGGGtgttttgggg                                                                                                                                                                                          >  1:189722/1‑223 (MQ=255)
                                 gtcggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTTGCGGGCGCGGGGGT‑C‑CGTTATCCTCCCCGACAGCAGCGGGGT‑GGGGTTAAAGGGGTGTCGCGTGCAACACCGGGGttttttgggggg                                                                                                                                                                                        >  1:379966/1‑221 (MQ=255)
                                    ggtcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGAT‑C‑GATATTCCTGCCGGAATGCAGAGCGGG‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGGGTTTTTGGGGGTtggt                                                                                                                                                                                     >  2:448242/1‑228 (MQ=255)
                                      tcgTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTTGATCTGCCGGTGGCTAAATTTATTACGCCGGATA‑TTATCGGCTGCCCGCTGTGG‑CCCCCGTTTTTATTCCTGCCCGGGATTCTGGCGGGCCCGGGGAT‑T‑GCTATTTCTGCCCGAATCGCGAACGGG‑GGGTTTTAATGGTTGTTTCATGCACCCACGGGGTTTTTTTGGGTtggggg                                                                                                                                                                                   >  1:175273/1‑228 (MQ=255)
                                                  cccGACGGGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGAT‑C‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCt                                                                                                                                                                        <  2:432211/227‑1 (MQ=255)
                                                              gtgtgggccatgggtggGGGGAATG‑GTGGATCGGCCGGTGGCTAATTTTATTACGCCGAATA‑TTGTCGGCTGCCGGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGAT‑C‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgcg                                                                                                                                                           <  2:175273/216‑1 (MQ=255)
                                                              gtgggTGCAAATGGGGGCGGGAAGG‑CTGGATCGGCCGGTGGCTAAATTTATTACGCCGAATA‑TTGTCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGAT‑C‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgcg                                                                                                                                                           <  1:488268/225‑1 (MQ=255)
                                                              ggtGGTGCC‑ATGGTGGGGGGTATGTCTGGTCCTGCGGGTGCGAAAATTTTTTA‑GCCCAATATTTATGGGGTGCC‑GTTGTGGCCCGCGGTTTTTCTTCCGGCCAGGGATTCTGGCGGGCGCGGCGAT‑C‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGgcg                                                                                                                                                           <  1:345358/226‑1 (MQ=255)
                                                                                        tgtataagggccagggCTAAATTTATTCCGCCGAATA‑TTATCGGGTGCCTGCGGGGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGAT‑C‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGGTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGttt                                                                                                                                  <  2:531130/213‑1 (MQ=255)
                                                                                        ggggtacTGCCGGGGGCTAAATTTTTTGCGCCGAATA‑TTATCGGCTGCCTGCTGGGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGAT‑C‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGttt                                                                                                                                  <  2:209057/221‑1 (MQ=255)
                                                                                        cTGGATCTGCCGGTGGCTAAATTTATTACGGCGAATA‑TTAGCGGCTGCCTGGTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGAT‑C‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGttt                                                                                                                                  <  2:452875/227‑1 (MQ=255)
                                                                                                         tatttttATTACGCAGAATT‑GTAGCGGCTGCCTGCTGGGG‑CCGCAGTTTTACTTTCTGCCAGTGGTTCTGGCGGGCGCGGCGAT‑C‑GATACTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGATGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCcg                                                                                                                <  1:330623/224‑1 (MQ=255)
                                                                                                           aatttcagccggcGAATT‑TTATCGGCTGACGGGTGGGG‑CCGCCGTTTTACTGCCTGCCCGGGATGCTGGCGGGCGCGGCGAT‑C‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGGTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCcg‑cg                                                                                                             <  2:555668/216‑1 (MQ=255)
                                                                                                            aTTTATTACGCCGAATT‑TTACCGGCTGCCGGCGGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGAT‑C‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCcg‑cg                                                                                                             <  2:341547/227‑1 (MQ=255)
                                                                                                                 ttaCGCCGAATA‑TTAGCGGCTGCCTGCTGTGG‑CCGGCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGAT‑C‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCCG‑CGGTCGt                                                                                                        <  2:586567/227‑1 (MQ=255)
                                                                                                                                        cctgcTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGAT‑C‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTTTGGGGCAGCGGGTAA‑GCGGCGGACCGGTTTAGTCATTATTTTTGCCC‑CGGGGGGTTGTTGTGGGTG‑GCGCCCtttgt                                                                                >  2:520821/1‑225 (MQ=255)
                                                                                                                                                           ttACCTCCTGCCGGGGATTCTGGCGGGCGCGGCGAT‑C‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCCG‑CGGTCGTTTGTTGTGGCTG‑ACGCCGTTGATTTCTGCCATCGGCgtggt                                                              <  2:201832/228‑1 (MQ=255)
                                                                                                                                                            ttctgccggCCAGGGATTTTGGGGGGCGGGGGTAT‑C‑GATTGTCCTGCCGGAATGTAGAGGGGG‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTGGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCCG‑CGGTCGTTTGTTGTGGCTG‑ACGCCGTTGATTTCTGCCATCGGCgtggt                                                              <  1:380458/219‑1 (MQ=255)
                                                                                                                                                                  cTGCCAGGGATTCTGGCGGGCGCGGCGAT‑C‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCCG‑CGGTCGTTTGTTGTGGCTG‑ACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCtg                                                       <  2:806/228‑1 (MQ=255)
                                                                                                                                                                       aGGGATTTTTGGTGGCGCGGGGGG‑T‑GATTTTCCGGCCGGAATGCAGGGGGGT‑GAGGTTAAATGGTTGCTGCTGGCAAGAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCCG‑CGGTCGTTTGTTGTGGCTG‑ACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTggtggt                                                  <  1:104262/228‑1 (MQ=255)
                                                                                                                                                                                        cggcgAT‑C‑GATATTGCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCCG‑CGGTCGTTTGTTGTGGCTG‑ACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGt                                 <  2:379966/228‑1 (MQ=255)
                                                                                                                                                                                          cggtgatggtttttCCTGCCGGAATGCAGAGCGGT‑TGGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCCG‑CGGTCGTGTGTTGTGGCTG‑ACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGt                                 <  1:566251/216‑1 (MQ=255)
                                                                                                                                                                                              t‑c‑gatATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTTGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTAATTGTCCCG‑CGGTCGTTTGTTGTGGCTG‑ACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATg                            <  2:428020/227‑1 (MQ=255)
                                                                                                                                                                                                        ccTGCGGGAATGGAGGGTGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTTGTTTGGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCCG‑CGGTCGTTTGTTGTGGCTG‑ACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGta                    <  2:608665/227‑1 (MQ=255)
                                                                                                                                                                                                          tGCCGGAATGCAGGGCGGT‑GATTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCCG‑CGGTCGTTTGTTGTGGCTG‑ACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGtata                  <  2:304889/227‑1 (MQ=255)
                                                                                                                                                                                                            ccGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGGGGTGTTTTTGTGGGTTGGGTGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCCG‑CGGTCGTTTGTTGTGGCTG‑ACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCg               <  2:189722/228‑1 (MQ=255)
                                                                                                                                                                                                                 aTGGAGAGTGGT‑GGGTTTAAATGGTTGCGTCGGGCAACAGGGGGTTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCCG‑CGGTCGTTTGTTGTGGCTG‑ACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATAtt          <  1:493559/228‑1 (MQ=255)
                                                                                                                                                                                                                      ggggggT‑GGGTTTAAATGGTTGGTGCTGGCAACAGGGGGGTTTTTGGGGGTTGGTGGCTGGTTGTGCGGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCCG‑CGGTCGTTTGTTGTGGCTG‑ACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCg     <  1:590194/224‑1 (MQ=255)
                                                                                                                                                                                                                       agCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTCTTTGGGCCGCGGGAAAGGCGTGGGACGGTTT‑AGTCCTTTTTTTTTCCGCCGCGGGTTGGTGGTGGCTGCGCGCGGGGGTTTTCT‑CCTCTGGCGCGGGGGGGGGGGTGGGGTTAATTTGCCACCCCTTGATGTCgggtgttagagcgtgtttgtg     >  1:185722/1‑209 (MQ=255)
                                                                                                                                                                                                                         cGGT‑GAGTTTAAATGGGTGCTGCTGGCAACAGCGGTGGTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCCG‑CGGTCGTTTGTTGTGGCTG‑ACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTaa  <  2:293514/228‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                 
GCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATG‑CTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATA‑TTATCGGCTGCCTGCTGTGG‑CCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGAT‑C‑GATATTCCTGCCGGAATGCAGAGCGGT‑GAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAA‑GCGACTGACCGTTTGAGTCATTATTTGTCCCG‑CGGTCGTTTGTTGTGGCTG‑ACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAA  >  NC_000913/71660‑72099

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