Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	127,606	C→G	8.7%	intergenic (+19/‑306)	aceF → / → lpd	pyruvate dehydrogenase, E2 subunit/lipoamide dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	127,606	0	C	G	8.7%	217.7 / 8.7	70	intergenic (+19/‑306)	aceF/lpd	pyruvate dehydrogenase, E2 subunit/lipoamide dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (36/27);  new base G (4/2);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCAC  >  NC_000913/127458‑127755                                                                                                                                                     |                                                                                                                                                      gTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAAcg                                                                                                                                                      >  1:324171/1‑150 (MQ=255)  ttCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAAcgg                                                                                                                                                     <  2:32102/150‑1 (MQ=255)    cGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAAcggcc                                                                                                                                                   <  1:588675/150‑1 (MQ=255)    cGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAAcggcc                                                                                                                                                   <  1:151610/150‑1 (MQ=255)        ccGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCt                                                                                                                                               <  1:955908/150‑1 (MQ=255)         cgcgTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCtt                                                                                                                                              >  1:66352/1‑150 (MQ=255)         cgcgTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCtt                                                                                                                                              <  2:18339/150‑1 (MQ=255)         cgcgTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAAAGCCGGCCGTTGGGCCGGCtc                                                                                                                                              >  2:884558/1‑149 (MQ=255)             gcTGATGCTGCCGATTTCTCTCTCCTTCGCCCCCCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAAAGCCGGCCGTTGGGCCGGCtctttt                                                                                                                                          <  1:1273455/149‑1 (MQ=255)              cTGATGGTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGACGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTc                                                                                                                                         <  2:142607/150‑1 (MQ=255)                    ctgcCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAAAGCCGGCCGTTGGGCCGGCtcttttacttaca                                                                                                                                   <  1:884558/150‑8 (MQ=38)                      gcCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAAtc                                                                                                                                 >  1:907979/1‑150 (MQ=255)                        cGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAAtctc                                                                                                                               >  2:357505/1‑150 (MQ=255)                         gATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCa                                                                                                                              >  1:755267/1‑150 (MQ=255)                           tttCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATg                                                                                                                            <  2:276306/150‑1 (MQ=255)                             tctctctcCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGaa                                                                                                                          >  1:1415498/1‑150 (MQ=255)                             tctctctcCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGaa                                                                                                                          >  1:1410691/1‑150 (MQ=255)                               tctctcCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATg                                                                                                                        <  1:272102/150‑1 (MQ=255)                                 tctcCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTa                                                                                                                      >  1:1401653/1‑150 (MQ=255)                                  ctcCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTAt                                                                                                                     >  1:1407359/1‑150 (MQ=255)                                    ccTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTg                                                                                                                   <  1:1358838/150‑1 (MQ=255)                                      ttCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTAGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAg                                                                                                                 >  1:62899/1‑150 (MQ=255)                                           ccaccGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGttat                                                                                                            >  2:1339886/1‑150 (MQ=255)                                            caccGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGTttatt                                                                                                           <  2:79918/150‑1 (MQ=255)                                                   ggagcGACGGTGCTGATGGTGCCCGTGTCATTACCATCATAAACCACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGGTATTAGCGAAt                                                                                                    <  2:1315736/146‑1 (MQ=255)                                                    gATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATa                                                                                                   >  2:1214669/1‑150 (MQ=255)                                                     aTCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAg                                                                                                  <  1:1103538/150‑1 (MQ=255)                                                          cGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACaaa                                                                                             <  2:66352/150‑1 (MQ=255)                                                               cTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGt                                                                                        >  1:717017/1‑150 (MQ=255)                                                                    ggTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCg                                                                                   >  1:854303/1‑150 (MQ=255)                                                                      tGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGtt                                                                                 >  2:1403639/1‑150 (MQ=255)                                                                         ccGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAAcgg                                                                                                                                                     >  1:1082243/1‑79 (MQ=255)                                                                         ccGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAAcgg                                                                                                                                                     <  2:1082243/79‑1 (MQ=255)                                                                         ccGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTttgt                                                                              <  1:945574/150‑1 (MQ=255)                                                                            tttCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTtgttgt                                                                           <  2:1415498/150‑1 (MQ=255)                                                                            tttCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTtgttgt                                                                           <  2:1410691/150‑1 (MQ=255)                                                                                 ttACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTaaa                                                                      <  1:1310861/150‑1 (MQ=255)                                                                                   aCCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTaaaaa                                                                    >  1:314253/1‑150 (MQ=255)                                                                                      atcatTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTg                                                                 <  1:720986/150‑1 (MQ=255)                                                                                      atcatTAACAACACGCTGTCTGACATTAGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAAttt                                                                 >  2:784472/1‑149 (MQ=255)                                                                                          ttAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTaa                                                             >  1:1386158/1‑150 (MQ=255)                                                                                          ttAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTaa                                                             <  1:164241/150‑1 (MQ=255)                                                                                            aacaacACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACa                                                           >  1:514041/1‑150 (MQ=255)                                                                                                acacGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAAttt                                                       >  2:1424876/1‑150 (MQ=255)                                                                                                  acGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTg                                                     <  1:126639/150‑1 (MQ=255)                                                                                                  acGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTg                                                     >  1:668364/1‑150 (MQ=255)                                                                                                               ttCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACg                                        >  2:407115/1‑150 (MQ=255)                                                                                                                tCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACgg                                       >  2:661792/1‑150 (MQ=255)                                                                                                                 caccGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGa                                      >  2:44352/3‑150 (MQ=255)                                                                                                                 cGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGa                                      >  2:71394/1‑150 (MQ=255)                                                                                                                    cggcTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAg                                   <  2:801438/147‑1 (MQ=255)                                                                                                                    cGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAg                                   <  1:661792/150‑1 (MQ=255)                                                                                                                        tGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACg                               >  2:1148102/1‑150 (MQ=255)                                                                                                                         ggTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGa                              >  2:1126316/1‑150 (MQ=255)                                                                                                                          gTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGAc                             <  2:755267/150‑1 (MQ=255)                                                                                                                          gTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGAc                             >  2:859972/1‑150 (MQ=255)                                                                                                                            gATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGAccc                           <  1:386930/150‑1 (MQ=255)                                                                                                                              tGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGTTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCgg                         <  1:1412940/150‑1 (MQ=255)                                                                                                                              tGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCgg                         >  1:299941/1‑150 (MQ=255)                                                                                                                               gtaagtaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggt                        <  1:1339886/150‑1 (MQ=255)                                                                                                                               gtaagtaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggt                        >  2:598960/1‑150 (MQ=255)                                                                                                                                 aagtaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtgg                      >  1:97651/1‑150 (MQ=255)                                                                                                                                  cagaaaaaaGCCGGCCGTTGGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggt                     >  2:255751/4‑150 (MQ=255)                                                                                                                                  agtaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTTAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATAGCGACGGATAGAACGACCCggtggt                     <  2:624748/150‑1 (MQ=255)                                                                                                                                  agtaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggt                     >  1:1121970/1‑150 (MQ=255)                                                                                                                                  agtaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggt                     <  2:493912/150‑1 (MQ=255)                                                                                                                                   agaaaaaaGCCGGCCGTTGGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggtg                    >  2:143497/3‑150 (MQ=255)                                                                                                                                    gaaaaaaGCCGGCCATTGGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggtgg                   >  2:1257169/2‑150 (MQ=255)                                                                                                                                      aaaGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGtt                 >  2:444327/1‑150 (MQ=255)                                                                                                                                      aaaGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGtt                 >  1:259078/1‑150 (MQ=255)                                                                                                                                       aaGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTa                >  2:461112/1‑150 (MQ=255)                                                                                                                                       aaGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTa                >  2:197485/1‑150 (MQ=255)                                                                                                                                          agCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAggg             <  1:197485/150‑1 (MQ=255)                                                                                                                                           gCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGt            >  2:1328322/1‑150 (MQ=255)                                                                                                                                                  cAATGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACtt     >  2:1422445/1‑150 (MQ=255)                                                                                                                                                    aCGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTca   <  2:514041/150‑1 (MQ=255)                                                                                                                                                     cggccggcTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGAt                                     <  2:588919/115‑1 (MQ=255)                                                                                                                                                     cggccggcTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGAt                                     >  1:588919/1‑115 (MQ=255)                                                                                                                                                     cggccggcTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTcac  >  2:464698/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCAC  >  NC_000913/127458‑127755

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