Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	195,331	G→A	6.1%	A143A (GCG→GCA)	ispU →	undecaprenyl pyrophosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	195,331	0	G	A	6.1%	212.7 / 3.8	65	A143A (GCG→GCA)	ispU	undecaprenyl pyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (33/28);  new base A (2/2);  total (35/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 = 1.27e-01

AGCGAAGTAAAAAGTCTGCACCGACATAACGTGCGTCTGCGTATTATTGGCGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGCTAAACCAGCATGTCTGTATGCATGAACTGGCCCCTGTAGATTTAG  >  NC_000913/195185‑195476                                                                                                                                                   |                                                                                                                                                  aGCGAAGTAAAAAGTCTGCACCGACATAACGTGCGTCTGCGTATTATTGGCGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGaa                                                                                                                                                 >  2:36340/1‑149 (MQ=255)        tAAAAAGTCTGCACCGACATAACGTGCGTCTGCGTATTATTGGCGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCATAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGt                                                                                                                                          >  2:324874/1‑149 (MQ=255)               tCTGCACCGACATAACGTGCGTCTGCGTATTATTGGCGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGtt                                                                                                                                   <  4:378317/149‑1 (MQ=255)                   cACCGACATAACGTGCGTCTGCGTATTATTGGCGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGa                                                                                                                               <  2:22591/149‑1 (MQ=255)                              cGTGCGTCTGCGTATTATTGGCGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAgg                                                                                                                    <  4:454194/149‑1 (MQ=255)                               gTGCGTCTGCGTATTATTGGCGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAggg                                                                                                                   >  2:19060/1‑149 (MQ=255)                                         gTATTATTGGCGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCaa                                                                                                         <  1:19060/149‑1 (MQ=255)                                          tattatTGGCGATACCAGTCGCTTTAACTCGCGTTTGCAAGGACGTATTCGTAAAGCTGGAGCGATAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAAATACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAAc                                                                                                        >  4:352050/1‑149 (MQ=255)                                          tattatTGGCGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAAc                                                                                                        <  4:262867/149‑1 (MQ=255)                                               ttGGCGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCt                                                                                                   >  2:323848/1‑149 (MQ=255)                                                  gCGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCt                                                                                                   >  3:50633/1‑146 (MQ=255)                                                  gCGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCt                                                                                                   <  4:50633/146‑1 (MQ=255)                                                  gCGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGaa                                                                                                >  3:447446/1‑149 (MQ=255)                                                    gATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGaaaa                                                                                              <  3:166006/149‑1 (MQ=255)                                                     aTACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAg                                                                                             <  1:145345/149‑1 (MQ=255)                                                         cAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTgc                                                                                         >  4:394757/1‑149 (MQ=255)                                                           gTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTgcag                                                                                       <  1:382125/149‑1 (MQ=255)                                                               cTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAg                                                                                                                          <  3:7932/110‑1 (MQ=255)                                                               cTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAg                                                                                                                          >  4:7932/1‑110 (MQ=255)                                                               cTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAg                                                                                   >  1:209733/1‑149 (MQ=255)                                                               cTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAg                                                                                   >  4:18545/1‑149 (MQ=255)                                                               cTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAg                                                                                   >  4:267011/1‑149 (MQ=255)                                                                tttAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGt                                                                                                                                          <  3:285546/93‑1 (MQ=255)                                                                tttAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGt                                                                                                                                          >  4:285546/1‑93 (MQ=255)                                                                tttAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAgg                                                                                  >  3:93822/1‑149 (MQ=255)                                                                tttAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAgg                                                                                  >  1:401643/1‑149 (MQ=255)                                                                tttAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAgg                                                                                  >  1:314837/1‑149 (MQ=255)                                                                tttAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAgg                                                                                  >  1:219071/1‑149 (MQ=255)                                                                 ttAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAg                                                                                             >  4:398383/1‑137 (MQ=255)                                                                 ttAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAg                                                                                             <  3:398383/137‑1 (MQ=255)                                                                    aCTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAAc                                                                              >  3:472985/1‑149 (MQ=255)                                                                      tCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCt                                                                            <  3:267011/149‑1 (MQ=255)                                                                          gTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCaa                                                                        <  2:219071/149‑1 (MQ=255)                                                                           tttGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAAc                                                                       <  3:394757/149‑1 (MQ=255)                                                                           tttGCAAGAACATATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAAc                                                                       >  1:250554/1‑149 (MQ=255)                                                                              gCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACcag                                                                    >  2:281532/1‑149 (MQ=255)                                                                              gCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACcag                                                                    >  4:423987/1‑149 (MQ=255)                                                                                aaGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCAAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACcagat                                                                  >  3:160458/1‑149 (MQ=255)                                                                                aaGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCAAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACcagat                                                                  >  3:158549/1‑149 (MQ=255)                                                                                  gAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACcagatca                                                                <  3:194968/149‑1 (MQ=255)                                                                                   aaCGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACcagatcag                                                               >  2:462539/1‑149 (MQ=255)                                                                                        aTTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCagataga                                                          <  3:16181/149‑1 (MQ=255)                                                                                         ttCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCagatagat                                                         <  4:472985/149‑1 (MQ=255)                                                                                         ttCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCagatagat                                                         <  1:459771/149‑1 (MQ=255)                                                                                            gTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGaa                                                      <  4:447446/149‑1 (MQ=255)                                                                                            gTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGaa                                                      >  4:94859/1‑149 (MQ=255)                                                                                                 tCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGGAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGAt                                                 >  1:249655/1‑149 (MQ=255)                                                                                                  cTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATg                                                <  1:36340/149‑1 (MQ=255)                                                                                                   tGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGc                                               >  3:244765/1‑149 (MQ=255)                                                                                                   tGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGc                                               <  2:209733/149‑1 (MQ=255)                                                                                                     aaGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGCTa                                             >  2:112602/1‑149 (MQ=255)                                                                                                            aaCAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGCTAAACCAGc                                      >  2:87023/1‑149 (MQ=255)                                                                                                                gCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGCTAAACCAGCATGt                                  <  1:112602/149‑1 (MQ=255)                                                                                                                    ggAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGCTAAACCAGCATGTCTGt                              >  4:35466/1‑149 (MQ=255)                                                                                                                     gAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGCTAAACCAGCATGTCTGTa                             <  3:18545/149‑1 (MQ=255)                                                                                                                        tACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGCTAAACCAGCATGTCTGTATGc                          <  1:323848/149‑1 (MQ=255)                                                                                                                             gTCTGACGCTGAATATTGCGGCAAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGCTAAACCAGCATGTCTGTATGCATGaa                     <  4:160458/149‑1 (MQ=255)                                                                                                                             gTCTGACGCTGAATATTGCGGCAAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGCTAAACCAGCATGTCTGTATGCATGaa                     <  4:158549/149‑1 (MQ=255)                                                                                                                                tGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGCTAAACCAGCATGTCTGTATGCATGAACTg                  <  3:444512/149‑1 (MQ=255)                                                                                                                                    gCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGCTAAAc                                          >  4:261510/1‑121 (MQ=255)                                                                                                                                    gCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGCTAAAc                                          <  3:261510/121‑1 (MQ=255)                                                                                                                                    gCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGCTAAACCAGCATGTCTGTATGCATGAACTGGccc              >  2:163714/1‑149 (MQ=255)                                                                                                                                      tGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAg                                                                                   >  2:39766/1‑78 (MQ=255)                                                                                                                                      tGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAg                                                                                   <  1:39766/78‑1 (MQ=255)                                                                                                                                            ttGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGCTAAACCAGCATGTCTGTATGCATGAACTGGCCCCTGTAGAt      >  2:237893/1‑149 (MQ=255)                                                                                                                                               cggcgAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGCTAAACCAGCATGTCTGTATGCATGAACTGGCCCCTGTAGATTTa   <  1:87023/149‑1 (MQ=255)                                                                                                                                                ggcgAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGCTAAACCAGCATGTCTGTATGCATGAACTGGCCCCTGTAGATTTAg  <  4:299668/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  AGCGAAGTAAAAAGTCTGCACCGACATAACGTGCGTCTGCGTATTATTGGCGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAACTGGCTGAAAAGGTGCAGCAAGGAAACCTGCAACCAGATCAGATAGATGAAGAGATGCTAAACCAGCATGTCTGTATGCATGAACTGGCCCCTGTAGATTTAG  >  NC_000913/195185‑195476

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