Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,495,322	G→T	A438S (GCC→TCC)	nirB →	nitrite reductase catalytic subunit NirB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,495,322	0	G	T	100.0%	302.7 / NA	85	A438S (GCC→TCC)	nirB	nitrite reductase catalytic subunit NirB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (33/52);  total (33/52)

AACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGTTAAC  >  NC_000913/3495185‑3495456                                                                                                                                          |                                                                                                                                       aaCTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTcc                                                                                                                                      <  2:638942/140‑1 (MQ=255)  aCTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCa                                                                                                                                     <  1:992694/140‑1 (MQ=255)  aCTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCa                                                                                                                                     <  1:355513/140‑1 (MQ=255)  aCTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCa                                                                                                                                     <  2:63616/140‑1 (MQ=255)   cTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCAt                                                                                                                                    >  1:1101368/1‑140 (MQ=255)      ccGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATcaa                                                                                                                                 >  1:602935/1‑140 (MQ=255)      ccGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATcaa                                                                                                                                 >  1:1579612/1‑140 (MQ=255)        ggAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATcaaca                                                                                                                               <  1:1437931/140‑1 (MQ=255)        ggAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATcaaca                                                                                                                               <  2:829497/140‑1 (MQ=255)          aaaaCCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACaaa                                                                                                                             <  2:1613164/140‑1 (MQ=255)          aaaaCCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACaaa                                                                                                                             <  2:1579612/140‑1 (MQ=255)            aaCCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAgg                                                                                                                           <  2:397403/140‑1 (MQ=255)             aCCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGc                                                                                                                          <  2:1101368/140‑1 (MQ=255)             aCCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGc                                                                                                                          <  2:914017/140‑1 (MQ=255)              cccGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACaaa                                                                                                                             <  1:1446693/136‑1 (MQ=255)              cccGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACaaa                                                                                                                             >  2:1446693/1‑136 (MQ=255)              cccGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCACCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCt                                                                                                                         >  1:841275/1‑140 (MQ=255)               ccGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTg                                                                                                                        <  1:1670194/140‑1 (MQ=255)                cGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGc                                                                                                                       >  1:923041/1‑140 (MQ=255)                cGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGc                                                                                                                       >  1:820946/1‑140 (MQ=255)                  gATTCCCTGATCCTGCCAGCACGCTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCca                                                                                                                     <  1:1082152/140‑1 (MQ=255)                     tCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCcacac                                                                                                                  <  2:1115876/140‑1 (MQ=255)                     tCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCcacac                                                                                                                  <  1:11735/140‑1 (MQ=255)                         tGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGtt                                                                                                              >  1:882512/1‑140 (MQ=255)                            tCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTgcg                                                                                                           <  2:343000/140‑1 (MQ=255)                                 ccAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCgct                                                                                                      >  1:70900/1‑140 (MQ=255)                                 ccAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCgct                                                                                                      >  1:1245004/1‑140 (MQ=255)                                 ccAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCgct                                                                                                      >  2:606737/1‑140 (MQ=255)                                   aGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCgctga                                                                                                    <  2:301394/140‑1 (MQ=255)                                     cacaCTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCgctga                                                                                                    <  1:1035414/138‑1 (MQ=255)                                     cacaCTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCgctga                                                                                                    >  2:1035414/1‑138 (MQ=255)                                     cacaCTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGCCTGCCACACAGTTGCGGCgctgaaa                                                                                                  >  1:666476/1‑140 (MQ=255)                                        acTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCgctgaaagct                                                                                               >  1:203747/1‑140 (MQ=255)                                          tCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCgctgaaagctga                                                                                             >  1:81207/1‑140 (MQ=255)                                              gTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAAcc                                                                                         <  1:624188/140‑1 (MQ=255)                                               tAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCa                                                                                        >  2:1063578/1‑140 (MQ=255)                                               tAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCa                                                                                        >  2:1647845/1‑140 (MQ=255)                                                  cGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAg                                                                                     <  2:670741/140‑1 (MQ=255)                                                   ggCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGc                                                                                    >  2:803035/1‑140 (MQ=255)                                                    gCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCg                                                                                   <  2:1252093/140‑1 (MQ=255)                                                    gCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCg                                                                                   >  1:78411/1‑140 (MQ=255)                                                    gCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCg                                                                                   >  1:1170101/1‑140 (MQ=255)                                                        gCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTa                                                                               >  1:599590/1‑140 (MQ=255)                                                        gCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTa                                                                               >  2:798045/1‑140 (MQ=255)                                                             cTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGc                                                                          <  2:1360529/140‑1 (MQ=255)                                                                   gtgtTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCggg                                                                    >  2:852722/1‑139 (MQ=255)                                                                   gtgtTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGt                                                                    >  1:903691/1‑140 (MQ=255)                                                                    tgtTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCgg                                                                     <  1:1198762/138‑1 (MQ=255)                                                                    tgtTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCgg                                                                     >  2:1198762/1‑138 (MQ=255)                                                                      ttGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGc                                                                 >  2:834050/1‑140 (MQ=255)                                                                        gATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTg                                                               <  1:1647845/140‑1 (MQ=255)                                                                              cTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATccc                                                         >  1:1015864/1‑140 (MQ=255)                                                                                 ccGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCt                                                      <  1:1006930/140‑1 (MQ=255)                                                                                    gACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGt                                                   <  2:1212896/140‑1 (MQ=255)                                                                                    gACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGt                                                   <  2:1170101/140‑1 (MQ=255)                                                                                    gACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGt                                                   <  2:923041/140‑1 (MQ=255)                                                                                    gACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGt                                                   <  1:939322/140‑1 (MQ=255)                                                                                    gACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGt                                                   <  1:606737/140‑1 (MQ=255)                                                                                    cccAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGATTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGt                                                   <  2:520040/138‑1 (MQ=255)                                                                                      cAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCa                                                 <  2:70900/140‑1 (MQ=255)                                                                                      cAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGACTGTATCCCGCTGGTCa                                                 <  2:1245004/140‑1 (MQ=255)                                                                                        gcgcgcAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACt                                               <  2:78411/140‑1 (MQ=255)                                                                                        gcgcgcAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACt                                               <  2:406213/140‑1 (MQ=255)                                                                                        gcacgcAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGTGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCGCt                                               <  2:1143211/137‑1 (MQ=255)                                                                                         cgcgcAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTc                                              >  2:981574/1‑140 (MQ=255)                                                                                         cgcgcAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTc                                              <  2:818991/140‑1 (MQ=255)                                                                                                     cTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAAcg                                  >  2:1661716/1‑140 (MQ=255)                                                                                                         tGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGa                              <  1:221147/140‑1 (MQ=255)                                                                                                         tGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGa                              <  2:841275/140‑1 (MQ=255)                                                                                                         tGCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGa                              <  2:81207/140‑1 (MQ=255)                                                                                                          gCTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGaa                             <  2:1361820/140‑1 (MQ=255)                                                                                                           cTTCGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAAc                            <  1:1063578/140‑1 (MQ=255)                                                                                                              cGACGTCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGTTGACACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTgg                         <  2:395923/140‑1 (MQ=255)                                                                                                                   tCACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGaaa                    <  1:213456/140‑1 (MQ=255)                                                                                                                    cACCAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAAc                   >  1:92192/1‑140 (MQ=255)                                                                                                                       cAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAgg                >  2:1428096/1‑140 (MQ=255)                                                                                                                       cAAAGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAgg                <  1:1576807/140‑1 (MQ=255)                                                                                                                          aGGTGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCa             <  2:1358134/140‑1 (MQ=255)                                                                                                                             tGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCg          <  1:1046990/140‑1 (MQ=255)                                                                                                                             tGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCg          <  1:490445/140‑1 (MQ=255)                                                                                                                             tGATCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCg          <  1:803035/140‑1 (MQ=255)                                                                                                                               aTCTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGaa        >  2:989444/1‑140 (MQ=255)                                                                                                                                 cTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGt      >  1:1189620/1‑140 (MQ=255)                                                                                                                                 cTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGt      <  1:1661716/140‑1 (MQ=255)                                                                                                                                 cTGATTGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGt      <  1:52700/140‑1 (MQ=255)                                                                                                                                     ttGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGTTaac  <  2:820946/140‑1 (MQ=255)                                                                                                                                     ttGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGTTaac  <  1:156438/140‑1 (MQ=255)                                                                                                                                     ttGCTTCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGTTaac  >  2:933843/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                       AACTGCCGGAAAACCCGGATTCCCTGATCCTGCCAGCACACTCGGGTAGCGGCAAGCCGTCTATCGGTGTTGATAAACTGCCGGACAGCGCGCAAATCTGCTCCTGCTTCGACGTCACCAAAGGTGATCTGATTGCTGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGTTAAC  >  NC_000913/3495185‑3495456

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