Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,407,931	G→T	12.1%	G426G (GGG→GGT)	rnr →	RNase R

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,407,931	0	G	T	12.1%	210.9 / 9.3	69	G426G (GGG→GGT)	rnr	RNase R
Reads supporting (aligned to +/- strand):  ref base G (26/34);  new base T (7/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 7.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAA  >  NC_000913/4407782‑4408080                                                                                                                                                      |                                                                                                                                                      tctgACGCGAGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTggg                                                                                                                                                       <  2:991274/150‑1 (MQ=255)    cACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGAtc                                                                                                                                                    <  1:1244606/150‑1 (MQ=255)           gTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTg                                                                                                                                             <  1:412741/150‑1 (MQ=255)               gACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTATATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGagag                                                                                                                                         >  1:868532/1‑150 (MQ=255)                       ccAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGaagaag                                                                                                                                 >  1:1371088/1‑150 (MQ=255)                       ccAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGaagaag                                                                                                                                 <  1:1464720/150‑1 (MQ=255)                       ccAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGaagaag                                                                                                                                 >  1:390597/1‑150 (MQ=255)                          aatgcTGGCATTTTCTGCAGGGCGATCATGATCTTGGCTAGCAGTACGCCCCGCTGGTTAAGGATCTCGAAGATTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGa                                                                                                                              <  1:708985/146‑1 (MQ=255)                                ggCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCa                                                                                                                        >  1:725783/1‑150 (MQ=255)                                  cATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtt                                                                                                                      <  2:97911/150‑1 (MQ=255)                                  cATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtt                                                                                                                      <  2:1370943/150‑1 (MQ=255)                                    tatTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtttt                                                                                                                    >  1:392833/1‑150 (MQ=255)                                        cTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAAc                                                                                                                <  1:802415/150‑1 (MQ=255)                                         tGCAGGGCGATCAGGATCTGCGCGAGCAGTACG‑CCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATACAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGc                                                                                                              >  1:1539508/1‑150 (MQ=255)                                          gCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGc                                                                                                              <  2:1267349/150‑1 (MQ=255)                                            aGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCAAAGAAGCGAAGTTCATTTTCAACGCTg                                                                                                            >  2:716001/1‑150 (MQ=255)                                             gggCGATCAGGATCTGCGCGATCAGTAAGCCCCGCTGGTTAAGCATCTCGAATAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGGTAATAACGCGGTGGGATCTCATTTTAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGa                                                                                                           <  2:1507276/150‑1 (MQ=255)                                             gggCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGa                                                                                                           >  1:1330978/1‑150 (MQ=255)                                               gCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAAc                                                                                                         <  2:725783/150‑1 (MQ=255)                                                cGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACg                                                                                                        >  1:1626063/1‑150 (MQ=255)                                                   tCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCg                                                                                                     <  2:1009207/150‑1 (MQ=255)                                                     aGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTa                                                                                                   >  1:154607/1‑150 (MQ=255)                                                       gATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTAtt                                                                                                 >  1:158163/1‑150 (MQ=255)                                                       gATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTAtt                                                                                                 <  2:1330978/150‑1 (MQ=255)                                                        aTCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACTCGGTGGTATCACAATAAATAAAAAAAAAAAGAAGTTCATTTTCAACGCTGAACGCCGTATTg                                                                                                >  2:791542/1‑150 (MQ=255)                                                        aTCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCTAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTg                                                                                                <  1:716001/150‑1 (MQ=255)                                                          cTGCGCGAACAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACTCGGTGGTATCACAATAAATAAAAAAAAAAAGAAGTTCATTTTCAACGCTGAACGCCGTATTGaa                                                                                              >  2:1318798/1‑150 (MQ=255)                                                           tGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAAc                                                                                             >  2:682775/1‑150 (MQ=255)                                                                  gCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGa                                                                                      >  2:443313/1‑150 (MQ=255)                                                                   cAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACTCGGTGGTATCACATTAAATAAAAAAAAAGAGAAGTTCATTTTCAACGATGAACGCCGTATTGAAAGTATCGaa                                                                                     >  2:826029/1‑150 (MQ=255)                                                                   cAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGaa                                                                                     <  2:773166/150‑1 (MQ=255)                                                                    aGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAAc                                                                                    <  2:1021396/150‑1 (MQ=255)                                                                    aGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAAc                                                                                    <  1:322704/150‑1 (MQ=255)                                                                     gTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACTCGGTGGTATCACATTTAATAGATAACAACCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACa                                                                                   >  2:22640/1‑150 (MQ=255)                                                                            ccGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAg                                                                            >  2:1114162/1‑150 (MQ=255)                                                                             cGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGc                                                                           <  2:372498/150‑1 (MQ=255)                                                                                    tAAGAATCTCGAAGAGTTGCATAACCTCTATAAATTGCTGGATAAAGCCCGTGAAGAACGCGGTGGTATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGATCGGCGTATTGAACGTATCGAACAGACCCAGCGTAacga                                                                    <  1:1632434/150‑1 (MQ=255)                                                                                      aGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAacgacg                                                                  >  2:390358/1‑150 (MQ=255)                                                                                       gCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGAcgc                                                                 <  2:1371088/150‑1 (MQ=255)                                                                                         aTCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACgcgc                                                               >  1:1286456/1‑150 (MQ=255)                                                                                           ctcGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCCGCGTAACGACGCGcac                                                             >  2:662204/1‑150 (MQ=255)                                                                                           ctcGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGcac                                                             <  2:1626063/150‑1 (MQ=255)                                                                                                     tGCATACCCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGTATCTCATTTGAGAGCGAAGAAGCGAAGTTCCTTTTCAACGCTGAACGCCGTATTGAACGTATCGAAAAGACCCAGCGTAACTACGCGCACAAATTAATTg                                                   <  1:1145752/150‑1 (MQ=255)                                                                                                     tGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTg                                                   <  2:392833/150‑1 (MQ=255)                                                                                                     tGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTg                                                   <  1:390358/150‑1 (MQ=255)                                                                                                     tGCATAACATCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGACGTATTGAACCAATCGAAAAGACCCAGAGTAACGAAGCGCACAAATTAATTg                                                   >  1:745561/1‑150 (MQ=255)                                                                                                      tcATAACCTCTATAAAGTGCTTGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTgt                                                  <  1:1318798/149‑2 (MQ=255)                                                                                                              tctATAAAGTGCTGGATAAAGCCCGTGAAGAACTCGGTGGTCTCACAATAAATAAAAAAAAAAAGAAGTTAATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCa                                          >  2:222286/1‑150 (MQ=255)                                                                                                                tataAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGAcgcg                                                                <  1:1656853/126‑1 (MQ=255)                                                                                                                tataAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGAcgcg                                                                >  2:1656853/1‑126 (MQ=255)                                                                                                                 ataAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGCCGCGCACAAATTATTGGAAGAGTGCATGa                                       >  1:750185/1‑150 (MQ=255)                                                                                                                 ataAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAACGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGa                                       >  2:1420737/1‑150 (MQ=255)                                                                                                                      gTGCTGGATAAAGCCCGTGTAGAACTCGGTGGTATCACAATAAATAAAAAAAAAAAAAAGTTAATTTTCAACGATGAACGACGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGAATGATTCTg                                  >  2:1311079/1‑150 (MQ=255)                                                                                                                        gCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGc                                <  2:1268389/150‑1 (MQ=255)                                                                                                                         cTGGATAAAGCCCGTGAAGATCTCGGTGGTATCACATTAAATATAAAAAAAACGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCg                               >  2:753252/1‑150 (MQ=255)                                                                                                                         cTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCg                               >  2:1431701/1‑150 (MQ=255)                                                                                                                         cTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCg                               <  1:1361885/150‑1 (MQ=255)                                                                                                                         cTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCg                               <  1:268962/150‑1 (MQ=255)                                                                                                                          tGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAagag                                              <  2:1409569/134‑1 (MQ=255)                                                                                                                          tGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAagag                                              >  1:1409569/1‑134 (MQ=255)                                                                                                                           ggATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAaattaat                                                     <  2:460381/126‑1 (MQ=255)                                                                                                                           ggATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAaattaat                                                     >  1:460381/1‑126 (MQ=255)                                                                                                                            gATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAat                            <  2:109419/150‑1 (MQ=255)                                                                                                                             aTAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAata                           >  1:626024/1‑150 (MQ=255)                                                                                                                             aTAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAata                           <  1:1431701/150‑1 (MQ=255)                                                                                                                                  gCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTcg                      <  2:839876/150‑1 (MQ=255)                                                                                                                                    ccGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTcggc                    <  2:868532/150‑1 (MQ=255)                                                                                                                                     cGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTcggcg                   >  2:1332192/1‑150 (MQ=255)                                                                                                                                                 ggTTGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCTTGATTCTGGCGAATATCTCGGCGGCGCGTTTCGtt       <  1:791542/150‑1 (MQ=255)                                                                                                                                                  gTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTg      <  1:1114162/150‑1 (MQ=255)                                                                                                                                                    gggATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTgag    >  1:554527/1‑150 (MQ=255)                                                                                                                                                    gggATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTgag    >  1:965248/1‑150 (MQ=255)                                                                                                                                                    gggATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTgag    <  1:826029/150‑1 (MQ=255)                                                                                                                                                      gATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGaa  >  2:1224618/1‑150 (MQ=255)                                                                                                                                                      gATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGaa  >  2:1106939/1‑150 (MQ=255)                                                                                                                                                      gATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGaa  >  2:1230503/1‑150 (MQ=255)                                                                                                                                                      gATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCAAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGaa  >  1:534364/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TCTCACGCGCGTCTGACCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAA  >  NC_000913/4407782‑4408080

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