Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	127,604	A→T	7.6%	intergenic (+17/‑308)	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,604	0	A	T	7.6%	226.7 / 4.4	75	intergenic (+17/‑308)	aceF/lpd	pyruvate dehydrogenase, E2 subunit/lipoamide dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (39/30);  new base T (4/2);  total (43/32)
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

GTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTT  >  NC_000913/127458‑127752                                                                                                                                                   |                                                                                                                                                     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:1410691/1‑150 (MQ=255)                             tctctctcCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGaa                                                                                                                       >  1:1415498/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:859972/1‑150 (MQ=255)                                                                                                                          gTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGAc                          <  2:755267/150‑1 (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                  <  2:493912/150‑1 (MQ=255)                                                                                                                                  agtaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggt                  >  1:1121970/1‑150 (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)                                                                                                                                                   |                                                                                                                                                     GTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTT  >  NC_000913/127458‑127752

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