Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	20	7,334	(T)7→6	intergenic (–/–)	– / –	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	20	7,328	0	T	.	98.8%	376.4 / ‑2.1	81	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (36/44);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGATTTGTGGGAT  >  20/7228‑7423                                                                                                     |                                                                                                ccacaccaTAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACt                                                                                                 >  2:787788/1‑101 (MQ=255) ccacaccaTAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACt                                                                                                 >  2:546504/1‑101 (MQ=255) ccacaccaTAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACt                                                                                                 >  1:1475294/1‑101 (MQ=255) ccacaccaTAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACt                                                                                                 <  2:1445402/101‑1 (MQ=255)   acaccaTAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttt                                                                                               <  2:121998/101‑1 (MQ=255)   acaccaTAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttt                                                                                               <  1:645365/101‑1 (MQ=255)     accaTAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttttt                                                                                             <  2:1168613/101‑1 (MQ=255)       catagcatgaGGGTCGTTGTCTATTGATAATGTTATGCAAAATTGTTTCGCTCTAGTCTTCTGGAATAAGGAGATATTAGTAACCTAATTCCGCtttttta                                                                                           <  2:569148/101‑2 (MQ=255)       cataacataaTGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACtttttta                                                                                           >  1:940855/1‑100 (MQ=255)        ataacataaTGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttttttat                                                                                          <  2:256304/101‑3 (MQ=255)        ataacataaTGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttttttat                                                                                          >  1:1125901/1‑99 (MQ=255)        ataacataaTGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttttttat                                                                                          <  2:264750/101‑3 (MQ=255)        ataacataaTGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACttttttat                                                                                          >  1:158849/1‑99 (MQ=255)           acataaTGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTtattat                                                                                      >  1:633191/1‑101 (MQ=255)            cataaTGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATtata                                                                                     <  2:726183/101‑1 (MQ=255)            cataaTGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATtata                                                                                     <  1:1700657/101‑1 (MQ=255)            cataaTGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATtata                                                                                     <  2:1529961/101‑1 (MQ=255)            cataaTGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATtata                                                                                     <  2:753905/101‑1 (MQ=255)                aTGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAAt                                                                                 <  2:31913/101‑1 (MQ=255)                aTGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAAt                                                                                 >  2:1495886/1‑101 (MQ=255)                aTGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAAt                                                                                 <  1:560122/101‑1 (MQ=255)                 tGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATc                                                                                <  2:797779/101‑1 (MQ=255)                  ggTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCt                                                                               <  2:491050/101‑1 (MQ=255)                  ggTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCt                                                                               >  2:264783/1‑101 (MQ=255)                    tCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCttt                                                                             <  2:648793/101‑1 (MQ=255)                    tCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCttt                                                                             <  1:1395998/101‑1 (MQ=255)                    tCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCttt                                                                             <  1:1267223/101‑1 (MQ=255)                     cTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTa                                                                            >  1:507145/1‑101 (MQ=255)                     cTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTa                                                                            <  1:1651322/101‑1 (MQ=255)                        tgtgTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTAtct                                                                         >  2:780978/1‑101 (MQ=255)                         gtgtATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTAtctt                                                                        <  1:450835/101‑1 (MQ=255)                          tgtATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTAtcttc                                                                       <  1:187092/101‑1 (MQ=255)                              ttGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTg                                                                   <  2:1013897/101‑1 (MQ=255)                               tGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTTATTATAAAATCTTTATCTTCTTTg                                                                   <  2:458116/101‑1 (MQ=255)                                   aaTGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCaa                                                              <  1:544737/101‑1 (MQ=255)                                    aTGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAAc                                                             >  2:859468/1‑101 (MQ=255)                                    aTGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAAc                                                             >  1:1307108/1‑101 (MQ=255)                                     tGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAAcc                                                            >  2:462998/1‑101 (MQ=255)                                       ttATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTa                                                          >  2:471006/1‑101 (MQ=255)                                       ttATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTa                                                          >  2:770933/1‑101 (MQ=255)                                       ttATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTa                                                          <  1:859473/101‑1 (MQ=255)                                       ttATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTa                                                          >  2:324895/1‑101 (MQ=255)                                        tATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTAc                                                         >  2:1015853/1‑101 (MQ=255)                                          ttCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCa                                                       <  2:1307102/101‑1 (MQ=255)                                          ttCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCa                                                       <  2:633186/101‑1 (MQ=255)                                          ttCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCa                                                       <  1:471011/101‑1 (MQ=255)                                          ttCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCa                                                       <  2:1355875/101‑1 (MQ=255)                                              aaaTTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTGTAAAATCTTTATCTTCTTTTGTCAACCTACCAGTTa                                                   <  2:840985/101‑1 (MQ=255)                                               aaTTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTAc                                                  <  2:1249062/101‑1 (MQ=255)                                               aaTTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTAc                                                  <  2:1618554/101‑1 (MQ=255)                                               aaTTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTAc                                                  <  2:158849/101‑1 (MQ=255)                                               aaTTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTAc                                                  <  2:354170/101‑1 (MQ=255)                                               aaTTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTAc                                                  <  2:884335/101‑1 (MQ=255)                                                aTTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACt                                                 <  2:726797/101‑1 (MQ=255)                                                aTTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACt                                                 <  2:1430363/101‑1 (MQ=255)                                                aTTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACt                                                 <  2:831803/101‑1 (MQ=255)                                                         cTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGcc                                        >  1:389234/1‑101 (MQ=255)                                                          tCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCa                                       >  2:1030429/1‑101 (MQ=255)                                                          tCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCa                                       <  1:1263424/101‑1 (MQ=255)                                                           ccAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCaa                                      >  1:882328/1‑101 (MQ=255)                                                           ccAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCaa                                      >  2:252588/1‑101 (MQ=255)                                                            cAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAac                                     <  1:546509/101‑1 (MQ=255)                                                             aGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAaca                                    >  2:683699/1‑101 (MQ=255)                                                              gTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAacac                                   >  2:1528588/1‑101 (MQ=255)                                                              gTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAacac                                   <  1:286105/101‑1 (MQ=255)                                                              gTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAacac                                   >  1:1396003/1‑101 (MQ=255)                                                                  tctGGAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGcc                               <  2:1125896/101‑1 (MQ=255)                                                                     ggAATAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGa                            >  2:1037304/1‑101 (MQ=255)                                                                       atctaGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACg                          <  2:1660355/97‑1 (MQ=255)                                                                       aaTAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACg                          >  2:930399/1‑101 (MQ=255)                                                                       aaTAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACg                          >  1:1435706/1‑101 (MQ=255)                                                                        aTAAGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGa                         >  1:724192/1‑101 (MQ=255)                                                                          aaGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGaaa                       <  1:1623974/101‑1 (MQ=255)                                                                           aGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGagat                      >  1:1660361/1‑98 (MQ=255)                                                                           aGGATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGaaaa                      >  2:484901/1‑101 (MQ=255)                                                                            ggATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAAc                     >  2:1334930/1‑101 (MQ=255)                                                                            ggATTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAAc                     <  1:787793/101‑1 (MQ=255)                                                                              aTTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGt                   >  1:366230/1‑101 (MQ=255)                                                                              aTTTTGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGt                   >  1:1402105/1‑101 (MQ=255)                                                                                  tGTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCa               <  1:1030434/101‑1 (MQ=255)                                                                                   gTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAg              <  1:1046349/101‑1 (MQ=255)                                                                                   gTAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAg              <  1:1548287/101‑1 (MQ=255)                                                                                    tAGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCCgg             >  2:852141/1‑101 (MQ=255)                                                                                     aGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGa            <  1:909205/101‑1 (MQ=255)                                                                                     aGTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGa            >  1:566885/1‑101 (MQ=255)                                                                                      gTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCATCCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGAt           <  1:264783/101‑1 (MQ=255)                                                                                      gTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGAt           >  2:1385999/1‑101 (MQ=255)                                                                                      gTTATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGAt           >  2:1681626/1‑101 (MQ=255)                                                                                        tATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGAttt         >  1:1027714/1‑101 (MQ=255)                                                                                        tATCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGAttt         >  2:1306227/1‑101 (MQ=255)                                                                                         aTCTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGATTtg        >  2:1131031/1‑101 (MQ=255)                                                                                           cTAATTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGATTtgtg      <  1:1334936/101‑1 (MQ=255)                                                                                             aaTTCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGATTTGTggg    <  1:1386005/101‑1 (MQ=255)                                                                                               ttCAAC‑TTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGATTTGTGGGAt  >  2:835216/1‑101 (MQ=255)                                                                                                     |                                                                                                CCACACCATAACATAATGGTCTTTGTGTATTGATAATGTTATTCAAAATTGTTTCACTCCAGTCTTCTGGAATAAGGATTTTGTAGTTATCTAATTCAACTTTTTTTATTATAAAATCTTTATCTTCTTTGGTCAACCTACCAGTTACTGCCGCAGCCAACACCGCCTGACGAAAACGTTTCAGGATTTGTGGGAT  >  20/7228‑7423

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