Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,128,229	A→T	7.5%	A1200A (GCA→GCT)	PP_1887 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,128,229	0	A	T	7.5%	180.3 / 3.0	62	A1200A (GCA→GCT)	PP_1887	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base T (3/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GAAGCGCCCGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGTGAT  >  NC_002947/2128083‑2128376                                                                                                                                                     |                                                                                                                                                     gAAGCGCCCGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGcac                                                                                                                                                     >  3:79183/1‑148 (MQ=255)       cccGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAg                                                                                                                                             >  4:142308/1‑149 (MQ=255)        ccGGGGACCCCTCCGGTCAAACCATTGCGAGGCAGATATTCGGCTTCGATGGCCTGGACAACGTCCTCTCTGTGATCACCTACAGACCGGTCACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑Ga                                                                                                                                              <  1:292918/149‑1 (MQ=255)              aCCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCTCATTT‑GTGCCGGATc                                                                                                                                      <  2:367060/149‑1 (MQ=255)                cccTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCAt                                                                                                                                    >  8:105887/1‑149 (MQ=255)                             ccATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATgcg                                                                                                                       >  3:132070/1‑149 (MQ=255)                                   ccAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTg                                                                                                                 <  1:336882/149‑1 (MQ=255)                                    cAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATCGAGCCGGATCATTCCGCCTGATGCGCTACt                                                                                                                  <  1:226599/148‑1 (MQ=255)                                      ggCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGcc                                                                                                                              >  6:156328/1‑133 (MQ=255)                                      ggCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGcc                                                                                                                              <  5:156328/133‑1 (MQ=255)                                       gCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGgaga                                                                                                             >  3:179564/1‑149 (MQ=255)                                            tatTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCTCCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTAt                                                                                                        <  5:320999/149‑1 (MQ=255)                                            tatTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTAt                                                                                                        <  3:142308/149‑1 (MQ=255)                                              ttCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGa                                                                                                      >  8:65764/1‑149 (MQ=255)                                                cGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACg                                                                                                    >  3:4206/1‑149 (MQ=255)                                                cGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACg                                                                                                    >  3:9487/1‑149 (MQ=255)                                                 ggCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGt                                                                                                   <  7:105887/149‑1 (MQ=255)                                                 ggCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGt                                                                                                   >  5:211287/1‑149 (MQ=255)                                                    ttCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGa                                                                                                >  3:92813/1‑149 (MQ=255)                                                    ttCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGa                                                                                                >  4:68583/1‑149 (MQ=255)                                                            ccTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGga                                                                                                               <  1:254926/126‑1 (MQ=255)                                                            ccTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGga                                                                                                               >  2:254926/1‑126 (MQ=255)                                                              tGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTg                                                                                      >  6:281167/1‑149 (MQ=255)                                                                gACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTg                                                                                      <  3:68583/147‑1 (MQ=255)                                                                          tctcGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑TAGCCGGATCATTCCGCCTGATCCGATAATGGAGACCTATGAAGTCGCTCACGACCTGTAGGAgttgtt                                                                           >  7:352362/1‑148 (MQ=255)                                                                          tctcGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAgttgtt                                                                           >  1:353787/1‑148 (MQ=255)                                                                          tctcGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAgttgtt                                                                           >  6:147296/1‑148 (MQ=255)                                                                              ggTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAg                                                                       <  4:92813/148‑1 (MQ=255)                                                                              ggTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGa                                                                      <  3:324464/149‑1 (MQ=255)                                                                                tGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTc                                                                          >  5:211653/1‑143 (MQ=255)                                                                                   tCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCTGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTc                                                                 <  4:79183/149‑1 (MQ=255)                                                                                    cACCTACAGACCGG‑‑ACAGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTc                                                                 >  7:40299/1‑148 (MQ=255)                                                                                     aCCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCa                                                                >  7:190349/1‑148 (MQ=255)                                                                                     aCCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCa                                                                >  5:56284/1‑148 (MQ=255)                                                                                       cTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAAt                                                             <  8:305083/149‑1 (MQ=255)                                                                                          cAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGAGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATAc                                                           >  8:279200/1‑148 (MQ=255)                                                                                                    aCGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAg                                                  <  7:65764/149‑1 (MQ=255)                                                                                                           gTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTa                                           <  6:211287/149‑1 (MQ=255)                                                                                                            tGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCa                                                                         >  6:294535/1‑118 (MQ=255)                                                                                                            tGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCa                                                                         <  5:294535/118‑1 (MQ=255)                                                                                                            tGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTAt                                          <  4:179564/149‑1 (MQ=255)                                                                                                             ggCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATc                                         >  1:327636/1‑149 (MQ=255)                                                                                                                 gCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATTACGTCGATCACGACCTGGAGTAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGAtt                                     <  6:56284/149‑1 (MQ=255)                                                                                                                  cGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGAttt                                    <  4:132070/149‑1 (MQ=255)                                                                                                                   ggACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTg                                   >  5:271570/1‑149 (MQ=255)                                                                                                                    gACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTg                                   >  6:143588/1‑148 (MQ=255)                                                                                                                      ccTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCAt                                <  2:327636/149‑1 (MQ=255)                                                                                                                          gTACGAATTCAAGTATAGCGATCCTGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTa                              <  4:9487/147‑1 (MQ=255)                                                                                                                          gTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTa                              <  4:4206/147‑1 (MQ=255)                                                                                                                            aCGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTacgac                          >  1:47835/1‑149 (MQ=255)                                                                                                                             cGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTacgacg                         >  7:80103/1‑149 (MQ=255)                                                                                                                             cGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTacgac                          <  3:221987/148‑1 (MQ=255)                                                                                                                               aaTTCAAGAATAGCGATCCGGCACTAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGCGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGc                       <  2:306387/149‑1 (MQ=255)                                                                                                                                 ttCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCaa                     >  6:192106/1‑149 (MQ=255)                                                                                                                                    aaGAATAGCGATCCGTCTCAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGTAGACCTATGACGTCGATCACGACCTGGAGTAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCTATgg                  <  6:376720/149‑1 (MQ=255)                                                                                                                                    aaGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGAtt                                     >  5:236507/1‑130 (MQ=255)                                                                                                                                    aaGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGAtt                                     <  6:236507/130‑1 (MQ=255)                                                                                                                                    aaGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATg                   <  7:279200/148‑1 (MQ=255)                                                                                                                                        ggaatgATCCGGCTCATTT‑GTGCCGGATCGTTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAAc              >  8:14961/6‑149 (MQ=255)                                                                                                                                             gATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGAt         <  5:147296/149‑1 (MQ=255)                                                                                                                                              aTCCGGCTCATTT‑GTGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGAt         >  7:142458/1‑148 (MQ=255)                                                                                                                                               tCCGGCTCATTT‑GTGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGAt         >  1:327092/1‑147 (MQ=255)                                                                                                                                               tCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCa       <  2:353787/149‑1 (MQ=255)                                                                                                                                                    cacaAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGtgat  <  5:192106/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GAAGCGCCCGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTTCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGG‑‑ACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAT‑GAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGTGAT  >  NC_002947/2128083‑2128376

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