Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	340,444	(A)5→4	coding (558/816 nt)	BW25113_RS01680 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	340,440	0	A	.	92.5%	322.9 / 28.4	93	coding (554/816 nt)	BW25113_RS01680	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (5/2);  new base . (39/47);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

GTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTA  >  NZ_CP009273/340352‑340525                                                                                         |                                                                                      gTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGaa                                                                                      >  1:583416/1‑90 (MQ=255) gTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGaa                                                                                      >  1:610207/1‑90 (MQ=255) gTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGaa                                                                                      >  2:366926/1‑90 (MQ=255)     gAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGaaaatg                                                                                  >  2:1340401/1‑88 (MQ=255)      aaCAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGaaaatgg                                                                                 >  2:1327440/1‑87 (MQ=255)         aGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCtt                                                                             <  2:405048/90‑1 (MQ=255)           aTTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTcg                                                                           >  1:150699/1‑90 (MQ=255)                aaaTGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGaaa                                                                      <  2:840432/90‑1 (MQ=255)                  aTGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATg                                                                    <  2:980278/90‑1 (MQ=255)                    gTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgtc                                                                  >  1:907424/1‑90 (MQ=255)                    gTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgtc                                                                  >  1:913766/1‑90 (MQ=255)                    gTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgtc                                                                  >  2:145673/1‑90 (MQ=255)                    gTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCT‑AAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgtcg                                                                 >  1:1237384/1‑90 (MQ=255)                      aTTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgtcgt                                                                >  2:932140/1‑90 (MQ=255)                       ttcttcTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgtcgtc                                                               <  1:461750/90‑1 (MQ=255)                       ttcttcTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgtcgtc                                                               <  2:22964/90‑1 (MQ=255)                       ttcttcTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgtcgtc                                                               <  2:283030/90‑1 (MQ=255)                           tcttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTAtt                                                           >  2:1453862/1‑90 (MQ=255)                           tcttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTAtt                                                           >  1:179747/1‑90 (MQ=255)                              tGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCg                                                        <  1:24804/90‑1 (MQ=255)                              tGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCg                                                        <  2:617916/90‑1 (MQ=255)                                aaTACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGga                                                      >  2:1437014/1‑90 (MQ=255)                                    ccATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGagaca                                                  >  2:540202/1‑90 (MQ=255)                                      atatGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgtcgtc                                                               >  1:935590/1‑75 (MQ=255)                                      atatGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgtcgtc                                                               <  2:935590/75‑1 (MQ=255)                                      atatGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacaga                                                >  2:1404824/1‑90 (MQ=255)                                      atatGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacaga                                                >  2:272934/1‑90 (MQ=255)                                      atatGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacaga                                                >  2:299813/1‑90 (MQ=255)                                      atatGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacaga                                                >  1:1028050/1‑90 (MQ=255)                                      atatGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacaga                                                >  2:981097/1‑90 (MQ=255)                                        atGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcg                                              >  1:1149468/1‑90 (MQ=255)                                         tGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgc                                             <  2:182898/90‑1 (MQ=255)                                          gCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgcg                                            <  1:1437014/90‑1 (MQ=255)                                          gCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgcg                                            <  1:1469446/90‑1 (MQ=255)                                          gCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgcg                                            <  1:540202/90‑1 (MQ=255)                                          gCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgcg                                            <  1:1453862/90‑1 (MQ=255)                                          gCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgcg                                            <  1:272934/90‑1 (MQ=255)                                          gCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgcg                                            <  2:1036708/90‑1 (MQ=255)                                            cTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgtcgtc                                                               >  2:1457215/1‑69 (MQ=255)                                            cTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATgtcgtc                                                               <  1:1457215/69‑1 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgcg                                            <  1:378005/86‑1 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgcg                                            >  2:378005/1‑86 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgcg                                            >  2:902389/1‑86 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgcg                                            <  1:902389/86‑1 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                        >  1:1409464/1‑90 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                        >  2:1359953/1‑90 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                        >  1:686446/1‑90 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                        >  1:273698/1‑90 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                        >  1:130381/1‑90 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                        >  1:111878/1‑90 (MQ=255)                                              ttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCa                                        >  1:386834/1‑90 (MQ=255)                                               tGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCAt                                       <  2:179747/90‑1 (MQ=255)                                                    agaATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGc                                  <  1:544798/90‑1 (MQ=255)                                                    agaATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGc                                  <  1:343176/90‑1 (MQ=255)                                                     gaATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCa                                 >  1:613650/1‑90 (MQ=255)                                                          ttCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGgcgc                            <  2:1058400/90‑1 (MQ=255)                                                           tCCCTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCt                           <  1:145673/90‑1 (MQ=255)                                                              cTAAAGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTgc                        <  2:150699/90‑1 (MQ=255)                                                                  aGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                    <  2:1038871/90‑1 (MQ=255)                                                                  aGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                    <  1:1359953/90‑1 (MQ=255)                                                                  aGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                    <  1:692879/90‑1 (MQ=255)                                                                  aGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                    <  1:302132/90‑1 (MQ=255)                                                                  aGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                    <  1:567316/90‑1 (MQ=255)                                                                  aGGGAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                    <  1:1086026/90‑1 (MQ=255)                                                                    ggAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATgg                  <  1:247310/90‑1 (MQ=255)                                                                     gAAAGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGa                 <  2:308389/90‑1 (MQ=255)                                                                       aaGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt               <  1:1404824/90‑1 (MQ=255)                                                                       aaGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt               <  2:454355/90‑1 (MQ=255)                                                                       aaGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt               <  2:1566606/90‑1 (MQ=255)                                                                       aaGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt               <  1:1340401/90‑1 (MQ=255)                                                                       aaGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt               <  1:1463115/90‑1 (MQ=255)                                                                       aaGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt               <  1:299813/90‑1 (MQ=255)                                                                       aaGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt               <  1:568031/90‑1 (MQ=255)                                                                       aaGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt               <  1:92166/90‑1 (MQ=255)                                                                       aaGAGTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGatt               <  1:981097/90‑1 (MQ=255)                                                                        agagTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTa                                     <  2:144400/68‑1 (MQ=255)                                                                        agagTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTa                                     >  1:144400/1‑68 (MQ=255)                                                                          agTGTATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattatt            >  2:711489/1‑90 (MQ=255)                                                                           gtgtATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTgcgc                      >  2:1310770/1‑79 (MQ=255)                                                                           gtgtATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTgcgc                      <  1:1310770/79‑1 (MQ=255)                                                                           gtgtATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattatta           >  1:686155/1‑90 (MQ=255)                                                                             gtATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCAt                                       >  1:1035758/1‑60 (MQ=255)                                                                             gtATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCAt                                       <  2:1035758/60‑1 (MQ=255)                                                                             gtATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattattatt         >  1:520365/1‑90 (MQ=255)                                                                             gtATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattattatt         >  1:237187/1‑90 (MQ=255)                                                                             gtATGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattattatt         >  2:1398688/1‑90 (MQ=255)                                                                               aTGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTac       <  2:454033/90‑1 (MQ=255)                                                                                tGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgc                                             >  1:1539720/1‑51 (MQ=255)                                                                                tGCATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGAcgc                                             <  2:1539720/51‑1 (MQ=255)                                                                                  cATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCAt                                       <  1:29253/55‑1 (MQ=255)                                                                                  cATGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCAt                                       >  2:29253/1‑55 (MQ=255)                                                                                   aTGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                    >  1:1382943/1‑74 (MQ=255)                                                                                   aTGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                    <  2:1382943/74‑1 (MQ=255)                                                                                   aTGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCGTGGATTATTATTACACtt   >  1:434258/1‑90 (MQ=255)                                                                                   aTGTGG‑AAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACtt   >  2:1579468/1‑90 (MQ=255)                                                                                     gtgGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTa  >  2:60134/1‑90 (MQ=255)                                                                                     gtgGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTa  >  2:56589/1‑90 (MQ=255)                                                                                     gtgGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTa  >  1:345640/1‑90 (MQ=255)                                                                                         |                                                                                      GTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTA  >  NZ_CP009273/340352‑340525

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