Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,830,940	A→G	T162A (ACG→GCG)	yicH →	AsmA famly protein YicH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,830,940	0	A	G	97.1%	230.8 / ‑2.7	69	T162A (ACG→GCG)	yicH	AsmA famly protein YicH
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (32/35);  total (33/36)
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.32e-01

GCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCC  >  NC_000913/3830804‑3831068                                                                                                                                         |                                                                                                                                 gcTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGa                                                                                                                               <  2:837637/140‑1 (MQ=255)   tACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAg                                                                                                                             >  1:575346/1‑140 (MQ=255)      cGTTCAAAGCCGATCGTCTG‑AACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGgcgc                                                                                                                         <  1:534115/140‑1 (MQ=255)       gTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGgcgc                                                                                                                         <  1:1043412/140‑1 (MQ=255)               ccGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAg                                                                                                                 >  2:95337/1‑140 (MQ=255)                 gATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGtt                                                                                                               <  1:886479/140‑1 (MQ=255)                     gTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGt                                                                                                           <  1:203001/140‑1 (MQ=255)                               gCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGAt                                                                                                    <  2:1303900/137‑1 (MQ=255)                               gCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGAt                                                                                                    >  1:1303900/1‑137 (MQ=255)                                cGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCt                                                                                                >  1:233434/1‑140 (MQ=255)                                cGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCt                                                                                                >  2:638460/1‑140 (MQ=255)                                 gTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCtt                                                                                               >  1:124107/1‑140 (MQ=255)                                 gTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCtt                                                                                               <  1:1269600/140‑1 (MQ=255)                                          cGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCt                                                                                          >  1:352320/1‑136 (MQ=255)                                          cGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCt                                                                                          <  2:352320/136‑1 (MQ=255)                                          cGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAat                                                                                      >  1:314408/1‑140 (MQ=255)                                          cGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAat                                                                                      >  1:508536/1‑140 (MQ=255)                                           gTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAatg                                                                                     >  2:1012727/1‑140 (MQ=255)                                           gTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAatg                                                                                     <  1:1126434/140‑1 (MQ=255)                                             ttAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAatgat                                                                                   <  2:1342495/140‑1 (MQ=255)                                               aaTAGCCCGAATAGCGACTGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTTGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGt                                                                                 >  2:567241/1‑140 (MQ=255)                                               aaTAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGt                                                                                 >  1:1192300/1‑140 (MQ=255)                                                  aGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTcc                                                                              <  2:233434/140‑1 (MQ=255)                                                    cccGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTg                                                                            >  2:884374/1‑140 (MQ=255)                                                     ccGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGc                                                                           <  2:1374779/140‑1 (MQ=255)                                                     ccGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGc                                                                           <  2:1358072/140‑1 (MQ=255)                                                     ccGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGc                                                                           <  2:634260/140‑1 (MQ=255)                                                     ccGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGc                                                                           <  1:357604/140‑1 (MQ=255)                                                      cGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGcc                                                                          >  1:261150/1‑140 (MQ=255)                                                      cGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGcc                                                                          >  2:137205/1‑140 (MQ=255)                                                       gAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGcca                                                                         <  1:638460/140‑1 (MQ=255)                                                       gAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGcca                                                                         <  2:124107/140‑1 (MQ=255)                                                        aaTAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGccac                                                                        >  1:886972/1‑140 (MQ=255)                                                        aaTAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGccac                                                                        >  2:100943/1‑140 (MQ=255)                                                           aGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCaa                                                                     <  2:850316/140‑1 (MQ=255)                                                             cGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATg                                                                   >  1:873230/1‑140 (MQ=255)                                                              gAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGt                                                                  >  2:878761/1‑140 (MQ=255)                                                                 tGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACt                                                               >  2:1240036/1‑140 (MQ=255)                                                                  ggAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTg                                                              >  1:681835/1‑140 (MQ=255)                                                                  ggAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTg                                                              >  1:754109/1‑140 (MQ=255)                                                                  ggAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTg                                                              >  1:756592/1‑140 (MQ=255)                                                                     aaCTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGAtt                                                           <  2:873230/140‑1 (MQ=255)                                                                     aaCTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGAtt                                                           <  2:681835/140‑1 (MQ=255)                                                                       cTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGa                                                         >  1:538905/1‑140 (MQ=255)                                                                         gAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAg                                                       <  2:754109/140‑1 (MQ=255)                                                                          aGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGACGGCGCAGATTCAGTTTAGTGCCGGATCGCTTCCGCTCAATGCTGTTCCTGCCCCCCATGTACTGATTGAAgg                                                      <  2:575346/140‑1 (MQ=255)                                                                          aGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAgg                                                      >  1:340704/1‑140 (MQ=255)                                                                          aGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAgg                                                      <  2:1192400/140‑1 (MQ=255)                                                                          aGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAgg                                                      <  2:1192300/140‑1 (MQ=255)                                                                          aGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAgg                                                      <  1:872811/140‑1 (MQ=255)                                                                          aGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAgg                                                      <  1:137205/140‑1 (MQ=255)                                                                           gcgcgcAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGc                                                     <  1:36459/140‑1 (MQ=255)                                                                            cgcgcAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCa                                                    >  1:660462/1‑140 (MQ=255)                                                                            cgcgcAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCa                                                    <  1:619321/140‑1 (MQ=255)                                                                            cgcgcAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCa                                                    >  2:1059741/1‑140 (MQ=255)                                                                             gcgcAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAg                                                   <  2:314408/140‑1 (MQ=255)                                                                               gcAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTa                                                 <  1:884374/140‑1 (MQ=255)                                                                                         aaTGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGa                                       <  2:756592/140‑1 (MQ=255)                                                                                               ggcgTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGt                                 >  1:493541/1‑140 (MQ=255)                                                                                                gcgTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGtt                                >  2:268017/1‑140 (MQ=255)                                                                                                  gTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTAc                              <  1:1240036/140‑1 (MQ=255)                                                                                                   tAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACg                             <  1:95337/140‑1 (MQ=255)                                                                                                      ttCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTg                          <  1:100943/140‑1 (MQ=255)                                                                                                       tCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGa                         <  1:878761/140‑1 (MQ=255)                                                                                                           tGGTCACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTaa                     >  1:140382/1‑140 (MQ=255)                                                                                                               cACCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTg                 <  2:140382/140‑1 (MQ=255)                                                                                                                aCCAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTgg                >  2:1222293/1‑140 (MQ=255)                                                                                                                  cAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGt              >  2:442881/1‑140 (MQ=255)                                                                                                                  cAAAAGCCGGTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGt              >  2:1352710/1‑140 (MQ=255)                                                                                                                          ggTAAAGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACAt      >  1:1327062/1‑140 (MQ=255)                                                                                                                              aaGTGCTGGGTGCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGcc  <  1:1059741/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                 GCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAGTGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTGACTAACCTGGGTGCCGACATCGCC  >  NC_000913/3830804‑3831068

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