Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,312,170	G→C	7.2%	L194V (CTG→GTG)	ompC ←	outer membrane porin C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,312,170	0	G	C	7.2%	222.8 / 7.5	69	L194V (CTG→GTG)	ompC	outer membrane porin C
Reads supporting (aligned to +/- strand):  ref base G (38/26);  new base C (3/2);  total (41/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

CAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCA  >  NC_000913/2312032‑2312315                                                                                                                                           |                                                                                                                                                  taGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCACTGCGTCaaaac                                                                                                                                        <  2:1326335/149‑5 (MQ=255) cAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGAc                                                                                                                                        <  2:63703/150‑1 (MQ=255) cAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGAc                                                                                                                                        <  1:1411219/150‑1 (MQ=255) cAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGAc                                                                                                                                        >  2:716254/1‑150 (MQ=255) cAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCACTGCGTCaaaac                                                                                                                                        <  2:32844/150‑5 (MQ=255)  aGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGAcc                                                                                                                                       >  1:305622/1‑150 (MQ=255)  aGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGAcc                                                                                                                                       >  1:767968/1‑150 (MQ=255)   gCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCg                                                                                                                                      >  1:1296775/1‑150 (MQ=255)      cGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCgttg                                                                                                                                   <  2:441319/150‑1 (MQ=255)        gTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCgttgtt                                                                                                                                 >  1:1411539/1‑150 (MQ=255)             cGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTgttagtta                                                                                                                            <  2:272447/150‑1 (MQ=255)               ttACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACg                                                                                                                          <  2:1296775/150‑1 (MQ=255)                tACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGc                                                                                                                         >  1:565317/1‑150 (MQ=255)                 aCCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGcc                                                                                                                        <  2:917099/150‑1 (MQ=255)                  ccGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGGGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACCCCGTCGCCGTTTTGCCGCAGTGCGTCACGACCGTTGTTAGTTACGccc                                                                                                                       >  2:1362839/1‑149 (MQ=255)                   cGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCAc                                                                                                                      >  1:714310/1‑150 (MQ=255)                    gATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTgttagtt                                                                                                                             <  2:1280284/142‑1 (MQ=255)                    gATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTgttagtt                                                                                                                             >  1:1280284/1‑142 (MQ=255)                        tAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTa                                                                                                                 >  2:873270/1‑150 (MQ=255)                          agcagcGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTaaa                                                                                                               >  1:1239363/1‑150 (MQ=255)                                ggTGTTCTGATCATCAGTTCGTTTGGAGCTGGAGATCTCACCACCGATAACGAAACCTTCGTAATCATAAGTGATAGAAGCGCCGACGCCGTCTCCGTTTTGAGGCTGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCttg                                                                                                         <  2:107910/150‑2 (MQ=255)                                ggTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGAGCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGGCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTc                                                                                                         <  1:901713/150‑1 (MQ=255)                                          gCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATggg                                                                                               <  2:530915/150‑1 (MQ=255)                                              cAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGc                                                                                           >  2:480407/1‑150 (MQ=255)                                                gTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCg                                                                                         <  1:1100864/150‑1 (MQ=255)                                                   cGTTTGGAGCTGTAGATCGCACCACCGATACAGAAACCTTAGTAATCATAAGTGATAGAACCGCCGACGCAGTCGCAGTTTTGACGAAGTTCGTCACGACCGTTGTTAGATAAGCAACTAGTAAAGCCTTCAACAGATGGGTTGCCGttt                                                                                      >  2:1086795/1‑150 (MQ=255)                                                   cGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGttt                                                                                      <  2:686965/150‑1 (MQ=255)                                                    gTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGtttt                                                                                     <  2:1411539/150‑1 (MQ=255)                                                    gTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGtttt                                                                                     <  1:631998/150‑1 (MQ=255)                                                     tttGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGttttt                                                                                    >  1:887924/1‑150 (MQ=255)                                                         gAGCTGGAGATCGCACCACCTATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTAccc                                                                                <  1:458100/150‑1 (MQ=255)                                                          aGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTg                                                                                            <  1:1349890/137‑1 (MQ=255)                                                          aGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTg                                                                                            >  2:1349890/1‑137 (MQ=255)                                                           gCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTg                                                                              >  2:318711/1‑150 (MQ=255)                                                                 gaTCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTg                                                                        >  2:791929/1‑150 (MQ=255)                                                                       accaccGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGc                                                                  >  2:417975/1‑150 (MQ=255)                                                                             gATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                            >  2:1093850/1‑150 (MQ=255)                                                                              aTACCGAAACCTTCTTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTc                                                           <  1:689255/150‑1 (MQ=255)                                                                                 ccGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAgg                                                        >  1:1196323/1‑150 (MQ=255)                                                                                 ccGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAgg                                                        >  2:1199488/1‑150 (MQ=255)                                                                                  cGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTAACCTGGTACTGAACAGCAAAGTTCAGGc                                                       <  2:1002971/150‑1 (MQ=255)                                                                                             tAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGa                                            >  2:1474018/1‑150 (MQ=255)                                                                                                 cATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCga                                        >  2:986865/1‑150 (MQ=255)                                                                                                 cATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCga                                        <  2:565317/150‑1 (MQ=255)                                                                                                   tAAGTGATAGAACCGCCGCCGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaag                                      >  1:490135/1‑150 (MQ=255)                                                                                                    aaGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaaga                                     >  2:1561848/1‑150 (MQ=255)                                                                                                     aGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaagaa                                    >  2:1054450/1‑150 (MQ=255)                                                                                                        gATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTc                                 >  2:295012/1‑150 (MQ=255)                                                                                                          tAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAg                               >  1:1394594/1‑150 (MQ=255)                                                                                                            gAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGttt                                                                                      >  1:1325117/1‑93 (MQ=255)                                                                                                            gAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGttt                                                                                      <  2:1325117/93‑1 (MQ=255)                                                                                                            gAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtg                             >  2:1594963/1‑150 (MQ=255)                                                                                                              aCCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGtt                           >  2:438867/1‑150 (MQ=255)                                                                                                                cgccgACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTTCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTAc                         <  2:1134349/150‑1 (MQ=255)                                                                                                                   cgACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGt                      <  1:480407/150‑1 (MQ=255)                                                                                                                       gccgtcgccgtTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGt                  >  1:209929/1‑150 (MQ=255)                                                                                                                        ccgtcgccgtTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTc                 >  1:1518118/1‑150 (MQ=255)                                                                                                                        ccgtcgccgtTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTc                 >  1:343457/1‑150 (MQ=255)                                                                                                                          gtcgccgtTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgc               <  2:714310/150‑1 (MQ=255)                                                                                                                          gtcgccgtTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgc               <  1:716254/150‑1 (MQ=255)                                                                                                                            cgccgtTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGa             >  1:267709/1‑150 (MQ=255)                                                                                                                               ggtcgtGACGCACTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGc          >  1:1513730/6‑150 (MQ=255)                                                                                                                               ggtcgtGACGCACTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGc          >  1:68620/6‑150 (MQ=255)                                                                                                                                   gtGACGCACTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGtt      >  1:258708/2‑150 (MQ=255)                                                                                                                                    tGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTa     >  1:1252225/1‑150 (MQ=255)                                                                                                                                     gACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc    >  2:1206362/1‑150 (MQ=255)                                                                                                                                      aCGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAcc   <  1:1474018/150‑1 (MQ=255)                                                                                                                                      aCGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAcc   <  1:353914/150‑1 (MQ=255)                                                                                                                                       cGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCa  >  2:917693/1‑150 (MQ=255)                                                                                                                                           |                                                                                                                                                  CAGCACGGTCGCCGTTACCGATGTAAGCAGCGGTGTTCTGAGCATCAGTACGTTTGGAGCTGGAGATCGCACCACCGATACCGAAACCTTCGTAATCATAAGTGATAGAACCGCCGACGCCGTCGCCGTTTTGACGCAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCA  >  NC_000913/2312032‑2312315

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