Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,792,826	C→T	*286* (TAG→TAA)	yibB ←	YibB family protein, function unknown

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,792,826	0	C	T	100.0%	293.7 / NA	89	*286* (TAG→TAA)	yibB	YibB family protein, function unknown
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (47/42);  total (47/42)

ACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATA  >  NC_000913/3792700‑3792960                                                                                                                               |                                                                                                                                       aCAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGaa                                                                                                                            >  2:577768/1‑139 (MQ=255)                  tAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAgcgc                                                                                                           <  1:577768/139‑1 (MQ=255)                        tgattgatTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAt                                                                                                     >  2:667297/1‑139 (MQ=255)                          atagcttttACAGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttt                                                                                                   >  1:1091372/6‑139 (MQ=255)                            tgattTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCAttttt                                                                                                 <  2:1695497/139‑1 (MQ=255)                                    cTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTaa                                                                                         <  1:667297/139‑1 (MQ=255)                                     tGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTaaa                                                                                        >  1:275940/1‑139 (MQ=255)                                       cATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTa                                                                                          >  2:1107881/1‑135 (MQ=255)                                       cATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTa                                                                                          <  1:1107881/135‑1 (MQ=255)                                         ttGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGa                                                                                                                             >  2:1761969/1‑98 (MQ=255)                                         ttGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGa                                                                                                                             <  1:1761969/98‑1 (MQ=255)                                         ttGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtgt                                                                                    >  2:1757781/1‑139 (MQ=255)                                                  aaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGa                                                                           <  1:177086/139‑1 (MQ=255)                                                      aaaaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGaaa                                                                                                                           <  1:1407624/87‑1 (MQ=255)                                                      aaaaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGaaa                                                                                                                           >  2:1407624/1‑87 (MQ=255)                                                      aaaaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGaa                                                                                                                       <  2:905486/91‑1 (MQ=255)                                                      aaaaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGaa                                                                                                                       >  1:905486/1‑91 (MQ=255)                                                      aaaaaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAagca                                                                       >  2:326244/1‑139 (MQ=255)                                                             gatgCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAAc                                                                <  2:1640490/139‑1 (MQ=255)                                                                gCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGGAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACaaa                                                             >  1:401218/1‑139 (MQ=255)                                                                   aaaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCg                                                          <  1:46054/139‑1 (MQ=255)                                                                     aGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGaa                                                        >  1:1637487/1‑139 (MQ=255)                                                                      gCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAAc                                                       <  2:275940/139‑1 (MQ=255)                                                                        atcatTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGGAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCa                                                     <  2:401218/139‑1 (MQ=255)                                                                                tCCTATCTCAGTTGGTATCAAAACGCTACGCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctc                                             <  1:1784090/139‑1 (MQ=255)                                                                                      ctcAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGa                                                                                                                             >  2:50011/1‑53 (MQ=255)                                                                                      ctcAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGa                                                                                                                             <  1:50011/53‑1 (MQ=255)                                                                                      ctcAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGGAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTc                                                                               >  2:1271148/1‑99 (MQ=255)                                                                                      ctcAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTc                                                                               <  1:1271148/99‑1 (MQ=255)                                                                                       tcAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAg                                      <  2:1637487/139‑1 (MQ=255)                                                                                           ttGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAg                                  >  1:5200/1‑139 (MQ=255)                                                                                            tGGTATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGt                                 <  1:1757781/139‑1 (MQ=255)                                                                                               tATCAAAACGCTACCCACTAATGGAGGGTAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATATTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACaaa                                                             >  2:1165232/1‑108 (MQ=255)                                                                                               tATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAGAGGAAAATTCTCAGCGcc                                                                                                          <  2:1238882/63‑1 (MQ=255)                                                                                               tATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGcc                                                                                                          >  1:1238882/1‑63 (MQ=255)                                                                                               tATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACaaa                                                             <  1:1165232/108‑1 (MQ=255)                                                                                               tATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAg                                      >  2:107773/1‑131 (MQ=255)                                                                                               tATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAg                                      <  1:107773/131‑1 (MQ=255)                                                                                               tATCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTa                              >  2:849827/1‑139 (MQ=255)                                                                                                aTCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCGTTTTCCTCTCCCCAGATAGTTTAg                             >  1:220373/1‑139 (MQ=255)                                                                                                aTCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTcc                                               >  2:973233/1‑121 (MQ=255)                                                                                                aTCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTcc                                               <  1:973233/121‑1 (MQ=255)                                                                                                aTCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTc                                                <  2:619121/120‑1 (MQ=255)                                                                                                aTCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTc                                                >  1:619121/1‑120 (MQ=255)                                                                                                aTCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGt                                 <  2:1110232/135‑1 (MQ=255)                                                                                                aTCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAg                             >  1:317991/1‑139 (MQ=255)                                                                                                aTCAAAACGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGAGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAAATAGt                                 >  1:1110232/1‑135 (MQ=255)                                                                                                   aaaaCGCTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGtt                          >  2:82566/1‑139 (MQ=255)                                                                                                         cTACCCACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAagag                    >  1:507707/1‑139 (MQ=255)                                                                                                            cccACTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATc                 <  1:1312523/139‑1 (MQ=255)                                                                                                               aCTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCggg              <  1:326244/139‑1 (MQ=255)                                                                                                               aCTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCggg              <  1:503948/139‑1 (MQ=255)                                                                                                               aCTAATGGAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCggg              <  2:1091372/139‑1 (MQ=255)                                                                                                                      gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTa                                                                                          <  2:1637375/56‑1 (MQ=255)                                                                                                                      gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTa                                                                                          <  2:218525/56‑1 (MQ=255)                                                                                                                      gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTa                                                                                          >  1:1637375/1‑56 (MQ=255)                                                                                                                      gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTa                                                                                          >  1:218525/1‑56 (MQ=255)                                                                                                                      gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctc                                             <  1:1487082/101‑1 (MQ=255)                                                                                                                      gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCctc                                             >  2:1487082/1‑101 (MQ=255)                                                                                                                      gAGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTg       >  2:91856/1‑139 (MQ=255)                                                                                                                       aGGGAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCGTTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGt      <  2:220373/139‑1 (MQ=255)                                                                                                                        gggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTa     >  2:621053/1‑139 (MQ=255)                                                                                                                         ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGa                                                                           >  1:1795290/1‑68 (MQ=255)                                                                                                                         ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGa                                                                           <  2:1795290/68‑1 (MQ=255)                                                                                                                         ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTcc                                          <  2:892419/101‑1 (MQ=255)                                                                                                                         ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTcc                                          >  1:892419/1‑101 (MQ=255)                                                                                                                         ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGa                                     <  2:789434/106‑1 (MQ=255)                                                                                                                         ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGa                                     >  1:789434/1‑106 (MQ=255)                                                                                                                         ggAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTaa    >  2:26158/1‑139 (MQ=255)                                                                                                                          gAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCt                                                                                           >  2:1880411/1‑51 (MQ=255)                                                                                                                          gAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCt                                                                                           <  1:1880411/51‑1 (MQ=255)                                                                                                                          gAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTaa    >  2:1094642/1‑138 (MQ=255)                                                                                                                          gAAACTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTaa    <  1:1094642/138‑1 (MQ=255)                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtgt                                                                                    <  2:1344794/57‑1 (MQ=255)                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAActgt                                                                                    >  1:1344794/1‑57 (MQ=255)                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGtt                                                                                   <  1:978954/58‑1 (MQ=255)                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGtt                                                                                   >  2:978954/1‑58 (MQ=255)                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAagca                                                                       <  2:1079209/70‑1 (MQ=255)                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAagca                                                                       >  1:1079209/1‑70 (MQ=255)                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAg                                      >  1:888363/1‑103 (MQ=255)                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAg                                      <  2:888363/103‑1 (MQ=255)                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGttt                               <  2:1555657/110‑1 (MQ=255)                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGttt                               >  1:1555657/1‑110 (MQ=255)                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCggg              >  2:675723/1‑127 (MQ=255)                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCggg              <  1:675723/127‑1 (MQ=255)                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATa  >  2:253758/1‑139 (MQ=255)                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATa  >  2:812011/1‑139 (MQ=255)                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATa  >  1:623092/1‑139 (MQ=255)                                                                                                                           aaaCTTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGAGGCGTTTGTAATa  >  1:1738904/1‑139 (MQ=255)                                                                                                                               |                                                                                                                                       ACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATA  >  NC_000913/3792700‑3792960

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