Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	327,213	G→C	5.9%	S12* (TCA→TGA)	betA ←	choline dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	327,213	0	G	C	5.9%	321.9 / 6.9	95	S12* (TCA→TGA)	betA	choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (44/45);  new base C (4/2);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGCCGGAG  >  NC_000913/327065‑327360                                                                                                                                                     |                                                                                                                                                    gTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGa                                                                                                                                                    <  1:401924/150‑1 (MQ=255) gTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGa                                                                                                                                                    >  1:1123298/1‑150 (MQ=255)    gTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGGCAAAGCGATAGGCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGcc                                                                                                                                                 <  2:1114293/150‑1 (MQ=255)     tAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCAGCTCAGCCg                                                                                                                                                >  2:1314678/1‑150 (MQ=255)       gCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGc                                                                                                                                              >  1:943414/1‑150 (MQ=255)          tttACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCAcc                                                                                                                                           >  1:1171143/1‑150 (MQ=255)          tttACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCAGCTCAGCCGGCaac                                                                                                                                           <  1:1314678/150‑3 (MQ=255)           ttACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTCAGCCGGCaacg                                                                                                                                          <  2:544614/150‑4 (MQ=255)             aCCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaat                                                                                                                                        >  1:739047/1‑150 (MQ=255)              cccTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGGTAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaata                                                                                                                                       >  1:683313/1‑150 (MQ=255)              cccTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCtata                                                                                                                                       <  1:56024/150‑1 (MQ=255)              cccTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaata                                                                                                                                       <  1:904402/150‑1 (MQ=255)              cccTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaata                                                                                                                                       <  1:838306/150‑1 (MQ=255)               ccTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaataa                                                                                                                                      >  1:1249741/1‑150 (MQ=255)                cTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCaataat                                                                                                                                     <  1:774117/150‑1 (MQ=255)                  gTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAatga                                                                                                                                   <  2:1360924/150‑1 (MQ=255)                  gTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAatga                                                                                                                                   >  1:7705/1‑150 (MQ=255)                    aGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAatgatg                                                                                                                                 <  1:909332/150‑1 (MQ=255)                        ggAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGt                                                                                                                             <  2:1005414/150‑1 (MQ=255)                         gAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTc                                                                                                                            >  2:1293566/1‑150 (MQ=255)                           aTGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCaa                                                                                                                          >  2:297998/1‑150 (MQ=255)                            tGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCaaa                                                                                                                         >  1:1439757/1‑150 (MQ=255)                             gCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAAt                                                                                                                        >  2:898902/1‑150 (MQ=255)                             gCCACGGCAGCGGGCATCTGGGTGTGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAAt                                                                                                                        >  1:1073650/1‑150 (MQ=255)                              ccAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAAtt                                                                                                                       <  1:151699/150‑1 (MQ=255)                              ccAGGGCAGCGGGAATCTGGCTGCGGAAGTCAAAGCGATAGTCACGGCCGCCCGCTTCAAGAAGCAGCACGTAGGTAATCGGATCTTCAGTCAGACGGGTAGCGAGAACCTTGCAGGCTGAGCCGGCACGAATAATGATGATGTCATAtt                                                                                                                       >  2:1022365/1‑150 (MQ=255)                                aGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGGCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGGTGAGCGGGCACCAATTATTATGTAGTCAAATTGc                                                                                                                     <  2:332918/150‑1 (MQ=255)                                  ggCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCaa                                                                                                                   <  1:1497641/150‑1 (MQ=255)                                    cAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAAt                                                                                                                 >  1:515319/1‑150 (MQ=255)                                      gCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATaa                                                                                                               <  1:1111077/150‑1 (MQ=255)                                              cTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTg                                                                                                       <  2:1179049/150‑1 (MQ=255)                                               tGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTgg                                                                                                      >  2:964027/1‑150 (MQ=255)                                                 ggTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGtt                                                                                                    <  2:1123298/150‑1 (MQ=255)                                                  gTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCTTATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTa                                                                                                   <  1:1461474/150‑1 (MQ=255)                                                  gTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTa                                                                                                   >  2:702244/1‑150 (MQ=255)                                                   tGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAg                                                                                                  <  2:807624/150‑1 (MQ=255)                                                    gCGGAATGCAAATGGTTAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGa                                                                                                 <  2:283414/150‑1 (MQ=255)                                                     cGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCAcc                                                                                                                                           >  2:1071943/1‑107 (MQ=255)                                                     cGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCAcc                                                                                                                                           <  1:1071943/107‑1 (MQ=255)                                                     cGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGaa                                                                                                >  2:1028094/1‑150 (MQ=255)                                                      ggAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAat                                                                                               <  1:900106/150‑1 (MQ=255)                                                       gAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAata                                                                                              <  2:513593/150‑1 (MQ=255)                                                         aGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATg                                                                                            <  2:1530127/150‑1 (MQ=255)                                                          gTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATGTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATgg                                                                                           <  2:502286/150‑1 (MQ=255)                                                           tCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCATCAGCACGGATGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGTCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGa                                                                                          <  2:338769/150‑1 (MQ=255)                                                            cAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGAc                                                                                         >  1:958926/1‑150 (MQ=255)                                                            cAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGAc                                                                                         <  2:958926/150‑1 (MQ=255)                                                            cAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGCTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGAc                                                                                         >  2:173172/1‑150 (MQ=255)                                                              aaGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTg                                                                                       >  2:1297374/1‑150 (MQ=255)                                                                   aTAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAAtt                                                                                  <  2:943414/150‑1 (MQ=255)                                                                    tAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAAttt                                                                                 >  2:264367/1‑150 (MQ=255)                                                                    tAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGCAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAAttt                                                                                 <  2:1345234/150‑1 (MQ=255)                                                                      gTCAGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGAATTCGGATCTTCAGTCAGACGCGTAGCGAGAACGTTGCCGGATGAGCCGGCACCAATAATGATGTAGTAAAATTGCAAATAAACCTCATGGTTAGAATATGGACTGGAATTTac                                                                               >  2:133930/1‑149 (MQ=255)                                                                       tCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGc                                                                              >  1:26534/1‑150 (MQ=255)                                                                        ccGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGCGCCGGCCCCAATAATGATGTAGTCAAATTGCAAATAAGCCTCCTGGTTAGAATATGGACTGGAATTTAGcc                                                                             >  1:1297138/1‑150 (MQ=255)                                                                         cGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCa                                                                            <  2:1167401/150‑1 (MQ=255)                                                                            gccgccCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCAtct                                                                         >  1:1412686/1‑150 (MQ=255)                                                                             ccgcccgcTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCAtctc                                                                        <  1:1332554/150‑1 (MQ=255)                                                                                    cTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACAtgg                                                                 >  2:753790/1‑150 (MQ=255)                                                                                         agcagcagcACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtggatgga                                                            <  2:1249741/150‑1 (MQ=255)                                                                                                gcaCGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTgg                                                                                      >  1:910134/1‑117 (MQ=255)                                                                                                gcaCGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTgg                                                                                      <  2:910134/117‑1 (MQ=255)                                                                                                gcaCGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtggatg                                                              <  2:258934/141‑1 (MQ=255)                                                                                                gcaCGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCtggatg                                                              >  1:258934/1‑141 (MQ=255)                                                                                                  aCGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTg                                                   <  2:548790/150‑1 (MQ=255)                                                                                                    ggAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTggg                                                 <  2:459734/150‑1 (MQ=255)                                                                                                       ggTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGgtgt                                              >  1:580310/1‑150 (MQ=255)                                                                                                        gttttCGGAGCTTCTGTCAGAGGGGTATCGATAACGTTTCCGGCAGAGCCGGAACCAATAATGTTGTAGGCATATTGCTAATATACTTCTTGGTTAGAATTTGGACTGGAATTTAGTCATCTCATCCTTGATGGACTTCACCTGGGTGTa                                             <  2:264977/147‑1 (MQ=255)                                                                                                          aTTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAAc                                           >  1:1111263/1‑150 (MQ=255)                                                                                                          aTTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAAc                                           >  1:305628/1‑150 (MQ=255)                                                                                                            tCGGATCTTCAGTCAGACGGGTAGCAAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAActc                                         >  2:259486/1‑150 (MQ=255)                                                                                                             cGGATCTTCAGTCAGACGGTTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAActct                                        <  1:133930/150‑1 (MQ=255)                                                                                                                   ttCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCg                                  >  1:1359786/1‑150 (MQ=255)                                                                                                                    tcagtcagACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGt                                 >  2:1535461/1‑150 (MQ=255)                                                                                                                     cagtcagACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGtc                                <  2:7705/150‑1 (MQ=255)                                                                                                                     cagtcagACCGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGtc                                <  2:1101904/150‑1 (MQ=255)                                                                                                                       gtcagCCGGGTAGCGAGAACGTTGCCGGCTGCGCCGGCACCAATAATGATGTAGTCAAATAGCAAATAAACCTCCTGGTTTGAATATGGACTGGAATTTAGCCATCTCAAACTGGATGGACTTCACCTCGGTGTAACTCTGGAGCGtcat                              >  2:1471968/1‑150 (MQ=255)                                                                                                                         cagACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGtcatca                            >  1:7542/1‑150 (MQ=255)                                                                                                                         cagACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGtcatca                            >  1:807522/1‑150 (MQ=255)                                                                                                                             cGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGcc                        >  1:863694/1‑150 (MQ=255)                                                                                                                                  aGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTtct                   <  2:1412686/150‑1 (MQ=255)                                                                                                                                   gCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTtctc                  >  1:525885/1‑150 (MQ=255)                                                                                                                                   gCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTtctc                  >  1:193614/1‑150 (MQ=255)                                                                                                                                    cGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTcg                 >  2:770192/1‑150 (MQ=255)                                                                                                                                     gagaACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTcgc                <  2:1171143/150‑1 (MQ=255)                                                                                                                                      agaACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTcgcg               <  1:297998/150‑1 (MQ=255)                                                                                                                                      agaACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTcgcg               <  2:842375/150‑1 (MQ=255)                                                                                                                                       tattggtGCCGGCTCAGCCGGCAACAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTcgcgc              >  2:1236819/7‑149 (MQ=255)                                                                                                                                           ggtGCCGGCTCAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCaa          >  1:807361/3‑150 (MQ=255)                                                                                                                                            gtGCCGGCTCAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGGGTCATCACGCCGTTCTCGCGACCAAt         >  1:17226/2‑150 (MQ=255)                                                                                                                                               ccGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGcc      <  1:1293566/150‑1 (MQ=255)                                                                                                                                                 ggttGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGCCgg    >  1:1412092/4‑150 (MQ=255)                                                                                                                                                 ggCTGAGCCGGCACCAATAATGTTGTATTCAAATTTCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGCCgg    <  2:683313/150‑1 (MQ=255)                                                                                                                                                 ggCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGCCgg    >  2:91381/1‑150 (MQ=255)                                                                                                                                                   cTGAGCCGGCACCAATAATCATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGCCGGac  <  1:403327/150‑2 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCGATAGTCCGGGCCGCCCGCTTCAAGCAGCAGCACGGAGGTATTCGGATCTTCAGTCAGACGGGTAGCGAGAACGTTGCCGGCTGAGCCGGCACCAATAATGATGTAGTCAAATTGCAAATAAACCTCCTGGTTAGAATATGGACTGGAATTTAGCCATCTCAACCTGGATGGACTTCACCTGGGTGTAACTCTGGAGCGTCATCACGCCGTTCTCGCGACCAATGCCGGAG  >  NC_000913/327065‑327360

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: ‑