Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BX571857	2,232,539	G→A	87.3%	noncoding (424/2923 nt)	unknown ←

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BX571857	2,232,539	0	G	A	87.3%	94.7 / 16.8	55	noncoding (424/2923 nt)	unknown
Reads supporting (aligned to +/- strand):  ref base G (2/5);  new base A (17/31);  total (19/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

TAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCATCTTTCCAGATGATTCGTCTAATGT  >  BX571857/2232495‑2232631                                             |                                                                                             tAGGAGATGGTCCTCCCAGATTCCGGCGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATc                                          <  2:687003/97‑1 (MQ=9)  aGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCa                                         <  1:871132/97‑1 (MQ=17)  aGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCa                                         <  1:702907/97‑1 (MQ=17)   ggAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCGTCGTACTCAGGATCCACTCAAGAGAGCCAACATTTTCGACTACAGGATTCTCAc                                        >  2:1071908/1‑97 (MQ=12)   ggAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCAc                                        <  2:990284/97‑1 (MQ=18)   ggAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCAc                                        <  2:490270/97‑1 (MQ=18)   ggAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCAc                                        <  1:1223905/97‑1 (MQ=18)    gagaTGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCAcc                                       <  1:1082976/97‑1 (MQ=18)     agaTGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACct                                      >  2:579389/1‑97 (MQ=21)     agaTGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACct                                      >  2:1414485/1‑97 (MQ=21)      gaTGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACctt                                     <  1:11988/97‑1 (MQ=34)      gaTGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACctt                                     <  1:1269564/97‑1 (MQ=21)      gaTGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACctt                                     <  1:264772/97‑1 (MQ=21)      gaTGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACctt                                     <  1:579389/97‑1 (MQ=21)       aTGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACcttc                                    >  2:560478/1‑97 (MQ=21)       aTGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACcttc                                    <  2:437716/97‑1 (MQ=21)       aTGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACcttc                                    >  1:1211815/1‑97 (MQ=21)       aTGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACcttc                                    <  2:664436/97‑1 (MQ=18)       aTGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACcttc                                    <  1:108710/97‑1 (MQ=21)        tGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACcttct                                   >  1:112867/1‑97 (MQ=21)        tGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACcttct                                   <  2:852443/97‑1 (MQ=21)         ggTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACcttctt                                  <  1:1212872/97‑1 (MQ=21)         ggTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACcttctt                                  >  1:894686/1‑97 (MQ=18)         ggTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACcttctt                                  >  2:1140274/1‑97 (MQ=18)          gTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCTCCTTCttt                                 <  1:182680/97‑1 (MQ=18)          gTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCttt                                 >  1:278630/1‑97 (MQ=34)           tcctccCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTg                                <  2:314294/97‑1 (MQ=18)            cctccCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGa                               >  1:496760/1‑97 (MQ=21)            cctccCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGa                               <  1:67881/97‑1 (MQ=21)             ctccCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGAt                              <  1:26825/97‑1 (MQ=21)             ctccCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGAt                              >  2:797262/1‑97 (MQ=21)             ctccCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGAt                              >  2:263228/1‑97 (MQ=21)              tccCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGAtt                             <  1:1080216/97‑1 (MQ=21)              tccCAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGAtt                             <  2:1038565/97‑1 (MQ=21)                ccAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCa                           <  1:1139888/97‑1 (MQ=21)                 cAGATTCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCAt                          <  1:1410109/97‑1 (MQ=21)                    aTTCCGACGGAATTTCACGTGCTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCATCtt                       <  1:540911/97‑1 (MQ=34)                     ttCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCATCttt                      <  2:975426/97‑1 (MQ=21)                     ttCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCATCttt                      >  2:58142/1‑97 (MQ=18)                      tCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCATCTTTc                     <  2:744497/97‑1 (MQ=21)                      tCCGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCATCTTTc                     >  1:713095/1‑97 (MQ=18)                       ccGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCATCTTTcc                    >  1:1394029/1‑97 (MQ=21)                       ccGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCATCTTTcc                    >  2:410719/1‑97 (MQ=18)                        cGACGGAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGGGGGACAACATTTTCGACTGCAGGATTATCCCCTTCTTTGATTCATCTTTCCa                   <  2:677170/97‑1 (MQ=11)                           cGGAATTTCACGTGCTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCATCTTTCCAgat                <  1:336398/97‑1 (MQ=34)                            ggAGTTTCACGTGCTCCATCGTGTTCAGGTTCCACTCAAGAGAGACAACATTTTCGCCTGCAGGATTATCGCCCTCTTTGATTCATCTTTCCAgatg               <  2:57781/97‑1 (MQ=11)                            ggAATTTCACGTGCTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCATCTTTCCAgatg               <  2:1178211/97‑1 (MQ=33)                            ggAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCATCTTTCCAgatg               <  2:171538/97‑1 (MQ=21)                             gAATTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCATCTTTCCAgatga              <  1:37550/97‑1 (MQ=21)                              aaTTTCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCATCTTTCCAgatgat             <  2:1327966/97‑1 (MQ=21)                                tttCACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATCATCCCCTTCTTTGATTCATCCTTCCAGATGATTc           >  2:764596/1‑97 (MQ=11)                                   cACGTGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCATCTTTCCAGATGATTCGTc        >  2:315337/1‑97 (MQ=18)                                       tGCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCATCTTTCCAGATGATTCGTCTAAt    <  1:575948/97‑1 (MQ=21)                                        gCTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCATCTTTCCAGATGATTCGTCTAATg   <  1:233407/97‑1 (MQ=18)                                         cTCCATCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGCTTATCACCTTCTTTGATTCATCTTTCCAGATGATTCGTCTAATgt  >  2:1240246/1‑97 (MQ=255)                                             |                                                                                             TAGGAGATGGTCCTCCCAGATTCCGACGGAATTTCACGTGCTCCGTCGTACTCAGGATCCACTCAAGAGAGACAACATTTTCGACTACAGGATTATCACCTTCTTTGATTCATCTTTCCAGATGATTCGTCTAATGT  >  BX571857/2232495‑2232631

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑