Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,786,729	G→A	pseudogene (333/2253 nt)	ygaQ →	uncharacterized protein, putative enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,786,729	0	G	A	100.0%	314.9 / NA	95	pseudogene (333/2253 nt)	ygaQ	uncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (47/48);  total (47/48)

AAAGTTATTAATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTA  >  NC_000913/2786598‑2786862                                                                                                                                    |                                                                                                                                      aaaGTTATTAATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATc                                                                                                                                <  1:579021/139‑1 (MQ=255)         taATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGAt                                                                                                                        >  2:823558/1‑139 (MQ=255)              aTCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGAt                                                                                                                        <  1:324462/134‑1 (MQ=255)              aTCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGAt                                                                                                                        >  2:324462/1‑134 (MQ=255)               tCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGa                                                                                                                 >  1:1020234/1‑139 (MQ=255)               tCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGa                                                                                                                 >  1:300377/1‑139 (MQ=255)                     aaaaaaCGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                           <  1:1155277/139‑1 (MQ=255)                     aaaaaaCGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                           <  2:67733/139‑1 (MQ=255)                     aaaaaaCGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                           <  2:511554/139‑1 (MQ=255)                             aTTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCaa                                                                                                                              <  1:193171/113‑1 (MQ=255)                             aTTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCaa                                                                                                                              >  2:193171/1‑113 (MQ=255)                                      tttCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGgt                                                                                          <  2:300377/139‑1 (MQ=255)                                       ttCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGgtg                                                                                         <  2:1166287/139‑1 (MQ=255)                                          ccTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                      <  2:1020234/139‑1 (MQ=255)                                            tttCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTAt                                                                                    <  2:890049/139‑1 (MQ=255)                                              tctcCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAaga                                                                                                         >  2:1035926/1‑116 (MQ=255)                                              tctcCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAaga                                                                                                         <  1:1035926/116‑1 (MQ=255)                                               ctcCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                      >  2:244271/1‑134 (MQ=255)                                               ctcCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                      <  1:244271/134‑1 (MQ=255)                                                 ccAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTa                                                                               <  2:548962/139‑1 (MQ=255)                                                        cAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAagag                                                                                                        >  2:1143379/1‑107 (MQ=255)                                                        cAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAagag                                                                                                        <  1:1143379/107‑1 (MQ=255)                                                         aaCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATc                                                                                                                                <  1:1022712/83‑1 (MQ=255)                                                         aaCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATc                                                                                                                                >  2:1022712/1‑83 (MQ=255)                                                          aCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGc                                                                      >  2:59811/1‑139 (MQ=255)                                                            ttACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTg                                                                    >  2:766826/1‑139 (MQ=255)                                                              aCCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGaa                                                                  >  2:697872/1‑139 (MQ=255)                                                               ccGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCg                                                                                            >  1:386200/1‑112 (MQ=255)                                                               ccGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCg                                                                                            <  2:386200/112‑1 (MQ=255)                                                                   gagTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACAttt                                                             >  2:1045403/1‑139 (MQ=255)                                                                       tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                         <  2:105965/139‑1 (MQ=255)                                                                       tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                         <  1:1045403/139‑1 (MQ=255)                                                                                   cacaGTGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATgg                                                                            >  2:671154/1‑108 (MQ=255)                                                                                   cacaGTGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATgg                                                                            <  1:671154/108‑1 (MQ=255)                                                                                   cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCg                                                                                            <  2:1350573/92‑1 (MQ=255)                                                                                   cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCg                                                                                            >  1:1350573/1‑92 (MQ=255)                                                                                   cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAc                                                          >  2:558562/1‑126 (MQ=255)                                                                                   cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAc                                                          <  1:558562/126‑1 (MQ=255)                                                                                        cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAccc                                                        >  2:729700/1‑123 (MQ=255)                                                                                        cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAccc                                                        <  1:729700/123‑1 (MQ=255)                                                                                        cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCcaca                                                     <  1:662392/126‑1 (MQ=255)                                                                                        cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCcaca                                                     <  2:1133433/126‑1 (MQ=255)                                                                                        cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCcaca                                                     >  1:1133433/1‑126 (MQ=255)                                                                                        cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCcaca                                                     >  2:662392/1‑126 (MQ=255)                                                                                        cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACaa                                        >  2:3775/1‑139 (MQ=255)                                                                                        cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACaa                                        >  2:1073910/1‑139 (MQ=255)                                                                                        cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACaa                                        >  2:830869/1‑139 (MQ=255)                                                                                        cGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACaa                                        >  1:77913/1‑139 (MQ=255)                                                                                                 aCAGCGCCAGACAGGATAAGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTcccc                               <  1:129548/139‑1 (MQ=255)                                                                                                  cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                      >  2:353732/1‑83 (MQ=255)                                                                                                  cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                      <  1:353732/83‑1 (MQ=255)                                                                                                  cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                         >  1:1313825/1‑112 (MQ=255)                                                                                                  cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                         <  2:1313825/112‑1 (MQ=255)                                                                                                  cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCa                              >  1:1343881/1‑139 (MQ=255)                                                                                                  cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCa                              <  1:766826/139‑1 (MQ=255)                                                                                                  cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCa                              <  1:59811/139‑1 (MQ=255)                                                                                                  cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCa                              >  2:328308/1‑139 (MQ=255)                                                                                                    gcgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                           >  1:512850/1‑60 (MQ=255)                                                                                                    gcgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                           <  2:512850/60‑1 (MQ=255)                                                                                                    gcgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGgtg                                                                                         <  1:346419/78‑1 (MQ=255)                                                                                                    gcgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGgtg                                                                                         >  2:346419/1‑78 (MQ=255)                                                                                                     cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCg                                                                                            <  2:273889/74‑1 (MQ=255)                                                                                                     cgcCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCg                                                                                            >  1:273889/1‑74 (MQ=255)                                                                                                       ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATc                                                                                                                                <  1:350049/37‑1 (MQ=255)                                                                                                       ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATc                                                                                                                                >  2:350049/1‑37 (MQ=255)                                                                                                       ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                         <  2:159357/107‑1 (MQ=255)                                                                                                       ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                         >  1:159357/1‑107 (MQ=255)                                                                                                       ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACaaa                                       >  1:1021750/1‑125 (MQ=255)                                                                                                       ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACaaa                                       <  2:1021750/125‑1 (MQ=255)                                                                                                       ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCa                              <  1:549068/134‑1 (MQ=255)                                                                                                       ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCa                              >  2:549068/1‑134 (MQ=255)                                                                                                       ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGa                         >  2:1256602/1‑139 (MQ=255)                                                                                                       ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGa                         >  1:507309/1‑139 (MQ=255)                                                                                                          gACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATaa                      <  1:3775/139‑1 (MQ=255)                                                                                                              ggATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAatat                  <  1:1256602/139‑1 (MQ=255)                                                                                                                   ggAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGc             <  2:1343881/139‑1 (MQ=255)                                                                                                                    gagaAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCa            >  2:534697/1‑139 (MQ=255)                                                                                                                     agaaagaaAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCat           <  2:77913/139‑1 (MQ=255)                                                                                                                       aaagaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCatat         >  2:576572/1‑139 (MQ=255)                                                                                                                         agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                      <  2:872838/60‑1 (MQ=255)                                                                                                                         agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGttt                                                                                      >  1:872838/1‑60 (MQ=255)                                                                                                                         agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaa                             <  1:548080/117‑1 (MQ=255)                                                                                                                         agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaa                             >  2:548080/1‑117 (MQ=255)                                                                                                                         agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGa                         >  1:429950/1‑121 (MQ=255)                                                                                                                         agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGa                         <  2:429950/121‑1 (MQ=255)                                                                                                                         agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc       <  1:414771/139‑1 (MQ=255)                                                                                                                         agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc       <  1:1073910/139‑1 (MQ=255)                                                                                                                         agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc       >  2:1252946/1‑139 (MQ=255)                                                                                                                         agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc       >  2:414771/1‑139 (MQ=255)                                                                                                                          gaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                         <  2:705327/88‑1 (MQ=255)                                                                                                                          gaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                         >  1:705327/1‑88 (MQ=255)                                                                                                                            aaaaCATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTAt    <  1:534697/139‑1 (MQ=255)                                                                                                                              aaCATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGa                                                                                                    <  1:810667/41‑1 (MQ=37)                                                                                                                              aaCATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGa                                                                                                    >  2:810667/1‑41 (MQ=37)                                                                                                                              aaCATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTa  >  2:382685/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                      AAAGTTATTAATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCTATTA  >  NC_000913/2786598‑2786862

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