Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	36,231	T→C	intergenic (‑69/+40)	caiD ← / ← caiC	carnitinyl‑CoA dehydratase/putative crotonobetaine/carnitine‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	36,231	0	T	C	98.0%	352.5 / ‑4.5	99	intergenic (‑69/+40)	caiD/caiC	carnitinyl‑CoA dehydratase/putative crotonobetaine/carnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (42/55);  total (43/56)
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

TTGCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCT  >  NC_000913/36134‑36327                                                                                                  |                                                                                                 ttGCGGGTCAGATGTAATGATTCACTCATTCCTTTGCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAgtgt                                                                                                <  2:472315/100‑1 (MQ=255)   gCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGa                                                                                              >  1:1949513/1‑100 (MQ=255)   gCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGa                                                                                              <  2:1498062/100‑1 (MQ=255)   gCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGa                                                                                              >  1:1129257/1‑100 (MQ=255)    cGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGat                                                                                             >  1:1654031/1‑100 (MQ=255)    cGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGat                                                                                             >  2:2036226/1‑100 (MQ=255)    cGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGat                                                                                             <  1:1212446/100‑1 (MQ=255)    cGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGat                                                                                             <  1:596487/100‑1 (MQ=255)     gggTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGatt                                                                                            >  1:17175/1‑100 (MQ=255)      ggTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGatta                                                                                           <  2:897770/100‑1 (MQ=255)      ggTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACAGCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGatta                                                                                           <  2:1098429/100‑1 (MQ=255)        tCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGattatt                                                                                         <  2:1571486/100‑1 (MQ=255)         cAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACAACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGattattn                                                                                        <  2:1459088/100‑2 (MQ=255)          aGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAg                                                                                       <  2:492821/100‑1 (MQ=255)           gATGTAATGATTCACTCATTCCTTTTCTCCATTTTTTTTTTTTTTTTTACGACAACATCCCTGCAAAAAATTCATATTGTTTTAGAGCGTGATTATTAGc                                                                                      >  2:1016053/1‑100 (MQ=255)            aTGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTGTTTTATGTACGACAAAAAAAAAAAAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCt                                                                                     <  2:1569680/100‑1 (MQ=255)              gTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCGTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCttg                                                                                   <  1:1381438/100‑3 (MQ=255)               tAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGc                                                                                  >  1:251531/1‑100 (MQ=255)                 aTGATTTACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAg                                                                                <  1:1849637/100‑1 (MQ=255)                  tGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAgg                                                                               <  1:1241203/100‑1 (MQ=255)                   gATTCACTCATTCCTTTTCTCCATTTTTTTTTTTTTTTTGACGACAACATCCCTGCAAAAAATTCATATTGTTTTAGAGCGTGATTATTAGCTGTCAg                                                                                >  2:1239074/1‑98 (MQ=255)                   gATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAggg                                                                              <  1:562707/100‑1 (MQ=255)                    aTTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGt                                                                             <  1:2036244/100‑1 (MQ=255)                     ttCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTa                                                                            <  2:1949495/100‑1 (MQ=255)                     ttCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTa                                                                            <  2:1860075/100‑1 (MQ=255)                      tCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAg                                                                           >  1:122353/1‑100 (MQ=255)                       cACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTTTTTTAGCCCGTGATTATTATCTGGCAGGGTAGt                                                                          >  2:1702651/1‑100 (MQ=255)                       cACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGt                                                                          <  2:876596/100‑1 (MQ=255)                       cACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGt                                                                          >  2:1732403/1‑100 (MQ=255)                       cACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGAGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGt                                                                          <  2:1078231/100‑1 (MQ=255)                           cATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTccc                                                                      <  2:945602/100‑1 (MQ=255)                             ttCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCctg                                                                    <  2:1592792/100‑1 (MQ=255)                              tCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCctgc                                                                   <  1:1934554/100‑1 (MQ=255)                              tCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCctgc                                                                   <  2:159537/100‑1 (MQ=255)                               ccTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCctgct                                                                  >  1:1965982/1‑100 (MQ=255)                                cTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCctgctg                                                                 >  1:1604452/1‑100 (MQ=255)                                cTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCctgctg                                                                 <  2:395222/100‑1 (MQ=255)                                cTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGAACTACCctgcca                                                                 <  1:1301715/100‑3 (MQ=255)                                cTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTAGCAGGGAACTGCCctgcca                                                                 <  1:1577071/100‑3 (MQ=255)                                   ttCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGttt                                                              <  1:1011494/100‑1 (MQ=255)                                     ctcCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCa                                                            >  1:665198/1‑100 (MQ=255)                                     ctcCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCa                                                            <  1:527552/100‑1 (MQ=255)                                     ctcCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCa                                                            >  2:1879520/1‑100 (MQ=255)                                      tcCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCat                                                           >  1:1428535/1‑100 (MQ=255)                                       ccATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCatt                                                          >  1:979313/1‑100 (MQ=255)                                       ccATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCatt                                                          <  1:1239080/100‑1 (MQ=255)                                        cATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCattt                                                         >  1:1049786/1‑100 (MQ=255)                                          tttttGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCatttat                                                       <  1:1886528/100‑1 (MQ=255)                                           ttttGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCatttatt                                                      <  1:413391/100‑1 (MQ=255)                                           ttttGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCatttatt                                                      >  1:1230336/1‑100 (MQ=255)                                            tttGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCatttattt                                                     <  2:122352/100‑1 (MQ=255)                                            tttGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCatttattt                                                     >  1:1876141/1‑100 (MQ=255)                                             ttGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGTATATTGTTTTAGAGCGTGACTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTc                                                    <  2:1118063/100‑1 (MQ=255)                                             ttGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTc                                                    <  1:368347/100‑1 (MQ=255)                                               gCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAg                                                  >  1:680031/1‑100 (MQ=255)                                               gCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAg                                                  >  2:185857/1‑100 (MQ=255)                                               gCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAg                                                  <  2:1654018/100‑1 (MQ=255)                                                cTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGa                                                 >  1:1595030/1‑100 (MQ=255)                                                 tttttCAGGGACGACAAAATAACTGAAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGa                                                 <  2:17175/99‑1 (MQ=255)                                                  ttttCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAtt                                               >  2:607960/1‑100 (MQ=255)                                                   tttCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAttc                                              <  1:1059191/100‑1 (MQ=255)                                                     tCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAttctt                                            >  2:1419710/1‑100 (MQ=255)                                                      cAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAttcttt                                           <  1:1879537/100‑1 (MQ=255)                                                        gggACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAttctttct                                         <  1:232692/100‑1 (MQ=255)                                                        gggACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAttctttct                                         <  2:502879/100‑1 (MQ=255)                                                          gacgacAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTaa                                       <  1:1732416/100‑1 (MQ=255)                                                          gacgacAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTaa                                       >  1:471838/1‑100 (MQ=255)                                                           acgacAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAat                                      >  2:213275/1‑100 (MQ=255)                                                            tgttgtcgtCCCTGAAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAatt                                     <  2:1105063/92‑1 (MQ=255)                                                             gacAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAatta                                    <  1:1635850/100‑1 (MQ=255)                                                              acaacaTCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAattat                                   >  1:373011/1‑100 (MQ=255)                                                              acaacaTCCCTCCAAAAGATACAGATTGTTCTAGAGCGTGACTCTTAGCTGGCAGGGTAGTTACCTGCTGTTTCATGTATTTCAGATTCTTTCTAattat                                   >  2:1945958/1‑100 (MQ=255)                                                                  caTCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTc                               >  2:1368696/1‑100 (MQ=255)                                                                  caTCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTc                               <  1:1011784/100‑1 (MQ=255)                                                                    tCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTccc                             >  1:1052827/1‑100 (MQ=255)                                                                     cccTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTcccc                            <  1:1002757/100‑1 (MQ=255)                                                                     cccTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTcccc                            >  2:313009/1‑100 (MQ=255)                                                                      ccTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCTCCg                           >  2:2000960/1‑100 (MQ=255)                                                                       cTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCCAATTATTTTCCCCGa                          <  1:1069430/100‑1 (MQ=255)                                                                         gCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGc                        <  1:1419719/100‑1 (MQ=255)                                                                          cAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCa                       <  2:1118731/100‑1 (MQ=255)                                                                                  gCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTg               <  1:72530/100‑1 (MQ=255)                                                                                  gCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTg               <  2:1129252/100‑1 (MQ=255)                                                                                   cATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTgg              >  1:1746739/1‑100 (MQ=255)                                                                                    atatTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTATCCCCGAGCAATTACGTGGc             <  2:1251008/100‑1 (MQ=255)                                                                                    atatTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCTCCGAGCAATTACGTgg              >  2:964081/1‑99 (MQ=255)                                                                                     tatTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCa            <  1:97606/100‑1 (MQ=255)                                                                                      atTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAg           >  2:2044910/1‑100 (MQ=255)                                                                                      atTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAg           >  2:127440/1‑100 (MQ=255)                                                                                       ttGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGa          <  1:1751030/100‑1 (MQ=255)                                                                                       ttGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGa          >  2:399067/1‑100 (MQ=255)                                                                                       ttGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGa          <  2:1230330/100‑1 (MQ=255)                                                                                        tGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGAt         >  1:277166/1‑100 (MQ=255)                                                                                        tGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGAt         >  2:580919/1‑100 (MQ=255)                                                                                          ttttAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCt       <  1:54432/100‑1 (MQ=255)                                                                                          ttttAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAA‑TACGTGGCAGATCtt      >  1:1642736/1‑100 (MQ=255)                                                                                           tttAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCtt      <  1:575451/100‑1 (MQ=255)                                                                                            ttAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCttt     <  2:1049781/100‑1 (MQ=255)                                                                                               gagCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCt  >  1:1336877/1‑100 (MQ=255)                                                                                               gagCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCt  >  1:1012273/1‑100 (MQ=255)                                                                                                  |                                                                                                 TTGCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCT  >  NC_000913/36134‑36327

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