Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	186,439	C→T	A738T (GCA→ACA)	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	186,439	0	C	T	100.0%	280.8 / NA	70	A738T (GCA→ACA)	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/32);  total (38/32)

GCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTAAATCATG  >  NC_000913/186292‑186583                                                                                                                                                    |                                                                                                                                                 gCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGtgtc                                                                                                                                               <  1:57566/151‑1 (MQ=255)      ggCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGa                                                                                                                                           >  2:776272/1‑150 (MQ=255)        cTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAac                                                                                                                                         >  2:198390/1‑150 (MQ=255)          gCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACt                                                                                                                                      <  1:209842/151‑1 (MQ=255)            cAGCTACGTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTa                                                                                                                                    <  1:234083/151‑1 (MQ=255)                    ttGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTgc                                                                                                                              <  2:834446/149‑1 (MQ=255)                     tGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTgcggc                                                                                                                           >  2:241998/1‑151 (MQ=255)                     tGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTgcgg                                                                                                                            <  2:105355/150‑1 (MQ=255)                      gCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTgcggc                                                                                                                           <  2:828617/150‑1 (MQ=255)                          cAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAGATTGCGGCGGTCt                                                                                                                      <  1:386277/151‑1 (MQ=255)                          cAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCt                                                                                                                      >  2:396199/1‑151 (MQ=255)                              aCTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAAt                                                                                                                   <  2:1063928/150‑1 (MQ=255)                                 tGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTc                                                                                                                 <  1:1057607/149‑1 (MQ=255)                                 tGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCgg                                                                                                               >  1:437866/1‑151 (MQ=255)                                  gCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCgg                                                                                                               <  2:715942/150‑1 (MQ=255)                                     ccAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGcaca                                                                                                           <  1:661282/151‑1 (MQ=255)                                     ccAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGcac                                                                                                            >  2:643105/1‑150 (MQ=255)                                     ccAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGcac                                                                                                            >  2:498151/1‑150 (MQ=255)                                      cAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGcac                                                                                                            >  2:984769/1‑149 (MQ=255)                                      cAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGcac                                                                                                            >  1:876703/1‑149 (MQ=255)                                        gACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGAc                                                                                                        <  2:1171766/151‑1 (MQ=255)                                        gACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGAc                                                                                                        <  2:1142931/151‑1 (MQ=255)                                            aaaCCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCCCAGAcggcg                                                                                                    >  1:658982/1‑151 (MQ=255)                                            aaaCCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGAcggcg                                                                                                    >  2:844087/1‑151 (MQ=255)                                             aaCCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGAcggcgg                                                                                                   <  2:492836/151‑1 (MQ=255)                                              aCCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGAcggcgg                                                                                                   >  1:985382/1‑150 (MQ=255)                                                cGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCa                                                                                                 >  1:257156/1‑150 (MQ=255)                                                      aCCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGAt                                                                                          <  2:463019/151‑1 (MQ=255)                                                                gATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCgg                                                                                  >  2:966007/1‑149 (MQ=255)                                                                 aTCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTcc                                                                               >  1:970474/1‑151 (MQ=255)                                                                    tGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCggg                                                                            >  2:979895/1‑151 (MQ=255)                                                                     gCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCggg                                                                            <  1:465049/150‑1 (MQ=255)                                                                      cGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCt                                                                          <  1:394071/151‑1 (MQ=255)                                                                       ggACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTc                                                                         >  2:1068277/1‑151 (MQ=255)                                                                         aCAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTcc                                                                        >  2:1094101/1‑150 (MQ=255)                                                                                  tGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATaaaaa                                                               <  2:464965/150‑1 (MQ=255)                                                                                      aTCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCg                                                          >  2:741411/1‑151 (MQ=255)                                                                                       tCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCgg                                                         >  1:875245/1‑151 (MQ=255)                                                                                        cGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCgg                                                         <  2:1162802/150‑1 (MQ=255)                                                                                         gggTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTg                                                       <  2:986129/151‑1 (MQ=255)                                                                                         gggTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTg                                                       >  1:158277/1‑151 (MQ=255)                                                                                           gTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGcc                                                     >  1:717253/1‑151 (MQ=255)                                                                                             tCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCt                                                    <  2:338052/150‑1 (MQ=255)                                                                                               cAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTcc                                                  >  2:1060411/1‑150 (MQ=255)                                                                                               cAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCa                                                 <  2:93820/151‑1 (MQ=255)                                                                                                    cGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCAcgcg                                             >  2:621160/1‑150 (MQ=255)                                                                                                         aaGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCt                                       <  1:318321/151‑1 (MQ=255)                                                                                                          aGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTg                                      >  2:111510/1‑151 (MQ=255)                                                                                                              aTCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCgg                                   <  1:741401/150‑1 (MQ=255)                                                                                                                ccATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGc                                 <  1:97706/150‑1 (MQ=255)                                                                                                                  aTCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTa                               >  1:9813/1‑150 (MQ=255)                                                                                                                  aTCGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTa                               >  1:626443/1‑150 (MQ=255)                                                                                                                    cGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCCCAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGACTGCGGGCtga                              >  1:449593/1‑147 (MQ=255)                                                                                                                    cGCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAg                             <  2:851438/150‑1 (MQ=255)                                                                                                                     gCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGc                            >  2:163327/1‑150 (MQ=255)                                                                                                                     gCCATACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGc                            >  2:72823/1‑150 (MQ=255)                                                                                                                         tACCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATAc                       >  1:1015401/1‑151 (MQ=255)                                                                                                                          aCCGTCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATAcc                      >  2:238899/1‑151 (MQ=255)                                                                                                                              tCGCGAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCg                  >  2:571683/1‑151 (MQ=255)                                                                                                                               cgcgAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCgg                 <  2:624460/151‑1 (MQ=255)                                                                                                                                 cgAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGGGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGAAATAACAGCg                  >  1:263470/1‑148 (MQ=255)                                                                                                                                 cgAGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGtt               <  1:925158/151‑1 (MQ=255)                                                                                                                                   aGTGGTGAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTg             <  2:6867/151‑1 (MQ=255)                                                                                                                                         gAAAATTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTaa        <  2:202697/150‑1 (MQ=255)                                                                                                                                         gAAAATTGGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTaaa       <  1:931693/151‑1 (MQ=255)                                                                                                                                           aaaTTTGTGTGTCGTGAACACTTAAATTGCGGCGGTGTAATTCGGGACAGACGGCGGCAAACAGGTAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAACCTGCGGGCTAAGAAATAACAGCGGGTTGATTaaa       >  1:567860/1‑149 (MQ=255)                                                                                                                                           aaaTTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTAAAt      >  1:964977/1‑150 (MQ=255)                                                                                                                                            aaTTTGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTAAATCa    <  2:333488/151‑1 (MQ=255)                                                                                                                                              tttGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTAAATCAt   >  1:1102692/1‑150 (MQ=255)                                                                                                                                              tttGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTAAATCATg  >  2:1114845/1‑151 (MQ=255)                                                                                                                                               ttGTGTGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTAAATCATg  <  1:967118/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 GCGGTGGCTGCCAGCTACTTTGCGTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTAAATCATG  >  NC_000913/186292‑186583

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