Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,356,201	0	T	G	46.4%	84.0 / 12.8	56	T304P (ACC→CCC)	sapA	antimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (8/22);  new base G (26/0);  total (34/22)
Fisher's exact test for biased strand distribution p-value = 3.66e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGC  >  NC_000913/1356054‑1356350                                                                                                                                                     |                                                                                                                                                        cGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGtt                                                                                                                                                      <  2:576864/151‑1 (MQ=255)       cATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGCTTTGGGGGGGTGaaa                                                                                                                                                   >  1:67518/1‑149 (MQ=255)        aTAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGTCGTACCATAATAGGGGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTCCGAGCGCATGGCGGACAGCGGGATCATTTAGCGGCGG‑CTTGGCGGGGGTaaagat                                                                                                                                                >  1:82277/1‑147 (MQ=255)        aTAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGcc                                                                                                                                               <  1:186073/151‑1 (MQ=255)             cccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑GTTGGCGGGGGTAAATGcc                                                                                                                                               >  1:695999/1‑146 (MQ=255)              ccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAg                                                                                                                                         <  1:483177/151‑1 (MQ=255)               cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAgg                                                                                                                                        <  2:490492/151‑1 (MQ=255)                aCGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGc                                                                                                                                       <  2:469733/151‑1 (MQ=255)                aCGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGc                                                                                                                                       <  1:26427/151‑1 (MQ=255)                 cGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGGTAAATGCCAGAAAAGCg                                                                                                                                      >  2:876830/1‑151 (MQ=255)                  gAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGGTAAAAGCCAGATAGGCg                                                                                                                                      >  2:970743/1‑150 (MQ=255)                          cgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCAT‑‑TAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATccc                                                                                                                          <  1:1012431/151‑1 (MQ=255)                           ggTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGGTAAATGCCAGATAAGCGAC‑GGTCATcc                                                                                                                           >  2:1017486/1‑151 (MQ=255)                                aaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTTAAAGGCAGAAAGGCGACGGTTCTTcccgg                                                                                                                        >  1:833107/1‑148 (MQ=255)                                          gccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTAAATGCCAGATAAGCGGC‑GGTCATCCCCGGGCGCAGCGtt                                                                                                             >  2:104406/1‑150 (MQ=255)                                           ccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTaa                                                                                                           <  1:557167/151‑1 (MQ=255)                                              tttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAATGACAGAAAAGCGAC‑GTTCATCCCCGGGCGCAGGGTGaaa                                                                                                          >  1:198882/1‑149 (MQ=255)                                                tCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTAAATGCCAGAAAGGCGAC‑GGTCATCCCAAGACGCAGCGTTaaa                                                                                                          >  2:110997/1‑147 (MQ=255)                                                      gTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTaaa                                                                                                          >  2:508073/1‑141 (MQ=255)                                                         ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTTAATTCCAGAAAGGCGAC‑GTTTATCCCCGGAAGCAGCGCTAAACGCAAACGCggg                                                                                              >  1:564275/1‑150 (MQ=255)                                                             aataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgtcgtc                                                                                         <  1:861904/151‑1 (MQ=255)                                                             aataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTAAAAGCCAGAAAGGCGAC‑GTTTATCCCCGGACGCAGCGTTAAAACCAAGCGCGGggcgtc                                                                                         >  1:833415/1‑151 (MQ=255)                                                                 gATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGc                                                                                     <  2:163889/151‑1 (MQ=255)                                                                         gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGa                                                                             <  2:1051387/151‑1 (MQ=255)                                                                          cATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAATGCCAGAAAGGCGAC‑GTTCATCCCAGGGCGCAGCGTTAAACGCAAGCGCGGGGCGTCACGCAAAAAgg                                                                              >  2:503128/1‑149 (MQ=255)                                                                            tCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTTAAAGCCCGATAGGCGGC‑GTTCATCCCCGGACGCAGCGTTAAAAGCAAGCGCGGGGCGGCACGCaaaa                                                                                 >  1:234081/1‑144 (MQ=255)                                                                             cAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgc                                                                         <  2:124860/151‑1 (MQ=255)                                                                               ggCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctg                                                                       <  2:134490/151‑1 (MQ=255)                                                                                 cgcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctggc                                                                     <  1:259102/151‑1 (MQ=255)                                                                                   cTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAATTCCAGATAGGCGAC‑GTTCATCCCCGGGCGCAGCGTTTAACGCAAGCGCGGGGCGTCACGCAAAATGGGTAgctg                                                                       >  2:307291/1‑147 (MQ=255)                                                                                     ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGTGTTTAAAGCCAGATAGGGGGC‑GGTCATCCCACGACGCAGCGGTAAACCCAAACGGGGGGCGGCACGCAAAAAGGATAGctggctgg                                                                  >  1:405688/1‑150 (MQ=255)                                                                                         attaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcagca                                                             <  2:163775/151‑1 (MQ=255)                                                                                           taATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGcagca‑g                                                           <  1:551098/150‑1 (MQ=255)                                                                                            aaTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTTAAAGCCACATACGCGAC‑GTTCATCCCCGGGCGCCGCGGTAAACGCAAGCGCGGGGCGGCACCCAAAAAGGATAGGTGGCTGgcag                                                               >  1:1041763/1‑146 (MQ=255)                                                                                            aaTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGca‑ggc                                                         <  2:747088/151‑1 (MQ=255)                                                                                             aTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTTAATGCCAGATAAGCGAC‑GGTCATCCCAGGACGCAGCGTTAAAAGCAAGCGCGGGGCGGCACGCAAAAAGGATAActtgctgg                                                                  >  2:432695/1‑142 (MQ=255)                                                                                              tCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGGGGGGTTAAAAGCCCGATAAGCGAC‑GTTCATCCCAGGACGCAGCGTTAAAAGCAAACGCGGGGCGTCACGCAAAATGGATAGCTGGCTGGGAGca‑ggcc                                                        >  2:360112/1‑150 (MQ=255)                                                                                                gCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGGTAAATGCCCGATAAGCGAC‑GTTTATCCCCGGGCGCCGCGGTTAACGCAAAAGCGGGGCGTCACGCAAAATGGATAGCCGGCCTGCAGcc‑ggccc                                                       >  1:720537/1‑148 (MQ=255)                                                                                                gCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTAAAAGCCAGAAAAGCCAC‑GGTCATCCCCGGACGCAGCGGTAAAAGCAAGCGCGGGGCGGCACGCAAAAAGGGTAaatggctg                                                                   >  1:372282/1‑138 (MQ=255)                                                                                                 ccAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGTTAAATGCCAGATAGGCGAC‑GTTCATCCCCGGACGCCGCGGTAAAAGCAAGCGCGGGTCGGCACGCAAAATGGATAGCTGGCTGGGAGca‑ggcccgg                                                     >  2:457453/1‑150 (MQ=255)                                                                                                 ccAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAAAGCCAGATAGGCGAC‑GTTCATCCCCGGACGCAGCGTTAAAAGCAAACGCGGGGCGGCACGCAAAATGGATAGCTGGCTGGCAGca‑ggcccgg                                                     >  2:277657/1‑150 (MQ=255)                                                                                                     tGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGTGGTAAATGCCAGAGAAGGGAC‑GGTCATCCCCAGACGCAGCGGTAAACGCAAGCGCGGGGCGGCCCGCAAAATAGATAGATGGCTGGCAGCagggccaggc                                                    >  1:638895/1‑148 (MQ=255)                                                                                                           cgcATGGCGGACAGCGGGATTATTTAGCGGCGG‑ATTTGCGGGGTTAAAAGCCACAAAGGCGAC‑GTTCATCCCAAGACGCAGCGGTAAACGCAAGCGCCGGGCGTCACGCCAAAAGGGTAACTGGCTGGCCGCA‑GGCCCCGCCAGaa                                               >  1:885998/1‑146 (MQ=255)                                                                                                             cATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGCGGTGTTAAATGCCAGATAAGCGGC‑GTTCATCCCCGGACGCAGCGTTTAAAGCAAGCGCGGGGCGGCCCGCAAAATAGATAGCTGGCTGGCAGCC‑GGCCCGGCCCGaa                                               >  1:343952/1‑144 (MQ=255)                                                                                                                ggCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCAtt                                      <  1:56795/150‑1 (MQ=255)                                                                                                                    gACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGTTAAATGCCACATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAACGCGGGGCGTCACGCCAAATGGATAGGTGGCTGGCAGCA‑GGCCCGGCCAGAACGTCGCATTcccc                                  >  2:569136/1‑150 (MQ=255)                                                                                                                     aCAGCGGGATTATTTAGCGGCGG‑TTTTGCGGGGGTAAATGCCACAGAAGTGGC‑GTTTATCCCCAGAACCCGCGTTGAAAACCAGAGCGGGTGGTCACGCAAAAAGGGTAGCCAGCTGGCCACC‑AGCCACGCCCGAACGGCCCATAcccccg                                >  2:506953/1‑149 (MQ=255)                                                                                                                         cGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTAAATTCCAGATAGGGGAC‑GGTTATTCCAGGACGCAGCGTTAAACGCAAACGCGGGTCGTCACGCAAAAAGGATAGCTCGCTGGCAGCC‑AGGCAGGGCAGAACGTCGCCTTcccccg                                >  1:744749/1‑145 (MQ=255)                                                                                                                               tattTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGtt                       <  2:619097/150‑1 (MQ=255)                                                                                                                                attTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGttt                      <  1:979286/150‑1 (MQ=255)                                                                                                                                attTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGttt                      <  2:378758/150‑1 (MQ=255)                                                                                                                                     gcggcgG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTc                     >  2:125731/1‑146 (MQ=255)                                                                                                                                          gg‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTc           >  2:46970/1‑151 (MQ=255)                                                                                                                                               tGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCgg        >  1:1041271/1‑150 (MQ=255)                                                                                                                                                 gCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCg         <  1:523674/147‑1 (MQ=255)                                                                                                                                                   ggTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgcc    >  1:215186/1‑150 (MQ=255)                                                                                                                                                    gtgtTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccgc  >  2:227151/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                        CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGC  >  NC_000913/1356054‑1356350

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