Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,695,244	T→C	L100L (CTA→CTG)	purL ←	phosphoribosylformyl‑glycineamide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,695,244	0	T	C	95.3%	220.0 / 1.3	64	L100L (CTA→CTG)	purL	phosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (29/32);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ACAACTGCTCTGCATCATCTAAAGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGTAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGC  >  NC_000913/2695096‑2695391                                                                                                                                                     |                                                                                                                                                    acaacTGCTCTGCATCATCTAAAGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAg                                                                                                                                                   <  2:733085/151‑1 (MQ=255)  caacTGCTCTGCATCATCTAAAGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGc                                                                                                                                                  >  1:284178/1‑151 (MQ=255)   aacTGCTCTGCATCATCTAAAGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGcc                                                                                                                                                 >  1:282269/1‑151 (MQ=255)       gCTCTGCATCATCTAAAGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGc                                                                                                                                              <  2:389357/150‑1 (MQ=255)       gCTCTGCATCATCTAAAGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGc                                                                                                                                              <  2:6764/150‑1 (MQ=255)         tctGCATCATCTAAAGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGt                                                                                                                                           >  2:983183/1‑151 (MQ=255)            gCATCATCTAAAGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTGCCTGTTGCAGCCCGCAGTTa                                                                                                                                         <  1:442151/150‑1 (MQ=255)            gCATCATCTAAAGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGtt                                                                                                                                          <  2:891868/149‑1 (MQ=255)            gCATCATCTAAAGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGtt                                                                                                                                          <  1:24603/149‑1 (MQ=255)            gCATCATCTAAAGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTAt                                                                                                                                        >  1:1107667/1‑151 (MQ=255)              atcatcTAAAGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATgg                                                                                                                                      <  2:401671/151‑1 (MQ=255)              atcatcTAAAGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATgg                                                                                                                                      >  1:1104598/1‑151 (MQ=255)                 atcTAAAGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGc                                                                                                                                    >  1:373492/1‑150 (MQ=255)                   cTAAAGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAat                                                                                                                                 >  2:526781/1‑151 (MQ=255)                     aaaGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAatat                                                                                                                               <  1:925116/151‑1 (MQ=255)                     aaaGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAata                                                                                                                                <  2:795790/150‑1 (MQ=255)                            aaaaGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGc                                                                                                                        >  1:1174550/1‑151 (MQ=255)                              aaGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCtt                                                                                                                      <  2:185556/151‑1 (MQ=255)                               aGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCttt                                                                                                                     >  2:1179037/1‑151 (MQ=255)                                gACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTc                                                                                                                    >  1:1003257/1‑151 (MQ=255)                                      ttCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGa                                                                                                               <  1:568705/150‑1 (MQ=255)                                       tCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGAc                                                                                                              >  1:432063/1‑150 (MQ=255)                                          atcatGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAgg                                                                                                          <  2:109580/151‑1 (MQ=255)                                            catGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGga                                                                                                        <  2:198172/151‑1 (MQ=255)                                                  gTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATg                                                                                                  <  1:457405/151‑1 (MQ=255)                                                     gTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGTAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTg                                                                                               >  2:526949/1‑151 (MQ=255)                                                        cAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGcc                                                                                             >  1:84862/1‑150 (MQ=255)                                                        cAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCa                                                                                            <  2:1006975/151‑1 (MQ=255)                                                            tCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGa                                                                                         <  2:1126747/150‑1 (MQ=255)                                                                  gTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCgggg                                                                                  <  2:500364/151‑1 (MQ=255)                                                                     aCCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCa                                                                               <  2:1083165/151‑1 (MQ=255)                                                                         gctgcCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACca                                                                            <  2:1104626/150‑1 (MQ=255)                                                                         gctgcCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACca                                                                            <  2:1154132/150‑1 (MQ=255)                                                                           tgcCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGTAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACcagc                                                                          >  2:137429/1‑150 (MQ=255)                                                                           tgcCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACcagc                                                                          >  1:633575/1‑150 (MQ=255)                                                                            gcCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACcagcag                                                                        <  2:754466/151‑1 (MQ=255)                                                                            gcCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACcagca                                                                         <  2:70190/150‑1 (MQ=255)                                                                            gcCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACcagca                                                                         <  1:771343/150‑1 (MQ=255)                                                                              cATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAgg                                                                       >  1:179095/1‑150 (MQ=255)                                                                              cATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAgg                                                                       >  1:260840/1‑150 (MQ=255)                                                                              cATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAgg                                                                       >  1:188882/1‑150 (MQ=255)                                                                                  gTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGttt                                                                  <  1:412077/151‑1 (MQ=255)                                                                                     cATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGcc                                                               >  2:386699/1‑151 (MQ=255)                                                                                          gTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGTAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCtt                                                             <  2:1009153/148‑1 (MQ=255)                                                                                          gTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCg                                                          >  2:254023/1‑151 (MQ=255)                                                                                            cAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCggg                                                        >  2:346790/1‑151 (MQ=255)                                                                                                gTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGt                                                    <  2:942120/151‑1 (MQ=255)                                                                                                  aCCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGtgg                                                  <  2:1184323/151‑1 (MQ=255)                                                                                                    cGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGtggct                                                <  2:522112/151‑1 (MQ=255)                                                                                                         tCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGtggctggc                                             <  1:1046850/149‑1 (MQ=255)                                                                                                         tCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCg                                            <  2:970662/150‑1 (MQ=255)                                                                                                             tataGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGc                                       >  1:875645/1‑151 (MQ=255)                                                                                                             tataGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGc                                       >  2:875602/1‑151 (MQ=255)                                                                                                              ataGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGcc                                      >  1:102557/1‑151 (MQ=255)                                                                                                              ataGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGcc                                      >  1:183006/1‑151 (MQ=255)                                                                                                                     gCAACGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGGGGCTGGCGAGTGCCGGGCCa                                >  1:645918/1‑150 (MQ=255)                                                                                                                        aCGCCGCGCTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATAtt                            >  1:679198/1‑151 (MQ=255)                                                                                                                             gcgcTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTcagc                       >  2:728505/1‑151 (MQ=255)                                                                                                                              cgcTCAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTcagca                      <  2:163997/151‑1 (MQ=255)                                                                                                                                  cAAGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGc                   <  2:855819/150‑1 (MQ=255)                                                                                                                                   aaGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGt                 >  1:1067709/1‑151 (MQ=255)                                                                                                                                   aaGGCGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGt                 >  1:907797/1‑151 (MQ=255)                                                                                                                                       cGGTTTACCTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCaa             <  2:252738/151‑1 (MQ=255)                                                                                                                                              ccTGTTGCAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTgtgc      >  2:973962/1‑151 (MQ=255)                                                                                                                                                  ttgcAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTgtgcgtgc  >  1:787074/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    ACAACTGCTCTGCATCATCTAAAGCAAAAAAGACCGTTTCCATCATGCGGTCGTGCAGTTCAGCGGTAACCTGCTGCCATTGTTCATTGGTCAGCGTACCGGCTTCTATATAGTAAGCAACGCCGCGCTCAAGGCGGTTTACCTGTTGTAGCCCGCAGTTATGGGCAATATCGGTCGCTTTCGAAGACCAGGGAGAGATGGTGCCAGGACGCGGGGTCACCAGCAGGAGTTTGCCTTGCGGGGCGTGGCTGGCGAGTGCCGGGCCATATTTCAGCAGGCGTTCAAGTTGTGCGTGC  >  NC_000913/2695096‑2695391

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑