Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,441,538	C→T	*213* (TAG→TAA)	msrA ←	methionine sulfoxide reductase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,441,538	0	C	T	100.0%	290.0 / NA	89	*213* (TAG→TAA)	msrA	methionine sulfoxide reductase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (48/41);  total (48/41)

CTATGCCGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTG  >  NC_000913/4441403‑4441667                                                                                                                                        |                                                                                                                                  cTATGCCGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTAt                                                                                                                                <  2:375167/139‑1 (MQ=255)   aTGCCGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGc                                                                                                                              >  1:1602532/1‑139 (MQ=255)         aTCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCcgg                                                                                                                        >  2:41109/1‑139 (MQ=255)                cTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcaga                                                                                                                 <  2:312378/139‑1 (MQ=255)                cTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcaga                                                                                                                 >  1:312378/1‑139 (MQ=255)                                atatatTGAATTTGGACGGTTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                                 <  1:288571/139‑1 (MQ=255)                                  atatTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgcca                                                                                               >  1:816832/1‑139 (MQ=255)                                  atatTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgcca                                                                                               >  1:1753739/1‑139 (MQ=255)                                      tGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                           <  2:816832/139‑1 (MQ=255)                                      tGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                           >  1:508647/1‑139 (MQ=255)                                             ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                                 <  2:634665/126‑1 (MQ=255)                                             ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                                 >  1:634665/1‑126 (MQ=255)                                             ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTa                                                                                    >  2:1390733/1‑139 (MQ=255)                                             ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTa                                                                                    >  1:514510/1‑139 (MQ=255)                                                  aTTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCAt                                                                               <  2:1602532/139‑1 (MQ=255)                                                     ttGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACg                                                                            >  2:633690/1‑139 (MQ=255)                                                       gTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACggg                                                                          <  2:508647/139‑1 (MQ=255)                                                          ggCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGttt                                                                       <  2:514510/139‑1 (MQ=255)                                                             cTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                           >  2:331329/1‑116 (MQ=255)                                                             cTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                           <  1:331329/116‑1 (MQ=255)                                                             cTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTa                                                                    >  2:502810/1‑139 (MQ=255)                                                             cTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTa                                                                    >  1:1228758/1‑139 (MQ=255)                                                                gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCATTTCCGCCAATTCCACAGTAACCATACGGGTTTTTATg                                                                  >  2:1001239/1‑138 (MQ=255)                                                                gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATg                                                                  <  1:1001239/138‑1 (MQ=255)                                                                gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGc                                                                 >  1:1933882/1‑139 (MQ=255)                                                                gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGc                                                                 >  1:428198/1‑139 (MQ=255)                                                                gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGc                                                                 >  2:601470/1‑139 (MQ=255)                                                                     tCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGata                                                            >  1:1531882/1‑139 (MQ=255)                                                                     tCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGACGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGata                                                            >  1:305962/1‑139 (MQ=255)                                                                        tttACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATg                                                                  >  2:1026001/1‑130 (MQ=255)                                                                        tttACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATg                                                                  <  1:1026001/130‑1 (MQ=255)                                                                          tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCcgg                                                                                                                        >  1:159745/1‑74 (MQ=255)                                                                          tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCcgg                                                                                                                        <  2:159745/74‑1 (MQ=255)                                                                          tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCt                                                       >  2:1205523/1‑139 (MQ=255)                                                                            cgccgcATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCtgg                                                     <  2:1531882/139‑1 (MQ=255)                                                                               cgcATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcaga                                                                                                                 >  1:1357761/1‑76 (MQ=255)                                                                               cgcATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcaga                                                                                                                 <  2:1357761/76‑1 (MQ=255)                                                                               cgcATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCtggtgg                                                  <  1:1205523/139‑1 (MQ=255)                                                                                    cAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATc                                             <  1:502810/139‑1 (MQ=255)                                                                                      ggCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCtt                                           <  2:1933882/139‑1 (MQ=255)                                                                                           aGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCtt                                           >  2:1917334/1‑134 (MQ=255)                                                                                           aGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCtt                                           <  1:1917334/134‑1 (MQ=255)                                                                                           aGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCa                                      <  2:1753739/139‑1 (MQ=255)                                                                                             aCAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                        <  1:728212/135‑1 (MQ=255)                                                                                             aCAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                        >  2:728212/1‑135 (MQ=255)                                                                                                aaGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCtggtgg                                                  >  1:340264/1‑122 (MQ=255)                                                                                                aaGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCtggtgg                                                  <  2:340264/122‑1 (MQ=255)                                                                                                 agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCcg                                                                                                                         >  2:1835271/1‑50 (MQ=255)                                                                                                 agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCcg                                                                                                                         <  1:1835271/50‑1 (MQ=255)                                                                                                 agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATa                                                                              >  2:486598/1‑93 (MQ=255)                                                                                                 agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATa                                                                              <  1:486598/93‑1 (MQ=255)                                                                                                 agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                        <  1:582306/131‑1 (MQ=255)                                                                                                 agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                        >  2:582306/1‑131 (MQ=255)                                                                                                 agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGATTCCGGTGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCAtaataa                                >  1:880189/1‑139 (MQ=255)                                                                                                  gagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                                 >  2:1938615/1‑73 (MQ=255)                                                                                                  gagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                                 <  1:1938615/73‑1 (MQ=255)                                                                                                  gagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATaaa                               <  2:428198/139‑1 (MQ=255)                                                                                                     cACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTc                                                                                                                           <  2:296200/44‑1 (MQ=255)                                                                                                     cACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTc                                                                                                                           >  1:296200/1‑44 (MQ=255)                                                                                                     cACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                           <  2:1506276/76‑1 (MQ=255)                                                                                                     cACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaattc                                                                                           >  1:1506276/1‑76 (MQ=255)                                                                                                      aCTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACgg                           >  2:1915480/1‑139 (MQ=255)                                                                                                       cTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGgt                          >  1:1256653/1‑139 (MQ=255)                                                                                                        tttATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCGATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCa                                                                <  2:1253028/100‑1 (MQ=255)                                                                                                        tttATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCa                                                                >  1:1253028/1‑100 (MQ=255)                                                                                                        tttATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGgtg                         >  2:461817/1‑139 (MQ=255)                                                                                                             caacaaACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTC‑‑CAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGtt                  <  1:41109/139‑1 (MQ=255)                                                                                                                  aaCTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGc               >  2:389790/1‑139 (MQ=255)                                                                                                                  aaCTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGc               >  1:118241/1‑139 (MQ=255)                                                                                                                   aCTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATa                                                                              >  2:1199675/1‑75 (MQ=255)                                                                                                                   aCTACATATGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATa                                                                              <  1:1199675/75‑1 (MQ=255)                                                                                                                     tACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCTTAATAAAACGGTGTGGCGTTAGCGAt            <  1:1915480/139‑1 (MQ=255)                                                                                                                          ttGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCg       <  1:461817/139‑1 (MQ=255)                                                                                                                             tACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaatt                                                                                            <  1:1033001/51‑1 (MQ=255)                                                                                                                             tACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaatt                                                                                            >  2:1033001/1‑51 (MQ=255)                                                                                                                              aCCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtggt   >  1:738636/1‑139 (MQ=255)                                                                                                                              aCCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGCTGTGGCGTTAGCGATTTCCgtggt   <  2:1256653/139‑1 (MQ=255)                                                                                                                               cccGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtggtg  <  1:633690/139‑1 (MQ=255)                                                                                                                               cccGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtggtg  >  1:756507/1‑139 (MQ=255)                                                                                                                                 cGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                                 >  2:1848891/1‑42 (MQ=255)                                                                                                                                 cGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                                 <  1:1848891/42‑1 (MQ=255)                                                                                                                                 cGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCAt                                                                               >  2:1508496/1‑60 (MQ=255)                                                                                                                                 cGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCAt                                                                               <  1:1508496/60‑1 (MQ=255)                                                                                                                                  gTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGt                                                                                     <  2:1702628/53‑1 (MQ=255)                                                                                                                                  gTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGt                                                                                     >  1:1702628/1‑53 (MQ=255)                                                                                                                                   tAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                        >  2:1523443/1‑97 (MQ=255)                                                                                                                                   tAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                        <  1:1523443/97‑1 (MQ=255)                                                                                                                                    aaCGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaat                                                                                             >  1:1892439/1‑43 (MQ=255)                                                                                                                                    aaCGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaat                                                                                             <  2:1892439/43‑1 (MQ=255)                                                                                                                                        |                                                                                                                                  CTATGCCGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTG  >  NC_000913/4441403‑4441667

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