Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,720,772	T→C	R31R (CGT→CGC)	insJ →	IS150 transposase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,720,772	0	T	C	100.0%	273.7 / NA	76	R31R (CGT→CGC)	insJ	IS150 transposase A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (43/33);  total (43/33)

ACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCT  >  NC_000913/3720636‑3720901                                                                                                                                         |                                                                                                                                  aCTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGctt                                                                                                                                 <  2:769069/139‑1 (MQ=255)    gCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGctttcg                                                                                                                              >  2:1086904/1‑137 (MQ=255)    gCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTgg                                                                                                                              >  1:36877/1‑139 (MQ=255)     cACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGgt                                                                                                                             >  1:782940/1‑139 (MQ=255)     cACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGgt                                                                                                                             >  1:769448/1‑139 (MQ=255)        ccATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTc                                                                                                                          >  2:744866/1‑139 (MQ=255)          aTTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGgt                                                                                                                             <  2:38098/134‑1 (MQ=255)          aTTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGgt                                                                                                                             >  1:38098/1‑134 (MQ=255)                 tgGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAAcc                                                                                                                 <  1:418996/139‑1 (MQ=255)                  ggACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAAcc                                                                                                                 >  1:1208992/1‑138 (MQ=255)                  ggACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAAcc                                                                                                                 <  2:1208992/138‑1 (MQ=255)                      gATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                            <  1:1806/139‑1 (MQ=255)                        tGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCa                                                                                                          >  2:811150/1‑139 (MQ=255)                        tGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCa                                                                                                          >  1:73015/1‑139 (MQ=255)                            aTGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGAc                                                                                                      <  2:769448/139‑1 (MQ=255)                            aTGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGAc                                                                                                      <  1:820821/139‑1 (MQ=255)                                  tAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGt                                                                                                <  2:127166/139‑1 (MQ=255)                                  tAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGt                                                                                                >  1:133617/1‑139 (MQ=255)                                    gCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                            >  2:605997/1‑125 (MQ=255)                                    gCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                            <  1:605997/125‑1 (MQ=255)                                     ccccTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTTATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                            <  1:875690/124‑1 (MQ=255)                                     ccccTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                            >  2:875690/1‑124 (MQ=255)                                     ccccTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATggg                                                                                                 <  2:1219964/135‑1 (MQ=255)                                     ccccTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATggg                                                                                                 >  1:1219964/1‑135 (MQ=255)                                      cccTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGcc                                                                                            <  1:744866/139‑1 (MQ=255)                                       ccTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGccc                                                                                           >  2:211101/1‑139 (MQ=255)                                        cTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCct                                                                                          <  1:811150/139‑1 (MQ=255)                                          aaTATGTCAACGCCAAAATCCACTTTTAAAAAGCGCCTTGAAGTCGTGAATCACTACTGCACAACTGATGATGGTTACAGGATCATCTAGTCACGCGTTGGTGTCCCGCGAACAGAGGTCAGGACAAGGGTTGCCctct                                                                                        >  2:971664/1‑139 (MQ=255)                                          aaTATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                            >  2:733008/1‑119 (MQ=255)                                          aaTATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                            <  1:733008/119‑1 (MQ=255)                                          aaTATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCctct                                                                                        >  2:1090232/1‑139 (MQ=255)                                                tCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                            >  1:501612/1‑113 (MQ=255)                                                tCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                            <  2:501612/113‑1 (MQ=255)                                                   aaGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTccccc                                                                                                                      >  1:134181/1‑100 (MQ=255)                                                   aaGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTccccc                                                                                                                      <  2:134181/100‑1 (MQ=255)                                                         aaaTACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGgag                                                                         >  2:1002181/1‑139 (MQ=255)                                                          aaTACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCct                                                                                          <  1:250720/121‑1 (MQ=255)                                                          aaTACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCct                                                                                          >  2:250720/1‑121 (MQ=255)                                                                         aGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGcc                                                                                            >  1:564440/1‑104 (MQ=255)                                                                         aGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGcc                                                                                            <  2:564440/104‑1 (MQ=255)                                                                          gcgcCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTccc                                                        <  2:605067/139‑1 (MQ=255)                                                                           cgcCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCa                                                       <  1:175032/139‑1 (MQ=255)                                                                             ccTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGcc                                                                                            >  1:781182/1‑100 (MQ=255)                                                                             ccTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGcc                                                                                            <  2:781182/100‑1 (MQ=255)                                                                             ccTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCaaa                                                     >  2:588274/1‑139 (MQ=255)                                                                                        tGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGt                                          >  1:36412/1‑139 (MQ=255)                                                                                              actactTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGc                                    <  1:971664/139‑1 (MQ=255)                                                                                              actactTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGc                                    <  1:1086904/139‑1 (MQ=255)                                                                                              actactTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGc                                    <  1:211101/139‑1 (MQ=255)                                                                                               ctactTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACAt                                                                                                    <  2:155825/74‑1 (MQ=255)                                                                                               ctactTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACAt                                                                                                    >  1:155825/1‑74 (MQ=255)                                                                                                  ctTCACAACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGAt                                >  1:29271/1‑139 (MQ=255)                                                                                                     cacaACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAAc                                                                               >  1:1157054/1‑89 (MQ=255)                                                                                                     cacaACTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAAc                                                                               <  2:1157054/89‑1 (MQ=255)                                                                                                        aaCTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCagag                          >  1:1048821/1‑139 (MQ=255)                                                                                                        aaCTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCagag                          <  2:29271/139‑1 (MQ=255)                                                                                                         aCTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGt                         <  2:36412/139‑1 (MQ=255)                                                                                                          cTGATGATGGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGtt                        >  2:184186/1‑139 (MQ=255)                                                                                                             atgatgGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCg                     >  2:918664/1‑139 (MQ=255)                                                                                                             atgatgGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCg                     >  1:700901/1‑139 (MQ=255)                                                                                                             atgatgGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGAC‑‑‑GGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTaaa                                               >  2:539770/1‑110 (MQ=255)                                                                                                             atgatgGTTACAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGAC‑‑‑GGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTaaa                                               <  1:539770/110‑1 (MQ=255)                                                                                                                       cAGGATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCGGATCCAGAGTTGCGTATTAAGGTc           >  1:722423/1‑139 (MQ=255)                                                                                                                          gATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTg        <  2:1228843/139‑1 (MQ=255)                                                                                                                          gATCATCTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTg        <  2:1048821/139‑1 (MQ=255)                                                                                                                             catcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                            >  1:138623/1‑36 (MQ=255)                                                                                                                             catcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGt                                                                                                            <  2:138623/36‑1 (MQ=255)                                                                                                                             catcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGaaa     >  1:996377/1‑139 (MQ=255)                                                                                                                              atcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACAt                                                                             >  1:240854/1‑66 (MQ=255)                                                                                                                              atcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACAt                                                                             <  2:240854/66‑1 (MQ=255)                                                                                                                              atcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAg    >  1:1245469/1‑139 (MQ=255)                                                                                                                              atcTCGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAg    >  2:935856/1‑139 (MQ=255)                                                                                                                               tctcGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGc   >  2:802483/1‑139 (MQ=255)                                                                                                                                ctcGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCt  >  2:130910/1‑139 (MQ=255)                                                                                                                                ctcGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGCTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCt  >  1:138459/1‑139 (MQ=255)                                                                                                                                ctcGGCACGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATAGGTTGCCCTGTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCt  >  2:796914/1‑139 (MQ=255)                                                                                                                                      aCGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTa                <  2:317859/119‑1 (MQ=255)                                                                                                                                      aCGCTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTa                >  1:317859/1‑119 (MQ=255)                                                                                                                                         |                                                                                                                                  ACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGTGTCCCCCGAACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCT  >  NC_000913/3720636‑3720901

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