Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,369,520	T→C	V43A (GTT→GCT)	pspD →	peripheral inner membrane phage‑shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,369,520	0	T	C	100.0%	186.3 / NA	56	V43A (GTT→GCT)	pspD	peripheral inner membrane phage‑shock protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/38);  total (18/38)

GTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGTTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAATGGATGGCA  >  NC_000913/1369381‑1369662                                                                                                                                             |                                                                                                                                                    gtgAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCg                                                                                                                                            <  2:511298/149‑1 (MQ=255)  tgAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGGTGGGGGATAAAATCAGCTGCTCGCCGA‑cc                                                                                                                                        >  1:907409/1‑151 (MQ=255)   gAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCg                                                                                                                                       <  2:82752/151‑1 (MQ=255)       aaaGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGGCGGGGGGTAAAAACACCTGCTCGCCCA‑cc                                                                                                                                        >  1:594253/1‑146 (MQ=255)       aaaGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGa                                                                                                                                   <  1:806405/151‑1 (MQ=255)        aaGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGGGGGAAAAACAGCTGCTCGCCCA‑cccc                                                                                                                                      >  1:702262/1‑145 (MQ=255)                   aCTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTgctg                                                                                                                        <  1:239351/150‑1 (MQ=255)                             aaCAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑Ct                                                                                                           <  1:232568/151‑1 (MQ=255)                              aCAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑g                                                                                                         <  2:223797/151‑1 (MQ=255)                               cAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGGGGGGGGGTAAAAAACGCCGCTCGCCCA‑CCGCCGGAAAAGGtggtgggcgggggggtgg                                                                                                           >  2:593309/1‑139 (MQ=255)                                aGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGa                                                                                                       <  2:740673/151‑1 (MQ=255)                                aGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGa                                                                                                       <  1:724431/151‑1 (MQ=255)                                 ggCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGaa                                                                                                      <  2:821900/151‑1 (MQ=255)                                     gggCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGGGTGGC‑GGGGTGGGGGGTAAAATCACCTGGTCGCCCC‑GCCCTGGAAATGGTGGTGGCTGG‑Gggg‑cgtg                                                                                                         >  2:453700/1‑142 (MQ=255)                                      ggCAAAAGGTAAAGCCTTGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTgctggct                                                                                                                     <  1:644842/134‑1 (MQ=255)                                          aaaGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTa                                                                                             <  2:87910/151‑1 (MQ=255)                                              gTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTc                                                                                         <  1:347274/151‑1 (MQ=255)                                                   gCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTTGGGGATAAAATCAGCTGCTCGCCGA‑CCGCCGGAAAAGTTGCTGGGTGT‑GGGG‑GC‑GGGAACGCTGTTTAGTCggg                                                                                      >  1:465812/1‑149 (MQ=255)                                                           tCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCtaata                                                                             <  1:74438/150‑1 (MQ=255)                                                              aaTTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACt                                                                         <  2:745783/151‑1 (MQ=255)                                                               aTTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTg                                                                        <  1:682924/151‑1 (MQ=255)                                                               aTTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTg                                                                        <  2:587414/151‑1 (MQ=255)                                                                 tAGCAGGCAAGCTGGTACTTCTTACCGCGCTGCGCTATGGCCCGGCGGGTGTGGC‑GCGCTGGTCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGTTGT‑GGCT‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTgg                                                                       <  2:614658/150‑1 (MQ=255)                                                                  aGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGca                                                                     <  2:121631/151‑1 (MQ=255)                                                                   gCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGcac                                                                    <  2:749013/151‑1 (MQ=255)                                                                     aGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGGGATAAAAACAGCTGCCCGCCCA‑CCGCCGGAAAAATTTCTGGGTGT‑GGGGCCG‑GGAACGCCGGTTAAG‑CGGGCCGCTAAAAAAATgg                                                                       >  1:164141/1‑146 (MQ=255)                                                                     aGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGcaca                                                                   <  1:736603/150‑1 (MQ=255)                                                                     aGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAg                                                                  <  1:805848/151‑1 (MQ=255)                                                                         aggCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑Gtt                                                             <  1:248767/149‑1 (MQ=255)                                                                         aaGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGGCTGGG‑GGTAAAAACAACTTCGCCCCCACCCCCTGAAAATG‑TGGTGGCCGTGGGGG‑TT‑GGAGCCGC‑GTTTAATCGGGGTGCTAATaaaactggacaaggg                                                               >  1:267841/1‑138 (MQ=255)                                                                          aGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCCTGGGGGTAAAAACAACCGCTCGCCCG‑CCCCCGGAAAAGTTTCTTGGTTG‑GGGG‑CTGGGAACCCTGGTTAAG‑CGGGGTGCCTATAAAACGGCACAGC‑gg                                                              >  2:103463/1‑148 (MQ=255)                                                                              ggTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTtata                                                          <  1:475423/149‑1 (MQ=255)                                                                              ggTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATaa                                                         <  2:371082/150‑1 (MQ=255)                                                                              ggTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATaa                                                         <  1:623432/150‑1 (MQ=255)                                                                               gTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATaaa                                                        <  2:455630/150‑1 (MQ=255)                                                                                        aCCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAAACAACTGCTCGCCGA‑CCCCTGAAAATGTTGCTGGGTGT‑GGCG‑CT‑GGGACCGCTGTTTAGGCGGGCTGCTAATAAACCGGCACAGC‑GGTTTAAAaggggagggg                                              >  1:745047/1‑151 (MQ=255)                                                                                         ccGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGGGGTAAAAACACCCGCTCGCCCC‑CCCCCCGAAATGTTTCTGGGGGG‑GGCG‑GC‑GGGAACGCCGGTAAATCGGGGTGGTAATAAAAAGGCACCGCGGTTTTAAAAgggggggg                                               >  2:208860/1‑150 (MQ=255)                                                                                                gcgcTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCg                                       <  1:438615/150‑1 (MQ=255)                                                                                                gcgcTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCg                                       <  2:164117/150‑1 (MQ=255)                                                                                                        gCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTg                              <  1:323776/151‑1 (MQ=255)                                                                                                            ggcggGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTaa                          <  2:748628/151‑1 (MQ=255)                                                                                                             gcggGTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCCA‑CCGCTGGAAATGTTGCTGGGTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAAATGGCACAGC‑GGTATAAAAGGTGGGGGGGGTTTTGCCAAAATTGttta                          >  2:802047/1‑148 (MQ=255)                                                                                                               gggTGTGGC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGG‑CCGCTGGAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GGTATAAAAGGGGAGGGGGGTTTTGCAAAAATTGTTaaa                         >  2:790112/1‑149 (MQ=255)                                                                                                                     ggC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAgg                  <  1:602560/150‑1 (MQ=255)                                                                                                                     ggC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGc                 <  1:632044/151‑1 (MQ=255)                                                                                                                     ggC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGc                 <  2:797418/151‑1 (MQ=255)                                                                                                                      gC‑GGGCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCg                <  1:819204/151‑1 (MQ=255)                                                                                                                          ggCTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGtat             >  1:48026/1‑151 (MQ=255)                                                                                                                            cTGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGtata            <  2:546916/150‑1 (MQ=255)                                                                                                                             tGGGCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAat          >  2:782337/1‑151 (MQ=255)                                                                                                                                gCGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAatgg        >  2:562758/1‑150 (MQ=255)                                                                                                                                 cGATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAatggat      <  1:434641/151‑1 (MQ=255)                                                                                                                                  gATAAAATCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAatggatg     >  1:598407/1‑151 (MQ=255)                                                                                                                                     aaaaTCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAATGGATGGCa  <  1:65954/151‑1 (MQ=255)                                                                                                                                     aaaaTCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAATGGATGGCa  <  2:339108/151‑1 (MQ=255)                                                                                                                                      aaaTCAGCTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAATGGATGGCa  >  2:394467/1‑150 (MQ=255)                                                                                                                                             |                                                                                                                                                    GTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGC‑GGGCTGGGCGATAAAATCAGTTGCTCGCCGA‑CCGCTGAAAATGTTGCTGGCTGT‑GGCG‑CT‑GGAACCGCTGTTAAGTCGGGCTGCTAATAAACTGGCACAGC‑GTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAATGGATGGCA  >  NC_000913/1369381‑1369662

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